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ABSTRACT 


Changes of length of a lead rod due to thermal expansion have been 
measured between 4-2°K and 10°kK by a heterodyne-beat method in which 
the changes modify the capacity of a condenser in a resonant circuit. The 
variation with temperature of the Griineisen constant y was deduced ; the 
accuracy was just insufficient to show up the jump in thermal expansion 
coefficient at the superconducting transition. 


$1. INTRODUCTION 


LittLe work has been done on the measurement of thermal expansion of 
solids at low temperatures in spite of the importance of such measure- 
ments in the study of the equation of state of solids. Most of the existing 
measurements are based on the use of optical interference, but because of 
the difficulty of using this principle at low temperatures and the small 
magnitude of the thermal expansion at low temperatures, large 
temperature intervals have usually been necessary to produce measurable 
changes of length. Thus expansion coefficients obtained in this way are 
averages and cannot exhibit the true temperature variation of this 
property. To obtain significant results at low temperatures it is 
necessary to use a technique capable of higher sensitivity than that of 
most of the previous methods so that relatively smaller temperature 
intervals can be used. 

This paper describes a method in which the change in the length of the 
specimen caused by heating modifies the capacity of a parallel plate 
condenser and this in turn produces a measurable change in the resonant 
frequency of a high frequency circuit from which the change of length can 
be deduced. A similar method has been used by Bijl and Pullan (1954, 1955) 
for measuring thermal expansion of copper and aluminium. The present 
method has been used for measuring the thermal expansion of lead in the 
range 4-10°k. Lead was chosen because it has a low Debye temperature 
and is therefore expected to have a fairly high coefficient of expansion in 
this range. Also it is superconducting below 7-22°K and it was hoped 
that the measurements might show the small discontinuity in the 
expansion coefficient at the transition and the small volume changes when 
the transition occurs in a magnetic field, as predicted by thermodynamics. 
Unfortunately, the sensitivity did not prove adequate to show up the 
former effect and technical difficulties prevented any reliable measure- 
ments of the latter effect. 
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A brief report of this work has already been published (Dheer and 
Surange 1957). 
§ 2. EXPERIMENTAL DETAILS 


2.1. The Cryostat Assembly 


The cryostat used in the measurements is shown in fig. 1. The fixed 
plate P, of the parallel plate condenser is rigidly attached to the top of a 
hermetically sealed brass can A through a glass metal seal G, which 
insulates it electrically from A. The specimen §, a cylinder 4-87 cm long 
and 0:7em diameter of spectroscopically pure lead (Johnson Matthey 
Ltd.) carries the other (earthed) condenser plate P,, and is mounted on a 


Fig. 1 


The Cryostat. 


copper base B supported in the can A by three thick tubes W of low 
thermal conductance (cupronickel). The two plates are adjusted to be 
in contact at room temperature and a gap of about 4mm is automatically 
obtained by the differential contraction of the specimen and the brass 
can A when the cryostat is cooled to liquid helium temperature. The 
electrical lead from P, is taken out of the cryostat through a shielding 
tube T, and a glass metal seal at the top of the cryostat. 
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A carbon resistor R (Allen Bradley 120hm, 4} watt) cemented inside a 
hole in the copper base B is used for measuring the temperature of the 
specimen. A search coil is wound on the specimen and serves to reveal 
the appearance of the superconducting state at the transition temperature 
by the abrupt fall in the magnetic permeability of the metal at this 
temperature. The electrical leads to the search coil, thermometer and 
heater are maintained in thermal contact with the helium bath by filling 
the metal capillary T, with beeswax over part of its length; the exit of 
the tube is sealed off with picein. 

On filling the cryostat with liquid helium the air inside the can freezes 
out and the specimen cools in about an hour from its initial temperature 
(S80°K) to 4:2°K by conduction through the tubes W. By passing current 
through the heater H, the specimen can be maintained at a temperature 
above the bath temperature, determined by the power input and the 
effective conductance between the specimen and the helium bath; the 
main contribution to this conductance comes from the tubes W. The 
time for attaining thermal equilibrium is of the order of a minute between 
4°K and 10°K. The power required to raise the temperature from 4°K to 
10°K was about 20mw, which is not excessive from the point of view of 
evaporating liquid helium (an extra 5em3/hour). There would probably 
be little difficulty in raising the temperature to about 50°K but the present 
experiments did not go above 10°K because of the rapidly falling sensitivity 
of the carbon thermometer used. 


2.2. Determination of Temperature 
For the measurement of the temperature of the specimen the resistance 
R of the carbon thermometer was measured by a Wheatstone bridge, and 
was assumed to vary with temperature according to the empirical 
equation of Clement and Quinnell (1952 a) 
loo hen lloe i= A Bie ee ae eet ee (2) 
where K, A and B are constants, which were determined by calibration 
between 2-4°K and 4:2°K against the vapour pressure of liquid helium and 
at 7-22°k by noting the resistance at the superconducting transition 
temperature of lead. It may be mentioned that the superconducting 
transition in the measuring field of about 4 gauss was sharp (width of order 
0-02°K), thus indicating that the lead was of good quality as regards 
superconductivity. Subsequently, it was verified by direct calibration 
against a gas thermometer that the temperatures determined in this way 
did not differ by more than 0-2° from the actual temperatures, in the 
range 4-10°k. The sensitivity of the bridge was such that temperature 
changes of order 10-3 degrees could be detected at about 10°K; at lower 
temperatures the sensitivity was greater. 


2.3. The Oscillator 


The frequency f of the oscillator was 7-5 Mc/sec; the Franklin oscillator 
circuit (see Ladner and Stoner 1950) was used for two reasons. First 
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because it requires a very loose coupling between the resonant circuit and 
the valve maintaining circuit. This reduces the effect of valve capaci- 
tances, anode voltage, etc., on the oscillation frequency which is therefore 
determined almost entirely by the constants of the resonant circuit alone 
so that a high degree of frequency stability can be attained. Secondly, 
one side of the resonant circuit can be earthed and the design of the low 
temperature part of the apparatus is thereby much simplified. Also by 
earthing the plate P, of the condenser any possible radio-frequency 
coupling of the thermometer and heater with the oscillator circuit was 
eliminated. 

To measure the change Af in the frequency of the oscillator associated 
with the expansion of the specimen, the oscillator was allowed to beat 
against a similar oscillator of fixed frequency and the frequency of the 
audible beat note was determined by observing Lissajous figures produced 
by combining the beat note with an audiofrequency note from a calibrated 
oscillator. Changes in f were then read directly on the dial of the audio- 
oscillator. A typical measurement is illustrated in fig. 2 showing the 
variation of the audiofrequency as the specimen is heated from 4-2°K to 
5:54°K. The heating current was switched on at B and the specimen 
attained a steady temperature at P (as indicated by the carbon thermo- 
meter) beyond which the oscillator is seen to follow its natural drift. At 
R the heating current was switched off and the specimen cooled to the 
bath temperature as is shown by the portion CD of the curve. The parts 
AB and CD were joined by a smooth curve and the mean change Af= X Y 
was estimated from the graph. At least three such measurements were 
made for each temperature 7’; each point plotted in fig. 3 is the average 
of such measurements. 

Figure 2 also illustrates the performance of the oscillator; it can be 
seen that the drift is of order 3 c/s per minute which may vary appreciably 
over several minutes, and the scatter of the points is presumably due to 
irregular changes of frequency of order 5c/s. It is estimated that Af can 
be determined to within about 10¢/s. 


2.4. Measurement of Changes of Length Al 

The change of length Al of the specimen is related to Af by the relation 

A 0+, : 
5,7 Cle Fo ee 
where A is the area and C the capacity of the parallel plate condenser and 
Cy the rest of the capacity in the resonant circuit (capacities in E.8.U.). 
(eceGon of Al thus requires a knowledge not only of A and f which 
present no difficulty, but also of C and C,. The capacity C, consists of 
C, (of order 15cm), the capacity of the jead from the plate P, (fig. 1) of 
ite condenser C with respect to the shielding tube T,; C, (of ardor 10 cm) 
the capacity of a small calibrated condenser and C, (of order 18 cm) the 
self capacity of the coil which serves as the md coranee of the resonant 
circuit together with the capacities introduced by the connecting wires. 
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The condenser C, is introduced deliberately to permit easy adjustment of 
the resonant frequency and for calibrating purposes. The various 
unknown capacities are measured as follows. 

The capacity C+ C, is determined by disconnecting the cryostat from 
the oscillator circuit and then measuring the amount by which the 
capacity of the variable condenser ©, has to be increased in order to bring 
the oscillator frequency to its original value. Now without disturbing 
the condenser lead the can A is opened and assembled again after 
removing the specimen 8 and the plate P, so that C is removed from the 
resonant circuit. C, is now determined in the same way with the help of 


Fig. 2 
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A typical measurement of Af. The temperature at P is 5-54°K ; Af=47 c/s 
which corresponds to Al//=0-3 x 10-®. 


the variable condenser, and the difference of the two determinations 
gives O. The total capacity C+C, is determined by introducing a small 
change 80 (about 0-5cm) in it by varying C,, and measuring the corres- 
ponding change $f in the oscillator frequency ; it is easily seen that 


Cee Ces (OC Opa a. se ee oe (3) 
Combining (2) and (3) we get 
Al = (A/42C?)(5C/8f) AF. . 2 ae 


Substitution of the values of A(5-06cm?), C(15cm) and 8C/8f(1-6 x 10~° 
emsec) gives Al/l=6-5x10-°Af. As already mentioned, Af can be 
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measured to within about 10c/s so that the accuracy of our measure- 
ments of Al/l is about 0:07 x 10-6. This is comparable to that attained 
by Rubin et al. (1954) in their optical interference method. 

The error introduced in the expansion measurements due to thermal 
strains in the supporting tubes W, produced by heating, was ascertained 
by a special experiment in which the can A was made much shorter, S was 
removed and P, was supported very close to P, which was attached directly 
to B. The total change in length produced by heating the base B up to 
10°K was found to be 3 x 10-7cm which is less than 1°% of the expansion 
of the specimen from 4:2 to 10°K. 


§ 3. RESULTS 
The temperature variation of (1,,—1,.99)/l is shown in fig. 3. For com- 
parison the corresponding readings of McLennan et al. (1931) obtained by 
an optical lever method are also shown. It can be seen that their points 


Fig. 3 


(12145)/ 0 x 10cm 


10°K 
Temperature variation of (1p— I4.09)/U. 


© Our observations. 
® McLennan et al. 
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show a greater scatter and lie rather above ours; the latter fact may 
perhaps be due to a negative expansion coefficient of the quartz supports 
used in their apparatus (see Bijl and Pullan 1955, p. 290). 

An attempt was made to find out if the results show any discontinuity 
in the coefficient of linear thermal expansion £, at the superconducting 
transition temperature 7-22°K. Two cubic equations of the form 


(lp —Iy.99) [b= A(T — 7-22) + B(T — 7-22)2 + O(T — 7-22)8 


were fitted separately to the data for 7<7-22°K and T>7-22°K; 
(lp —ly.22)/1 was determined by reading (l,.9)—J499)/1 from fig. 3 and sub- 
tracting it from the measured value (1,—J459)/l. The difference in the 
two values of the coefficients A thus obtained gives the value of any 
possible discontinuity £,,—f8, at 7-22°K. This procedure suggested a 
value 2-6 x 10~-’deg?. Since, however, the value of (1,95 —Jy.99)/ deter- 
mined from the graph is subject to experimental error, it was desirable to 
check the reliability of the procedure and the above analysis was repeated 
for each of two extreme possible estimates of the value of (1,95 —14.99)/1; 
this led to estimates of —2 x 10-7 and 7 x 10-7 deg~!.__ It is clear then that 
the data are not precise enough to establish definitely the existence of a 
discontinuity in 8 at 7’, though they are evidently not inconsistent with 
the discontinuity 8, —8,=0-66x10-*deg predicted by the thermo- 


dynamic formula 
Nog opel 0H, 
Bn—-Bs= a ak ap ie 


In the numerical evaluation of this formula we have assumed 


(0H,/0T)» = 223 gauss deg 
and 
(0H,/0p) p= — 11-1 x 10° gauss dyne-1 cm? 


(Olsen and Rohrer 1957). Recently we have learnt that Lasarev and 
Sudovstov have reported measurements of a discontinuity in 6 to the 
Moscow Low Temperature conference (July 1957). 

An attempt was also made to measure the changes of length when 
superconductivity is destroyed by a magnetic field. This could not be 
done directly because the frequency of the oscillator proved to be appre- 
ciably modified by the presence of a magnetic field, perhaps due to dimen- 
sional changes of magnetostrictive origin. To eliminate these effects the 
following indirect procedure was tried. First the change AJ, in the length 
of the specimen produced by heating it from a temperature 7’, below 7’, 
to a temperature 7’, just above 7’, was measured. The magnetic field was 
then switched on and the corresponding change A/, was measured when 
the specimen is heated from 7, to T,; (Al,—Al,)/I gives the relative 
change of length (/,—1,,)/l associated with the transition at 7. A few 
measurements made between 3 and 7°K agreed qualitatively with those of 
Olsen and Rohrer (1957) but showed a considerable scatter (up to a factor 2 
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in Aljl at any one temperature). Further investigation suggested that 
the results could not be regarded as reliable because different values of 
Al, and Al, were obtained even for 7',>7',; this indicated that the dimen- 
sional changes produced by the magnetic field do not remain constant as 
the temperature is varied. Evidently then the apparatus in its present 
form is not very suitable for magnetic measurements. 

It is of interest to compare the results on thermal expansion with the 
Griineisen relation. At low temperatures the electronic contribution to 
thermal expansion might be expected to be appreciable and the 
Griineisen’s formula can be expressed as (see, for instance, Bij] and Pullan 
1955) 

SB/V= VG ee, 2 ee ek ee eee 


where y is the compressibility and C, and C, the lattice and electronic 
specific heats per unit volume respectively. In order to avoid compli- 
cations due to the superconducting state, the measured values of 
Ln—lyoos Were converted to 1p, —ly29, by using Olsen and Rohrer’s (1957) 
data on 1/,—1,; the subscripts s and n denote the superconducting and 
normal state respectively. The correction was 2% at 10°kK and about 
10% below 7°x. The method of least squares was then used to fit an 
equation of the form 


(Lin —la.oon)/L= al T — 4-20) +0(T — 4-20)? 4+ o(7'—4-20)8 . . (6) 


to the data; the values a=16-2x10-%, b=4:2x10-% and c=2-4x10- 
were found. The observed values departed from the relation (6) by at 
most 2x 10~7 but were mostly within 0-7x 10~*. Differentation of eqn. 
(6) gives 

Bx 108=108-5—52-374+727T2, . . . . . (7) 


This formula is of course valid only as an interpolation formula for the 
actual data and cannot be used reliably for extrapolation far outside the 
range 5-9°K. In order to assess the reliability of the method of deter- 
mining 6, a formula similar to (6) but containing an extra term in 
(7'— 4-20)4 was also fitted to the data. Differentation gave 


B x 108§= — 82-34 41-37 —7-672+0-76T3 . . . . (8) 


which departs from (7) by at most about 3% except at the ends of the 
range (5°K and 10°k) where the departure reaches 7%. Thus (7) which 
is used in the discussion below is unlikely to be in error by more than a 
few per cent in the relevant range. 

The second term on the right-hand side of eqn. (5) was estimated using 
the value y,=1-7 determined by Olsen and Rohrer (1957) and the value 
C= 4:20 x 10-° T cal deg-tcm~3 as determined by Horowitz et al. (1952). 
We then have y,C,=0-71 x 10-4 T' cal deg-t cm-3;, this is about 1% of 0, 


at 10°K and becomes 5% at 4:5°K. OC, was determined by combining the 
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measurements of Horowitz et al. in the liquid helium and liquid hydro- 
gen range of temperatures and those of Clement and Quinnell (1952 b) near 
7:22°K. Since there seem to be no reliable measurements of specific heat 
other than those mentioned it is rather difficult to assess the reliability of 
the estimate of C|. This uncertainty is likely to be greatest (~5%) at 
10°K because of the large range of interpolation involved (7° to 14°K). 
Values of y, calculated by taking the room temperature value of y are 
plotted against 7'/@ in fig. 4; @=96-3 is the Debye characteristic 
temperature oflead. Owing to the various uncertainties the values of y, are 
probably not reliable to better than 0-1 near 10°K and are quite uncertain 


Fig. 4 


T/@ 


Variation of y, with 7'/@. The room temperature value of y=2-64 is shown 
by a line on the axis of ordinates. 


below 5°K where both f and C, are becoming small and unreliable. It is 
satisfactory that y, seems to flatten off reasonably well towards the high 
temperature value of 2-64 but there is no evidence for a flattening off at 
very low temperatures such as might be expected in the region cre the 
specific heat begins to follow a 7° law (this is probably below 4:2°K). The 
fall off in y, from its high temperature value is similar for lead to that 
found by Bijl and Pullan (1955) for copper and aluminium but occurs at 
an appreciably lower value of 7'/@. 
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ABSTRACT 


Linear discontinuities strictly in (1, 1, 1) planes have been observed on the 
natural octahedral faces of a diamond. The lines are studied optically. 
The characteristics of the lines suggest that they originate in slip which could 
well have occurred while the diamond was still growing and presumably hot. 
The slip hypothesis is the only reasonable one which accounts for the observed 
characteristics of the lines. 


$1. InTRODUCTION 
Diamonps are known to give good cleavages in octahedral planes. It 
is, however, not possible to induce slip in diamond by application of 
shearing stresses, at least not at room temperature. 

Not much is known about the behaviour of diamond at high tempera- 
tures. Silicon and germanium crystals which have the diamond 
structure are also brittle at room temperatures. At higher temperatures, 
however, both silicon and germanium become ductile. Slip bands have 
been observed on germanium crystals deformed at temperatures about 
500°o and for silicon about 900°c (Gallagher 1952, Maddin and Chen 
1954). These temperatures are about two-thirds of the way to the 
melting points of the corresponding crystals. Diamond therefore might 
possibly be expected to slip at high temperatures, which in this case 
might well be above 2000°c. Experiments with diamonds at such high 
temperatures would be formidably difficult, and in any case diamond is 
not stable at ordinary pressures at these temperatures (Robertson et al. 
1934, 1936). 

In the present paper observations are reported on some growth 
features on diamond faces which seem to imply that the diamond con- 
cerned might well have been deformed plastically, presumably then at 
a high temperature. 


§ 2. OBSERVATIONS 
The type of linear discontinuity observed here appears to be 
relatively rare in diamond. Nevertheless we have observed it on six 
diamonds out of many hundreds investigated. One diamond showed 
the linear discontinuities under consideration arranged in such a way 
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that it made their nature fairly clear. Our work has been concentrated 
mainly on this specimen which was investigated in detail. This stone 
was an almost perfect octahedron weighing 1-2 carats. Some 50 line 
markings were observed on 6 of the 8 faces of the octahedron (the other 
2 faces did not show any). All the lines were oriented strictly parallel 
to octahedral planes of the crystal. Moreover, many of the lines could 
be traced continuously from one face to the neighbouring faces in such 
a way as to define a plane which in all cases was parallel to octahedral 
planes. A photograph (x30) of one of the faces of the octahedron is 
shown in fig. 1, Pl. 28. The discontinuity lines will be shown later with 
much higher magnifications. Rows of growth trigons are seen and their 
centres are just on the lines, as if strung on them. 

An overall view on the orientation of the lines in general, and of their 
interrelation is given by fig. 2, Pl. 29. This shows a photograph of a 
cardboard model of the crystal made by combining the photomicro- 
graphs of all the 8 faces of the octahedron. The octahedral planes 
defined by the lines continue across the crystal and can clearly be fixed. 
as seen on the photograph. 

Figure 3, Pl. 29 shows part of a line using higher magnification 
(x 600). It shows clearly how smaller trigons are centred on the lines. 
The dark bold curved lines below the row of trigons are edges of growth 
sheets which emerged from the lower part of the crystal, below the 
picture spreading upwards, and stopped at the row of trigons. Figure 4, 
Pl. 30 shows (x 600) triple limes with only about 5 microns separating 
the individual lines. Again the smaller trigons are centred on one or 
another of the three lines. 

The topography of the surface in the neighbourhood of the lines is 
revealed by fig. 5, Pl. 30 and fig 6, Pl. 31. Figure 5 is a multiple-beam 
interferogram (x 250) taken across one of the lines. Clearly they are 
regions of discontinuity. On approaching the row of trigons from the 
lower part of the picture there is a sharp discontinuous drop from the 
level of the surface. On approaching the line from the upper side, i.e. 
from the direction which the apexes of the row of trigons are directed, 
there is a gradual fall until the line is reached when again an abrupt fall 
occurs. Figure 6 (x 350) shows a high-dispersion multiple-beam inter- 
ferogram. It can again be seen clearly how growth waves, which 
emerged from below, were stopped by a barricade formed by the bases 
of the trigons, while the growth waves on the upper part are seen to 
penetrate between the apexes of the trigons, so that a gradual slope is 
formed. Growth from this direction often crosses the discontinuity 
itself as can be seen in fig. 6 and in figs. 3 and 4. 


2.1. Step Heights 
These were measured by the methods of light-profile microscopy 
(Tolansky 1952) and by multiple-beam interferometry, the latter 
including fringes of equal chromatic order (Tolansky 1948). The results 
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may be summed as follows : (i) The steps across the discontinuity lines 
range in height from nearly zero up to about 2 microns, and differ not 
only from line to line, but also along the same line. In some cases the 
step changes progressively along the line. Usually, however, it changes 
irregularly along a given line. (ii) Without exception, all the steps are 
directed in such a way that the lower level is towards the bases of the 
trigons, while the upper level is towards the apexes of the trigons. 

In the following discussion it will be shown how the characteristics of 
the observed discontinuous lines can be accounted for by assuming that 
the discontinuity lines originate from disturbances in the growth of the 
crystal as a result of slip, which we postulate could have occurred while 
the crystal was still in the process of growth, and presumably at a fairly 
high temperature. 


$3. Discussion 


The high concentration of growth trigons along the lines implies that 
the observed discontinuities on the surfaces of the crystal were created 
whilst it was still growing. 

Misfits bound to octahedral planes of the crystal can be one of the 
following: (a) cracks in the crystal, (b) twinnings, (c) stacking fault 
surfaces, (d) slip in the crystal. We do of course reject outright the 
unfounded suggestion of Omar and Kenawi (1957) that the growth 
trigons might be etch figures. We consider (d) is the most probable 
explanation here for the following reasons. 

(a) Cracks: It is most improbable that the crystal could have 
cracked in so many places, in some cases into slices only a few microns 
thick (see fig. 4), and then subsequently healed itself to such a perfection, 
for the crystal as a whole is a good crystal. The possibility of cracks 
can be discarded in this case. 

(6) Twinning : It is known that the orientation of trigons changes 
on crossing a twinning boundary. As this does not happen when 
crossing a discontinuity line on our crystal, it can be concluded that the 
parts of the crystal on both sides of the discontinuity line are not 
twinned relative to each other. We have further proved this point by 
etching experiments. We have already firmly established in this 
laboratory that a diamond etched lightly in hot potassium nitrate 
produces small triangular etch pits oriented opposite to growth trigons. 
Such etch pits reverse over a twin boundary. Figure 7, Pl. 31 shows 
the etch pits (x 400) over one of our discontinuity lines. Etching is 
preferentially stronger on the line itself, but the main issue clearly seen 
is the identical orientation of the pits on either side. There remains to 
be considered the possibility that submicroscopic thin twinned slices 
divide both parts of the crystal. From our microscopic observations it 
is certain that the thickness of any such postulated twinned regions, if 
present, is always less than 0-2 microns. It seems most unreasonable 
to assume that the crystal contains so many twinned regions that for 
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some strange reason are all limited to submicroscopic thicknesses. The 
only reasonable limit would then be that all of them are monolayer 
twins, but this is effectively identical with stacking fault surfaces, which 
are to be discussed now. 

(c) and (d) Stacking fault surfaces and slip: It is more difficult to 
decide between these two possibilities for slip with a slip vector which 
corresponds to a translation from perfect lattice points to twin-lattice 
points will leave the crystal with a stacking fault surface. The only 
difference between (c) and (d) in this case would be in the process of 
formation of the stacking fault surface. Under (c) we understand the 
case in which the stacking fault surface was formed during the process 
of growth on the surface of the growing crystal ; while under (d) we 
consider slip of the growing crystal by the action of shearing forces 
parallel to its octahedral planes. 

On assumption (c) the distribution of the stacking fault surfaces in 
the crystal (some of them pass the crystal near its centre) would imply 
that some of them were formed in very early stages of growth of the 
crystal. It would then be very difficult to explain how these stacking 
fault surfaces extended in such a perfect way right out to the external 
faces of the crystal. This, further, seems much more improbable when 
taking into account the excellent crystallographic form of the octahedron 
shape of the crystal. 

Only the assumption of the occurrence of slip remains, and it will be 
shown now how this assumption accounts for the observed characteristics 
of the lines in a reasonably satisfactory manner. 

We assume that slip occurred when the crystal was still growing, and 
that the observed discontinuities are the result of overgrowth on the 
surfaces of the already slipped crystal. Nothing can therefore be said 
as to the original number of slip lines nor about the original amount 
of slip. On a growing crystal surface, perfect slip in integral lattice 
multiples would supply perfect steps for rapid advance of growth layers. 
Slip would soon be overgrown leaving no traces on the surface of the 
original slip lines. 

However, an imperfect step due to slip, which is not an integral of the 
lattice, would form a misfit boundary on the surface, at which growth 
emerging from either side would be retarded. Trigons should then be 
expected to be formed along the slip lines where the lattice was disturbed 
(Halperin 1954). 

Once the row of trigons has started to form, an asymetry is introduced 
with regard to the advance of the oncoming growth layers (from both 
sides) towards the slip line. On one side, where the bases of the trigons 
are arranged parallel to the slip line, they form a barricade at which 
growth is necessarily stopped. On the other side, however, growth 
layers are free to penetrate between the apexes of the trigons and stop 
just on the slip line. Thus, no matter what the direction and height of 
the slip step originally, the side towards the apexes of the trigons will 
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grow much faster and will finally become the higher level of the step, 
indeed exactly as we find experimentally. Growth, although retarded in 
rate, might continue across the step, and the final step height at a given 
point on the slip line will then depend on the degree of perfection of 
the step at that point as well as on other factors, which will determine 
the rate of penetration of growth across the slip line. This could account 
for the irregular changes in the step heights found along a given line. 
A striking example (x 375) is shown in fig. 8, Pl. 32 where a part of the 
upper slip line is completely covered by growth which penetrated across 
the slip lines, for the line is partially obliterated over some half of the 
picture, reappearing at the sides. 

It follows from these considerations that the rows of trigons formed 
along the slip lines actually preserve them by preventing growth on one 
side from reaching the lines. This is also demonstrated (x 160) in fig. 9, 
Pl. 32. Here, besides the four slip lines with trigons along them, a trace 
of another appears in the centre. The impression is given of a line 
obliterated everywhere except where there is a trigon, and here a small 
part of it has been preserved in the neighbourhood of this trigon. Atten- 
tion is incidently drawn to the fact that lines of the character analysed 
here have already long been reported in the literature. Sutton (1928) 
mentions observations of rows of trigons ‘ strung’ on straight lines and 
gives a description of such a line in one of his sketches. Williams (1932) 
also shows microphotographs which also demonstrate this phenomenon. 

It seems then that a slip hypothesis is the only one which accounts for 
the characteristics of the observed lines. This is also in accord with 
theoretical considerations (Jaswon and Dove 1955) by which potentially 
glissile dislocations exist in diamond, but are rendered immobile at low 
temperatures. 

We cannot of course even make a guess as to the temperature at which 
the conjectured slip would have taken place. 
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ABSTRACT 


The ionization produced by 27 electron—positron pairs with energies 
between 80 and 2000 Bev, recorded in stacks of emulsion, has been studied. 
The energies were estimated from the development of the soft cascades in 
an auxiliary experiment. Suppression of the ionization near the origin of 
the pairs was observed, and the measured effect found to be in satisfactory 
agreement with theoretical predictions. The observations allow an inde- 
pendent estimate of the energy of the individual pair from the initial 
ionization, and for pairs with energy ~300 Bey most of the values thus 
obtained were in good accord with those from other methods. 

There is a correlation between the suppression effect and the disparity in 
energy of the two particles of a pair. This suggests that ionization measure- 
ments give information similar to that provided by studies of the relative 
scattering of the two components of a pair. 


$1. INTRODUCTION 


It was suggested by King (private communication, see Perkins 1955), 
and later confirmed by Perkins (1955) that the ionization of electron— 
positron pairs of extremely high energy is reduced near their origin, owing 
to the destructive interference of the electromagnetic fields of two particles 
with equal charges opposite in sign. 

Roughly speaking, a pair with energy > 100 Bev is created with an 
opening angle <2.x10-°radians (Borsellino 1953), and travels a few 
hundred microns before the separation between the electron and positron 
exceeds the maximum effective impact parameter for the ionization of 
atoms (10-®-10~*em). For smaller values of the separation, the ionization 
is smaller than that for two independent, singly charged, particles, each of 
minimum ionization. 

Although the detailed features of the suppression effect may depend on 
the actual mechanisms of ionization, and grain-formation in emulsion, an 
application of the simple theory of energy loss is known to give a fairly 
satisfactory account of this phenomenon. In addition to its theoretical 
implications, the phenomenon has some practical importance, for it 
provides a useful method of determining the energy of pairs of extremely 
great energy—up to 10°-10%ey, 
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There are several established methods of determining the energy of 
electron pairs. At low energies (<10 Bev), the measurements either of 
the multiple scattering, or of opening angle have proved useful. For much 
higher energies, however (>100Bev), it is difficult to obtain reliable 
information by such measurements. For this energy region, a method of 
energy determination has been developed by Pinkau, based on measuring 
the lateral spread of the ensuing soft cascade (Pinkau 1956, 1957). This 
method has the advantage that several observations can be made at 
different depths in the cascade, so that errors due to statistical fluctuations 
are reduced. The disadvantages of the method are that it is laborious; 
and also that the development of the cascade has to be observed over 
several radiation units. This imposes certain geometrical restrictions 
which are important when the original pairs convert near the edge of an 
emulsion stack. On the other hand, the method based on the decrease of 
ionization near the origin of a pair requires only simple measurements over 
the first few millimetres of the track. 

At the present time, one of the great difficulties associated with the 
investigations of interactions of extremely high energy is the rarity of 
suitable events: it is therefore very important to extract as much informa- 
tion as possible from the limited material available. Because of the 
importance of the investigations, it seemed desirable to study the ionizations 
near the origin of pairs in more detail, and to establish a practicable method 
of energy determination based on this phenomenon. 

In the present paper, we shall show that the theory of the effect is in 
satisfactory agreement with observation, and that ionization measure- 
ments near the origin of pairs can be used to make reliable energy estimates 
on pairs with energy > 300 Bev. 


§ 2. THEORY 
The effect has been studied theoretically by several authors. The 
calculation of the ionization produced by a pair has been performed by 
Cudakov (1955), using a classical treatment. Yekutieli (1957), and 
Mito and Ezawa (1957), made improvements in the calculations by using 
a semi-classical treatment, whilst Burkhardt (1957) attempted to eluci- 
date the quantum theoretical aspects of the phenomenon: Burkhardt 
calculated the ionization produced by a wave-packet representing the 
motion of the two particles moving with a small angle of divergence. His 
wave-packet has a spread corresponding to the distribution of the opening 
angle of the pair produced by a photon of a definite energy. He showed, 
that a quantum-mechanical calculation of this type, gives the same 
results as the simple classical treatments. In the following, we shall use 
the results given by Mito and Ezawa, because it predicts a reasonable 
asymptotic value of the ionization for large separations of the particles. 
Following Mito and Ezawa, the ionization produced by a pair, Ypair 
with a separation ) (in cm) may be written 


dH 47Ne*| . ‘ b F 
Jpair — T= f [ 2@-24,(»0<) |} --« « (2.1) 


x mc” 
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whsts Vp = |= | 7:36 1018 seg7t oe ene) 
M 
is the plasma frequency of the emulsion, and K,(x) the modified Bessel 
function with imaginary argument, which has the asymptotic form, 
K,(«)——>0.  @ is a parameter which represents the magnitude of the 
a—>() 


ionization produced by a single relativistic particle, (go). It is defined by 


the relation 
Q = flog (4c?/y?v,,27,2), go=(4rNet/m.c?7)Q, . . . (2.3) 


where y= 1-781 (Euler constant). 7, represents the minimum value of the 
impact parameter, and is left as a free parameter. The ratio of the ion- 
ization produced by a pair, to that of a single charge relativistic particle, is 


then given by: 
2x Sostt —2[ 1-7 Ko( 2) | Ar ert oe 
Jo Q ¢ 


To evaluate the magnitude of the suppression effect, we must know the 
separation, b, between the two components of the pair, which is deter- 
mined by two principal factors: (i) the initial angle of divergence of the 
pair at creation: and, (ii) the subsequent multiple scattering of the 
electron and positron. 

The effect of multiple scattermg has been evaluated by many authors 
(Scott 1949). We shall consider two particles of a pair, with energies 1, 
and H_ (total energy H,= H_, + H_) respectively, travelling initially parallel 
from a common origin. Then the distribution of the spatial separation 6’ 
between them will be 


12 
P(b') db’ dé’ = exp ( 2 =) b'db'dd’ (d’: azimuthal angle) (2.5) 


at distance d from the origin. Here o is the root mean square separation : 


o=4/{b'2) = J3 Jats 73) uel (20) 


and 4 is defined in terms of the conventional scattering constant K for 
angles between successive chords: 


37 
A= = EE isis Aine Pate ane Mee ALS 
We adopt for K a value corresponding to small cell sizes, 


K =4-71 Mev cm—12, Se es me eS 


On the other hand, the most Se initial opening angle for the pair is 
estimated by Borsellino (1953), a 


4m ,c? , (EH 
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The corresponding separation b”, at a distance d from the origin, due only 
to the initial opening angle, is then given by 


rn 4nge* E.. 
i = Te | as(=) | Se ee ecu 


Here ¢(z) is a slowly varying function of x with a minimum ¢(4)=1. 
Comparing the dependence of (2.6) and (2.10) on d, we see that the 

effect of multiple scattering predominates after a certain distance, d, 

given by 


7 (m.c?)* 
d~ one 


for the case of the equi-partition of energy (H,=H_=4H,). If the energy 
is not equally distributed between the two components, this distance will 
be reduced by a factor 

: 84°(H./H)) 


CELTS VS 


~ 150 p, fee ee (21) 


(2.12) 


Therefore, in the cases of practical interest, where the suppression is 
appreciable over distances > 200, both effects contribute almost equal 
amounts to the relative separation of the pair. 

To reduce the expression (2.4) to a more practical form, we note that in 
the region of appreciable suppression (K, > 0-7), we can approximate (2-4) 


es ia ae: capes as = 1+ Glog218% 10% | (bin cm). (2.13) 
Q © 2¢ Q 


Then, substituting the expression 


A=) 2b! 2b"b cos. se 2204) 
for b, and averaging over 0’ and ¢’ with the distribution (2.5) we obtain the 
formula 1 1 pi 

~ = 2- Op a |e 2.15 
R | 1+ Glee 18 x 10° 6 30 i( =) | ( ) 


where Hi(x) is the exponential integral. The second term in the bracket 
is due to the initial opening angle, and the third term represents the 
additional effect of scattering. In the limit of small and large scattering, 
(2.15) reduces to: 


ine E ioe vey? | - E log 218 x 108 u | . (2.16) 


0) 2c 
/ 
ane R~| 14 glow | = E ple 1-63 x 10% |, 5 Wee) 


respectively. It can be seen in the last formula that the effective separa- 
tion due to the scattering is oy~!? instead of o. 

So far we have assumed that we know the individual energies of the two 
electrons of a pair. In most cases, however, we know only the total 
energy of both electrons as estimated by cascade measurements. We 
shall therefore next evaluate the average ionization of pairs with a given 
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total energy, taking into account the distribution in values of the dis- 
parity. For the distribution we shall use that resulting from the assumption 
of complete screening given by a conventional theory. The average 
value of logo will then be 


log o = log (const. (1/H.,? + 1/H_?)") =log (const. +/2(2/H,)), (2.18) 


where H,,, has been calculated by averaging over all values of the dis- 
parity up to a minimum value (#, /H,)=0-05, and is given by 


KR -006H,. oo tag, Oca 


Concerning the term containing 6”, we do not know the precise variation 
of this quantity with the disparity. For the sake of simplicity, we shall 
assume that the new effective value of b” (eqn. (2.10)), can be obtained by 
replacing H, by L,, in the factor outside the bracket, and setting ¢ equal 
to unity. 

Thus, in order to take into account the effect of the disparity in energies 
of the two electrons of the pair, we simply replace H, in eqn. (2.15) for the 
equi-partition case, by the quantity H,,=—0-56H,. In the following we 
shall call the case as the average partition case. 


§ 3. EXPERIMENTAL PROCEDURE 


The electron—positron pairs on which measurements were made were 
found in five large emulsion stacks, in systematic scanning for electro- 
magnetic cascades. The total energy of each pair was determined from 
the development of the cascade in a separate experiment by Pinkau 
(private communication). In the present experiment, pairs with energy 
2 80 Bev and length per emulsion >1mm were used. We believe that 
the selection of the events is unbiased. Twenty-seven pairs were found 
satisfying the above conditions. The energy region covered by these 
pairs ranges from 80 to 2000 Bev. 

The ionization of each pair was determined by using the gap and blob 
method of Fowler and Perkins (1956), A gap-length of 4:16. was used, 
and, in addition, to reduce the statistical fluctuations, gap-lengths of 
2-50 » and 5-83» were also used on the first 750 u of the track of the pair. 

The ionization parameter, g, an dits statistical fluctuation were eval- 
uated using the formula 


g V(Nplog B/H,’ 


where B is the blob density: H, is the gap density for gaps >1, and N ‘ 
is the total number of gaps counted. In those cases where the ionization 
was determined using two different gap-lengths, the ionization was eval- 
uated for each gap-length, and the weighted mean was adopted with the 
smaller one of the two statistical errors. 
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The ionization was compared with the plateau ionization g, of a singly 
charged particle in the region of the emulsion near the origin of the pair. 
The magnitude of the ionization was normalized in units of go, after being 
corrected for depth effects, when it was necessary. 

The mean value of gp is 0-52, and although it varies from stack to 
stack and from plate to plate, the variation is not by more than a factor 2. 
There is no indication in the present experiments of any dependence of 
the suppression effect on the degree of development. 

While the total energy of the pair, was commonly obtained from studies 
of the resulting cascade, significant measurements of relative scattering 
were also sometimes possiblet. In such cases, the correlation between 
the disparity in energy of the electrons and the magnitude of the sup- 
pression effect was investigated. 


§4. EXPERIMENTAL RESULTS 


Since the existence of the suppression effect has been established by the 
pioneer work of D.H. Perkins, it appears important at the present stage 
to study its energy dependence, and to test the possibility of using it to 
determine the energy of electron-pairs. 


4.1 


Figure 1 shows the plot of the mean ionization against the distance from 
the pair-origin: the sample has been divided into three energy groups: 


Group I 17 pairs 80<H,< 200zBev, #,= 138 Bev, H,,= 77 Bev, 
Group II 6pairs 300<H,< 500Bev, H,= 425zev, H,.=240 Bev, 
Group IIT 4 pairs 1000< #,< 2000 Bev, H,=1450 Bev, H,,=810 Bev, 


E,, is the total energy of the pair determined by cascade measurements ; 
and H#, is the geometrical mean of H,’s. The solid and dotted lines in fig. 1 
represent the theoretical prediction given by (2.15), calculated for the equi- 
partition case, and the average partition case. It should be noted that 
the intersection of the asymptotic curve given by (2.15) with the axis R=1 
is determined by the value of v,,, whereas the slope of the curve is determined 
by the choice of Q. The best overall fit to the data, using the dotted curves 
is given by choosing the parameter & = 6-45, corresponding to a value of 
ro=72x102%cm. The form of the theoretical curves is in qualitative 
agreement with the observed variation in grain density. There is an 
obvious dependence of the suppression effect on energy as expected. 
Group III is worthy of special consideration, in view of the high energies 
involved. We have four pairs in this energy group. ‘Three of them (1, 2 
and 3) were found in the core of a jet of type 0+ 4p, together with a pair 
of lower energy, ~30Bev. The total energy of pairs in the jet was found 


+The author is indebted to Mr. K. Pinkau for his kindness in making the 
h * . . . . . 
results of his measurements available prior to publication. 
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to be 5000 Bev by cascade measurements: considerations of the relative 
angle of divergence of the pairs give the energy estimates quoi in fig..1, 
under the assumption that they have originated from two 7°-mesons in 
the core of the jet. The fourth pair occurs at the origin of a purely electro- 
magnetic cascade produced by a photon which originates outside the 
stack. The energy was determined rather unambiguously by cascade 
measurements. The observed ionization suggests deviations from a 
monotonic behaviour, but they are nevertheless attributed to statistical 
fluctuations. It should be emphasized that no appreciable deviation of 
the asymptotic value of the ionization from twice plateau value was 
observed, 
Fig. 1 


Group I 
80 £ Ec £200 Bev 
Ee =138\Bev|_ — 


Res 
Group IL 
300 $E.* 500Bev 


eae 


= 
7 


R=O0-5 { 


° 
R:1-0 a fm , } 
| i= 0 an Ss — 
Group Il Zz wv 18 
1000 £ Ec * 2000 Bev (4 pairs) 
Bee 
"etione ork) - pair (1)~ 2000 Bev 


pair 2)~ 1500 Bev 
pair @)~ 1500 Bev 
pair (4)~ 1000 Bev 


mean of 4 pains 


2 5) oe 2 5) i00dne ee 5 


The ionization versus the distance from the pair origin. The solid and dotted 
lines are evaluated by (2.15) for the cases of the equi-partition and the 
average partition of the energy, respectively. For Group III the 
ionizations of the individual pairs are recorded together with the mean 
values. 


4.2 


On the basis of the agreement between the observations and theoretical 
predictions in the last section, it seems worthwhile to try to estimate the 
energy of pairs by ionization measurements near the origin. To this end, 
it would be preferable to deal with the average ionization up to certain 
distances rather than the variation of the ionization over these distances. 
The average ionization is plotted against total energy of the pairs in figs. 
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2(a) and 2(b). The solid and dotted curves represent the theoretical 
values as given by formula (2.15), for the equi-partition case and the 
average case, respectively. In averaging (2.15) a lower limit of the 


Fig. 2 (a) 


5 1002 Ber Toco, 2 5 
E, Bev 
The average ionization up to 200 versus the total energy H, estimated by 
cascade measurements. The solid and dotted lines are due to the 


theoretical prediction (2.15) corresponding to the equi-partition and 
average partition of energy. 


Fig. 2 (5) 


Group IL Group IIL 


5 100 2 5 JO00 2 5 
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The average ionization up to 100, 5004 and 1000 1. The solid and dotted 
curves represent the theoretical prediction (2.15). 
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distance was taken at 10» to avoid the false contributions from the portion 
too near the origin, Figure 2 (a) shows the average ionization up to 200 
for all the pairs. In fig. 2(b) various lengths are adopted for different 
energy groups; here also the agreement is satisfactory. 

From these figures we can obtain the independent estimates of the 
energies of individual pairs, using the solid theoretical curves assuming 
equi-partition. 


Fig. 3 


{000 


100 


(a) The ratio of the energy estimated by ionization (#;) to that estimated by 
cascade measurements (H,) versus H,. Arrow + means that the 
most probable estimate of H,/H, is outside the scale; > means that 
its upper limit is also outside the scale. 

(b) The distribution of H;/H,. Five points fall outside the scale. 


The energy values H, estimated in this way are plotted in fig. 3 against 
the energy values H, determined by cascade measurements. The average 
ionization over the distances 200 1, 500 « and 1000 u were used for groups 
I, If and III respectively. For several pairs in group I, no significant 
decrease of the ionization was observed, this being masked by large statis- 
tical fluctuations. For these cases we have shown only the upper limits 
of the energy estimates. In this figure arrow (+) indicates that the most 
probable value of #,/E,, is outside the seale, and (>) means that even its 
upper limit is outside the scale. The histogram in fig. 3(b) shows the 
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distribution of the values #;/H,. H, coincides with #, within a factor 3, 
in 19 cases out of 27. Only one point falls in the region #,/E,>3. The 
distribution has a spread of the order of a factor 3, in either direction, 
around the median value of H,/H,~ 0-67. This distribution seems too 
broad to be explained by statistical fluctuations in gap and blob densities. 
It suggests that there is some real fluctuation in the suppression effect 
itself. In fact, there are several conceivable causes of the fluctuation. 
Multiple scattering results in a distribution of the relative separation 
given by (2.5). There may also be an effect due to the fluctuation in the 
initial angle of divergence for a given total energy}. In addition to these 


H,/E, vs Hs/E,. 
E;=the energy estimated by ionization measurements ; 
HL, =the energy estimated by cascade measurements ; 
E,=the energy estimated by relative scattering measurements. 


there is a fluctuation due to the disparity in the energies of two particles 
ofapair. In§2we considered the average effect of the disparity in energy 
partition and found that the effective energy of the average pair is 
E.,= 0-56H, instead of H,. However, the disparity in the energy, oud 


+ Burkhardt (1957) has given an interesting discussion about this point, 
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thus the effective energy, should fluctuate from pair to pair, giving a spread 
in the observed suppression effect for the pairs with a given total energy. 

Fortunately, in 13 cases, the results of the measurements of the relative 
scattering were available, and thus gave an estimate of the degree of energy 
disparity. For these pairs, we have studied the correlation between the 
results of relative scattering measurements and ionization measurements. 
The results are shown in fig. 4. The ordinate shows #,/H,, and the abscissa 
shows H,/E,. The ‘scattering energy’ H, was deduced from the data of 
relative scattering, making ce os assumption that both the 
particles were of equal energy (H,=24/2K/dre1, where K is scattering 
constant, ére1 scattering angle). ray circle represents one pair and the 
figures show the corresponding value of energy H, in Bev estimated from 
cascade measurements. Arrows indicate that the points are only upper or 
lower limits. Obviously, two quantities H,/H, and H;,/H, correlate 
strongly. We can also see that the points are along the straight line 
E,=E,, which will be expected under the circumstances where the relative 
separation is determined mainly by multiple scattering. In the present 
experiment, due to the poor statistics, we cannot decide whether there is 
an appreciable change in the initial angle of divergences according to the 
disparity in energy. 

However, it is rather clear that in the measurements of the suppression 
of ionization, one obtains information on a quantity of almost the same 
physical nature as the relative scattering. 

In conclusion, it may be said that the ionization measurements give a 
fairly good energy estimate for pairs with energy >300Bey. On the 
other hand, this method cannot be expected to give significant results for 
pairs with energy < 200 Bev, because of the small amount of information 
to be gained from the ionization over the first few hundred microns of the 
track. 

Among previous experimenters, Perkins (1955) has measured the 
ionization of six pairs. Their energies were estimated to be between 80 
and 400 Bev with a mean ~180 Bev. Some of the pairs which are used in 
these experiments were re-measured and included in the present experi- 
ment. Wolter and Miesowicz (1957) have measured the ionization of a 
single pair with somewhat higher energy ~1000Bev. They gave the 
mean ionization over the first 2604 and the next 1300. Both these 
results are consistent with those of the present experiment. 


§ 5. SUMMARY AND CONCLUSIONS 


Summarizing the above discussions, we would like to emphasize the 
following points : 


(1) The theory for the suppression effect of ionization by high energy 
pairs is in agreement with observations. 


(2) It is confirmed that the suppression is energy dependent, a fact 


which provides the basis for a method of determining the energy of pairs 
by ionization measurements. 
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(3) A strong correlation between the suppression effect and the dis- 
parity in energy between the particles of a pair has been observed. 

(4) The formula (2.15) gives reasonable basis for estimating the energies 
of pairs, at least for values >300Bey. This serves as a simple method 
for the energy determination, and is expected to be extremely useful in 
those cases where the available length of the pair in the emulsion stack is 
very short, or where there are several pairs of which the soft cascades 
cannot be separated from one another. The method seems to give informa- 
tion of a nature similar to that derived from the relative scattering 
measurements. At low energies the latter are superior to the former. At 
higher energies, however, the multiple scattering method is not applicable, 
because effects due to multiple scattering are easily masked by errors of 
measurements, and by effects due to the creation of bremsstrahlung by one 
or both members of a pair before the two component tracks are resolved. 
Reliable measurements are then difficult to make, and ionization measure - 
ments are therefore preferable: they become increasingly useful as the 
pair energy increases. 

It may be remarked that the accumulation of the data on pairs with 
energy > 1000 Bev is highly desirable, in order to separate the effects of the 
multiple scattering from those due to the initial angle of divergence. The 
detailed study of the distance dependence of the suppression for these 
pairs may contribute to our understanding of the mechanism of grain 
formation in emulsion. The examination of the suppression in over- 
and under-developed emulsions may be useful for the same purpose. 
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ABSTRACT 


The contours of the slip bands produced by fatigue straining of metals 
are viewed by a process of taper-sectioning and optical microscopy which 
permits effective magnifications of some 20000 times. The observations 
indicate ways in which the slip bands may turn into cracks and start disintegra- 
tion of the metal surface. They also show that the mechanism of deformation 
and failure changes as the amplitude of the cyclic straining increases, the 
general form of the fatigue S/N curve arising from superposition of two distinct 
mechanisms. The observations support the view that fatigue cracking is 
primarily a simple geometrical consequence of to-and-fro fine slip movements 
in the fatigue slip-bands. 


§ 1. [yTRODUCTION 


Eary investigators, notably Stanton and Bairstow (see Gough 1924), 
showed that fatigue failure results from fluctuating strains which contain a 
plastic component of small amplitude. Gough and his co-workers (1933) 
then showed that this component proceded by slip, which in cubic metals 
concentrates in bands, and that failure normally begins as a fine crack in 
those bands where slip is particularly dense. They concluded that the 
cracking was due to excessive deformation there. Later workers have 
amplified these observations. The fine structure of the bands is clearly 
shown, for instance, in beautiful photographs of aluminium published 
by Laurent (1952), Bullen (1952), and Forsyth (1952). Localizing of 
deformation in the bands is proved in a specially interesting manner by 
the observations of Hanstock (1957) on the over-ageing of fatigue bands 
in aluminium alloys and of Forsyth and Stubbington (1954) on the extrusion 
of material from such bands. 

However, though the general conclusion of Gough and his co-workers 
has been established, there are as yet few observations to show exactly 
how the localized deformation changes a fatigue slip band into a fine crack 
or to test the theoretical explanations that have been put forward. One 
experimental difficulty is that normal metallographic methods, while 
adequate for resolving a fatigue band laterally, as illustrated by fig. 2, 
Pl. 33, do not so easily show the structure in depth. The present work 
attempts to overcome the difficulty by observing the slip bands in section. 

Preliminary tests (Wood and Segall 1957 a) showed that for this purpose 
sufficient magnification could be obtained by taper-sectioning a specimen 
at a very small angle to the surface and then viewing the polished section 
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under high optical magnification. Thus a taper angle of 2° or 3° could 
magnify a surface disturbance in section by some 20 times, the cosecant 
of the angle of section. This, multiplied by optical magnifications of 
500 to 1000, gave effective magnifications of some 10 000 to 20 000, which 
appeared likely to show in some detail the disturbances created where a 
band of slip planes emerged from the surface. 


§ 2. EXPERIMENTAL PROCEDURE 


The general procedure was to prepare fatigue test pieces with initially 
smooth surfaces, protect the slip bands produced in straining by an electro- 
deposited coating, and then section. In the preliminary tests taper- 
sectioning was carried out by standard techniques. Part of a test piece 
was mounted in plastic with its surface inclined at a small angle to the 


Specimen C is clamped in jig by nut A. Surface disturbance 6R is magnified 
in section tot=d5R coseca=6Rh. R/x. (Not to scale.) 


surface of the mount. Parallel grinding of the mount then exposed the 
required section. This method, however, proved tedious; also it did 
not lend itself to precise measurement of the very small angles needed 
for high magnification. 

Consequently, for the following tests a simpler technique was devised. 
The test length of a cylindrical specimen after straining and plating was 
clamped in a Perspex jig, as shown in fig. 1, and the effect of taper sectioning 
obtained simply by grinding on the specimen a narrow flat. It will be 
seen that the plane of the flat is fixed geometrically by its lengthwise direc- 
tion and by the ball point of the screw 8, and that a radial disturbance at 
the specimen surface is magnified in section by coseca=R/x. It was 
found in this way that sections at very small angle could readily be obtained 
by making x small, and the angle still accurately found, for x could be 
measured by microscope. The method had also the advantage that 
several sections could be made on the one specimen merely by turning it in 
the jig. Also the magnifying factor could be varied along the one section, 
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if desired, by clamping the specimen with its length at a slight angle to the 
plane of the jig and so grinding a flat with a width 2” varying along its 
length. 

The metals used in the experiments described below were annealed 
OFHO copper and commercially pure alpha-brass containing 70-06% zinc. 
Specimens were prepared with test length 1-75 in. and diameter 0-2 in., 
dimensions chosen because of other work on the fatigue properties of these 
metals, though a larger diameter would have made it easier to obtain a 
high magnifying factor R/x. The specimens were finished by fine turning 
and final grinding on 600 grade emery, and after annealing were again 
smoothed by electropolishing in phosphoric acid. They were then sub- 
jected to fatigue straining, protected by a heavy plating of silver, and 
sectioned by grinding on a 600 emery surface. The sections were prepared 
for metallographic examinations by standard mechanical polishing with 
emery, diamond dust, and alumina, and by final etching in 10% ammonium 
persulphate. In the photographs, later reproduced, the magnifications 
are written in the form ¢ x m, the first being the theoretical taper magnifica - 
tions of a surface disturbance and the second the magnifications by micro- 
scope. 

The fatigue straining was alternating torsion applied in a Chevenard- 
type testing machine (1942), the straining being chosen to show how the 
fatigue bands developed during the progressive cycling at typical 
amplitudes. These details are given below as the metallographic 
observations are described, the amplitudes being expressed in terms of 
the angles of twist on the standardized test lengths. 


§ 3. INTERNAL FatTIGuE BANDS 


Examination of the sections was facilitated because etching brought 
out the internal slip zones. This was to be expected for it is known that 
zones of heavy deformation can often be etched up in sections of strained 
metals. Samuels (1955) in this way has recently identified internal slip 
bands in copper and copper alloys after direct straining, and Kemsley (1956) 
in ordinary sections of copper after fatigue straining. The effect in the 
present taper-sections is illustrated by fig. 3, Pl. 83, obtained from brass 
after 5x 10° cycles at +5°. This photograph is taken at relatively low 
magnification so as to show several grains and to confirm that the bands 
have the characteristics of slip, being parallel in one grain and differing 
in direction from grain to grain. 

Some properties of these fatigue bands are described later. The 
immediate point of interest is the disturbances of surface contour seen at 
higher magnification where the fatigue bands meet the specimen surface. 
These are described below. 


$4. NotcH aND PEAK EFFECTS 


It was found that the fatigue bands caused striking changes in local 
surface contour, They are conveniently classed as follows, 
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The first type is formed by fatigue bands which at the surface open into 
sharp notches. A typical example is provided by fig. 4, Pl. 33, taken from 
brass after 10° cycles at +5°, only 1/10th of the estimated fatigue life. 
Further cycling appears to increase the number of fatigue bands and of 
such notches, rather than increase the depth of each, and to cause final 
fracture by a linking of these small fissures, the crack jumping from one 
set to another. Typical linking is shown in fig. 5, Pl. 34, taken from a 
brass specimen which failed in another part of the specimen after 10° cycles 
at +5°. 

The second type is a converse effect, where relative slip movements in 
the fatigue band build up a peak or ridge. Figure 6, Pl. 34, taken from 
brass after 5 x 10° cycles at +2-5°, just before the fracture stage at the 
amplitude, shows a typical peak. The effect apparently arises when for- 
ward slip movements concentrate on one set of planes in a band and 
backward movements on a neighbouring set, so that the region in between 
is left protruding. 

The third type is a combination of the foregoing, both notches and peaks 
forming at the head of the same slip band, so resulting in an irregular con- 
tour. A representative example is shown by fig. 7, Pl. 34, taken from copper 
after 5x 104 cycles at +2-5°, approximately 1/20th of the expected life. 

The following are special cases of interest. The first occurs when fatigue 
bands happen to lie almost parallel to the surface. It is found that such 
bands can open into fissures also parallel to the surface; an example is 
provided by fig. 8, Pl. 35, taken from brass after 10° cycles at +5°. The 
second is when a grain is capable of cross-slip. Then fissures can develop 
along two directions in the grain. This is shown by fig. 9, Pl. 35, taken 
from the same specimen. It will also be noticed that the two sets of 
fissures can intersect and thus loosen a block of metal. In fact attention 
has been drawn to both cases because they illustrate particularly effective 
methods in which the fatigue bands can disintegrate the metal surface. 

The above observations appear to support the view previously put for- 
ward (Wood 1955, 1957) that a fatigue band is built up by to-and-fro 
fine-slip movements that can produce a notch or peak or any intermediate 
contour according as they become distributed across the band during 
progressive cycling. Initiation of cracks, particularly at the notches, 
then becomes simply a geometrical consequence of the to-and-fro fine 
slip peculiar to fatigue straining. The observations would also support a 
similar but more detailed process proposed by Mott (1958) in terms of 
dislocation theory. 


§ 5. Buock MOVEMENTS BETWEEN SLIP BANnps 
It is clear from the above observations that the fatigue bands may divide 
a grain into blocks and that the boundaries of the blocks are loosened when 
the bands become fissures. It is to be expected that the blocks in conse- 
quence may move relatively to each other, This was found to be a common 
occurrence, 
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The beginning of such movement in a typical case is shown by fig. 10, 
Pl. 35, taken from brass after 5 x 10° cycles at +5°, approximately half 
the expected life. Figure 11, Pl. 36, taken from the same specimen, 
shows the same kind of movement at a more advanced stage. Figure 12, 
Pl. 36, taken from brass after 10° cycles at this amplitude, shows a still 
more advanced stage, with a crack jumping from one fissure to the next 
and at the same time undermining surface blocks to such an extent that 
they can drop out. 

If the fatigue bands happen to be very close together then the block 
movements can build up a pronounced peak. An example is provided by 
fig. 13, Pl. 36, taken from the same specimen as the last photograph. 
Breakdown of such a peak evidently becomes a further effective process 
of surface disintegration. 

A special case of interest, which also amounts to relative block movement, 
frequently occurs at the boundary of neighbouring grains in which the slip 
movements happen to be inclined to each other at a small angle. Then a 
fissure forms at the boundary and one grain rears over the other. A typical 
example is given in fig. 14, Pl. 37, also taken from brass after 10° cycles 
at +5° (which caused failure in another part of this specimen). 

Thus the relative block movements, though facilitated by the fissure 
formations, in effect provide a second mechanism whereby disintegration 
and failure of a specimen may occur. It will again be noticed that 
these processes again begin at an early stage of the expected specimen life, 
being usually observed easily at 1/10th the life, although they appeared 
to accelerate in the final stages. 


$6. FURTHER PROPERTIES OF THE FATIGUE BANDS 


The following observations are selected because they appear relevant 
to the general problem of fatigue failure. 

(1) First it will have been noticed from the photographs that the etched- 
up fatigue bands are often limited in length and uneven in strength. This 
implies that the intensity of deformation is not uniform along the slip 
planes, so that only portions of the bands respond to etching. The effect 
was shown clearly in grains which exhibited cross-slip. The heavy 
deformation could start at one point on the first slip system and then continue 
on the second system, so creating a pattern of ‘ crazy slip ’ which made the 
irregular deformation particularly obvious. An example is given by 
fig. 15, Pl. 37, taken from brass after 7 x 10° cycles at +5°. ; 

Though such non-uniformity is to be expected, it is of special interest 
because of the further observation that as the amplitude of cyclic strain to 
which a specimen is subjected is made larger the active lengths of the fatigue 
bands become smaller, with the result that when the amplitude is large enough 
no slip bands at all etch wp. (‘This accounts for an observation by Kemsley 
(1956) who found that his etching techniques brought out fatigue bands 
only in specimens tested at what he termed low stresses.) The point is 
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also of interest in confirming suggestions (Wood 1957) that the mechanism 
of cyclic deformation should change with increasing amplitude and pure 
fatigue occur only when the plastic component in the cyclic strain is small. 
Further, it is of interest because the way in which the deformation leads to 
fracture should also change with increasing amplitude. This change in 
fact showed itself as follows. Whereas cracks at low amplitudes resulted 
in effect from the opening up of the fatigue bands, so following regular 
crystallographic paths, cracks at high amplitudes showed no such regular 
behaviour. Instead they formed at isolated regions, which could be points 
to which the fatigue bands had contracted, and at once branched out in 
quite irregular formations. These features will be clear from figs. 16, 
Pl. 37, and 17, Pl. 38, taken from brass which failed after 104 cycles at 
+ 15° and 10° cycles at + 20° respectively. The first shows at the surface 
merely relics of fatigue bands, and cracks which, though beginning among 
these bands, soon branch out quite irregularly. The second shows the 
complete absence of the fatigue bands and again the irregularly branching 
of cracks typical of the large amplitudes. This transition in mechanism 
of deformation and fracture began at amplitudes of about +10° in the 
brass and +7-5° in the copper. 

(2) A second group of observations deals with the changes produced 
when a fatigued specimen is heat-treated before the plating and sectioning. 
These observations were made to find whether the fatigue bands were 
zones of deformation of a kind, such as lattice distortions or vacancy 
concentrations, that could be dispersed by heating. 

In a typical experiment two specimens of copper were subjected to 10° 
and 5 x 10° cycles at + 25°, respectively 1/10th and 1/2 of their expected 
life, and each cut in two. One part was heated in vacuo for 1 hour at 
800°c, then plated and sectioned in the usual way. The other part was 
plated and sectioned without heating in order to provide comparisons. 
The comparison specimens showed the bands and contour disturbances 
to be expected from the fatigue straining. The heated parts showed the 
following changes. First the fatigue bands instead of disappearing broke 
up into links of what appeared to be cavities. This change is illustrated 
by fig. 18, Pl. 38 from the first specimen. Close examination indicated 
that a large part of the effect was certainly due to the fact that the bands 
in many instances must be deeply penetrating fissures which during heat 
treatment sintered into cavities at intervals along their length. For many 
examples were found of the kind shown by fig. 19, Pl. 38, taken from the 
second specimen, where the end of the fissure was rounded off and separated 
from cavities that clearly seemed initially attached. 

The second observation was that though heating rounded off the root 
of many fissures it never removed them. This is of interest in view of the 
important finding made first by Sinclair and Dolan (1951) that heating a 
fatigue test-piece at intervals during cycling does not extend its life, and 
their conclusion that irreparable damage to the structure is done in the 
early stages of straining. The present observations that deep persistent 
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fissures are formed early in the life of a specimen would support that 
conclusion. 

(3) A final group of observations confirmed that the intense fatigue 
bands were lines of weakness at an early stage. A section showing fatigue 
bands on specimens subjected to 1/10th of their expected life was bent 
through a small angle along its length in a four-point bending device. 
It was found that many bands then opened up into obvious cracks. An 
example from brass before and after bending is clear from figs. 20 and 21, 
Pl. 39, despite some unavoidable blurring of the latter photograph because 
of the curvature of the specimen. In general, the bands opening up most 
were those which before the bending had already formed obvious fissures 
at the specimen surface. 


§ 7. CONCLUSIONS 


Though observations need still to be obtained on the behaviour of other 
materials and the effects of other modes of cyclic straining, some initial 
conclusions may be drawn. 

First, the mechanism of deformation and fracture by cyclic straining 
depends on amplitude. The higher amplitudes, producing few internal 
fatigue bands and therefore few bands of localized high distortion, are those 
corresponding to the early part of an S/N curve where the curve falls 
sharply towards the N-axis. The low amplitudes producing intense 
fatigue bands are those corresponding to the later part of the S/N curve 
where it tends to parallelism with the N-axis. Accordingly, the special 
shape of the S/N curve results from the superposition of two distinct 
mechanism. This accords with other work on these materials (Wood, 1957, 
1958, Wood and Segall 1957 b). 

A second conclusion is possible in view of the other work referred to, 
which showed that the larger amplitudes produced progressive strain- 
hardening. It would appear then that the mechanism of failure at large 
amplitudes is akin to that of failure in static deformation, being started by 
stress-concentrations at a point where dislocations are piled up. The 
previous work also showed that small amplitudes in contrast produced 
relatively little strain-hardening, and therefore permitted totals of plastic 
deformation that would amount to some thousands per cent. The 
mechanism of failure in this case can scarcely be due to high stress-concentra- 
tions. At the same time it cannot primarily be due to high concentrations 
of vacancies or interstitials, for these would disperse on heat treatment. 
It appears therefore that the prime cause must lie in the notches and peaks 
built up by the geometry of the alternating slip movements in the heavily 
deformed fatigue bands; although concentrations of vacancies or other 
non-hardening point defects could well be a secondary cause, for they 
would allow the surface fissures to penetrate rapidly into such weakened 
ZONES. 

Finally it appears from the observations that numerous small but sharp 
notches are formed by at least one-tenth of the estimated life of a specimen, 
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Macroscopic failure is evidently a matter of time, occurring when the 
fissures become numerous enough to form the basis for a large-scale crackt. 
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ABSTRACT 


As the inherent crystal anisotropy of iron tends to zero above 600°C the 
initial magnetization above this temperature is considered in terms of purely 
rotational movements of the domain vectors from preferential directions 
determined by applied and internal stresses. Experimental results on the 
influence of applied stress on the initial permeability are described. It is 
shown that the results are in fair agreement with the theoretical predictions. 


§ 1. InTRODUCTION 


THE initial magnetization in iron at ordinary temperatures is ascribed to 
domain boundary movements essentially. The reason for this is the 
existence of high crystal anisotropy which results in negligibly small 
rotational movements of the domain vectors in low fields (Becker and 
Doring 1943). Asthe temperature of iron is raised the anisotropy constants 
decrease and tend to approach zero above 600°c (see e.g. Kittel 1949). 
This should lead one to expect that at temperatures between 600°C and 
the Curie point the initial magnetization process should be of a rotational 
type similar to the case of nickel at room temperature (Becker and Doring 
1943). Moreover, with the decrease of anisotropy energy as the tempera- 
ture is raised one would expect an associated increase of initial perme- 
ability—a well known phenomenon. 

In order to test how far the magnetic behaviour in low fields can be 
accounted for in terms of rotational movements above 600°C a tensile 
stress was applied to an iron specimen and the resulting effect on initial 
permeability observed. It so happens that the (111) magnetostriction 
in iron also disappears above 600°o (‘Takaki 1937, Lee 1955), and with 
the remaining single magnetostriction value the calculation of stress 
effects on initial permeability becomes possible for a polycrystalline speci- 
men. At room temperature, for instance, the two magnetostrictions of 
opposite sign, together with the high anisotropy energy, result in such 
complicated stress effects that no clear conclusions could be drawn from 
such measurements. 


§ 2. EXPERIMENTAL MErHop 


The specimen (Johnson and Matthey specpure iron) was in the form of 
a wire 1 mm in diameter and 30cm long. After annealing for 30 hours 
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in an atmosphere of wet hydrogen it was mounted vertically in a specially 
designed tubular oven which provided the necessary elevated temperatures. 
Throughout the experiment the specimen was kept in an atmosphere of 
purified hydrogen to prevent oxidation. Weights from the lower end 
provided the tensile stress. The magnetic measurements were carried 
out by ballistic galvanometer. A typical B-H curve at an elevated 
temperature is shown in fig. 1. As the galvanometer arrangement was 


Fig. 1 
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Insert indicates that the linearity of the initial portion lends itself to an accurate 
determination of the initial permeability. 


extremely sensitive (9 maxwell per mm of deflection) the initial linear 
portion of the curve could be determined accurately and a reliable value 
obtained for the initial permeability. After the application of a stress 
the specimen was first thoroughly demagnetized by a continuously 
decreasing alternating field and then the d.c. magnetizing field applied. 
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From the dimensions of the specimen one sees that end effects can be 


neglected. ote 1 
A typical curve illustrating the increase of initial permeability with 


applied stress is shown in fig. 2. 


Fig. 2 
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Curve showing variation of initial permeability with applied stress at 610-8°o. 
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[go = 1227 c.g.s. units is the value at zero stress. 


§3. THEORY 


With no applied tension or field we picture the magnetization of a 
domain as being held by weak internal stresses due to magnetostriction, 
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lattice imperfections, ete. See fig. 3 (a). On applying a magnetic field 
the spins will tend to swing into the direction of the field giving an overall 
magnetizing effect. If a tension is applied first, the spins will tend to set 
themselves in the direction of the applied stress, because at this temperature 
there exists only a single positive magnetostriction. The effect of an 
applied stress will thus be to increase the initial permeability. 

One can see that the influence of stress will be limited, because when all 
the spins are parallel to the stress axis as a result of stress, the magnetization 
will proceed as a result of 180° Bloch wall movements. The inhibitions 
offered to this movement by traces of impurities, lattice imperfections, 
ete., will not be decreased by stress. The decreasing effect of stress as the 


latter is increased is borne out by the flattening out of the curve as seen in 
fig. 2. 


Fig. 3 


(0) 


Schematic presentation for iron above 600°c. 


Referring to fig. 3 (6) for calculation purposes : The total energy due 
to magnetic field H, internal stress o; and applied external stress o, is 
given by 

F =3)o; sin 7B + 3Ao sin? 6— HJ, cos 6 a0 Lar sae) 


where A, the (100) magnetostriction, is now treated as an isotropic magneto- 
striction for the polycrystalline specimen. As explained before, the 
anisotropy energy is negligible. Equation (1) must be minimized to 
establish the equilibrium condition and hence the expression for the 
permeability is obtained. We follow a mathematical technique employed 
by Rathenau and Fast (1955). The contributions of field energy and 
external strain energy are considered to be small compared with that due 
to the internal strains. This assumption is justified if the applied magnetic 
field and external stress approach zero. 
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Equation (1) leads, with $A0 < 8A, and HJ, < 3Ao;, neglecting quantities 
of the fourth order, to 
F = 820 sin® «— HJ, cos €+ $Aa,8? 


— 3Aa(B sin 2€ — B® cos 2e) 
—HJ,{B sin «—($P*) cos €]. 


The equilibrium angle follows from 07 /d0B=0, giving 


B= (ZAo sin 2e + HJ, sin e) (2) 
(3A0,; + 3Ac cos 2e + HJ, cos €) eds 

We consider a domain in which the direction of the spins, in the absence 

of field and stress, is given by angle «. Upon application of stress and 


field unit volume contributes J, cos @ to the magnetization along the 
direction of H. ‘To the intial susceptibility it contributes. 


K(€) = Baw cos («—B) | 


H-0 


= | sin («—f) sa if 


which, bearing in mind that we are dealing with initial susceptibility, for 
small f is 


e(e)=J,| (sin e—f cos sf ls «” pee ES 
Substituting (8),,,, and (08/0H),,,, from (2) one obtains 
AGS J 2 sin? € (1 z 3Ac Cos? € [1 —(o/o;) sin? €} 
3X0; 3Ao0; + 3Ao cos 2€/ [1+ (c/o;) cos Ze]? 


which for small c/o; reduces to 


_ J2sin® 
3,0; 


2 


K(e€) [1 \(a)o;) (l= 4 cos? <) aa meee) 


Assuming 3Ac independent of « and for an isotropic distribution of e: 


K(e)=4 ie K(e) sin € de 


pede lfo | lfo 
= Fat! 5(Z)f aa + 3(Z)h- amare 


Substitute the expression xy=2.J,2/9\c;, which is obtained from the 
above if c=0. 
Thus we can write 


oe - VWkoo 
Ko LOw * 
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and for permeabilities, since u=4ar«+1, p>1 


Boy DApiyo 
Lo 407J ,? 
=l+ko Spr fet Iabue A tay (8) 


and the graph of ju/4) vs o (for small values of c) should be a straight line, 
with slope &, cutting the ordinate at unity. 


§ 4. RESULTS AND DiscussIon 


Using the values of Takaki (1937) to plot a curve of magnetostriction 
versus temperature and reading off values corresponding to our tempera- 
tures, the values for k were calculated from eqn. (6) and compared with 
the measured values taken from curves such as in fig. 2. The results 
are given in the table below: 


k Measured k Cale. 
x 10° x 109 


2-432 2-014 
3183 2-412 


4-050 3-239 
4-999 3°82] 
3-841 5363 
3°913 6-531 


The fair agreement between theory and experiment in the region 600°C 
to 700°c indicates that the major part of the initial magnetization process 
in this region is accounted for. As the Curie point is approached the 
deviation between theory and experiment becomes more marked. This 
is not entirely unexpected as we have no adequate knowledge of the value 
and behaviour in this region of the magnetic constants involved. 

If, instead of being caused by spin rotation, we consider the magnetization 
as being due to reversible wall displacement, we must fit the readings 
with formulae of the type. 

k= Kol + $(o/04)} 
due to Rathenau and Fast (1955) or 
K=Ko{1+3-2(a/0;)} 


as found by Bozorth and Williams (1945) for 45 Permalloy. Both of these 
formulae give very much bigger differences between theory and experiment. 
Thus the original picture of spin rotation seems to fit the facts far better 
than that of wall movement. 
The formula for the initial permeability without zero stress, as obtained 
from eqn. (5) is y= (87J,2/9Ao,) +1. It is the same as that of Becker 
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and Déring (1943) obtained from direct derivation for rotational 
processes. Using the value of p, at 610-8°c the calculated value for 
o,= 12-4 x 107dynes/em?. Taking an estimated value for Young’s modulus 
of 1-4 x 10!2dynes/em? at this temperature (e.g. see Smithels 1955), the 
value of the magnetostrictive stress is approximately 4 x 10’ dynes/cm?. 
This indicates that the order of magnitude obtained for the internal 
stress is reasonable if we assume that the residual stresses normally 
present in a metal have been mostly relieved at the high temperature. 

Concerning the steady rise of initial permeability with temperature and 
the rather rapid increase as the Curie point is reached, Fahlenbrach (1948) 
proposes that the randomness of domain vector orientation increases with 
temperature and causes the observed rise. In the light of the foregoing 
discussion such a permeability rise can be pictured quite simply. With 
zero crystal anisotropy and the diminishing of internal magnetostrictive 
strain as the Curie point is approached less and less resistance is offered 
to reorientation of the domains by an applied magnetic field and con- 
sequently a rise in the initial permeability. 
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ABSTRACT 


A hole theory of fluids is described in which the potential of a molecule 
in a cell is not assumed to be spherically symmetric about the cell centre. 
Although this theory is incomplete in several respects, it has lead to the 
discovery of a function of pressure, density and temperature which when 
evaluated for real fluids, such as argon and krypton, appears to be a function 
of density only, and independent of the state, liquid or gas. It is shown 
how this property would lead to an equation of state valid for liquid and 
gas, and the results of a first attempt in this direction are given, which are 
in better agreement with experiment than those of other lattice or hole 
theories. 


§ 1. INTRODUCTION 


This paper describes a simple but as yet incomplete hole theory of fluids, 
in which the potential energy of the molecule in a cell is not assumed to 
be spherically symmetric about the centre of the cell. The cell dimension 
is obtained by minimizing the free energy. In the course of the develop- 
ment of this theory, a parameter was introduced, whose value as a 
consequence of the theory could be obtained from the experimental 
values of pressure, density and temperature. The value of this parameter 
was evaluated in this way for argon, and it was found to have the important 
property of being a function of density only, independent of the state, 
liquid or gas. That such a function of pressure, density and temperature 
should be a function of the density only for a real fluid should have 
important consequences in the theory of fluids. A knowledge of its 
dependence on density, for example, leads immediately to an equation of 
state valid for liquid and gas and consequently to the evaluation of other 
thermodynamic functions. A first attempt at such an equation of state 
is described, and represents an improvement in experimental agreement 
over others obtained from cell er hole theories. 

The present theory differs from previous hole theories (Cernuschi and 
Eyring 1939, Ono 1947, Peek and Hill 1950, Rowlinson and Curtiss 1951) 
in the expression used for the generalized free volume per cell (see e.g. 
Hirschfelder et al. 1954, p. 313) and in the manner in which this quantity 
is introduced. In these earlier theories the logarithm of the generalized 
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free volume is taken to be linearly dependent on w, the fraction of nearest 
neighbour sites about the central site which are vacant (see Hirschfelder 
et al. 1954, p. 314, or Rowlinson and Curtiss 1951) ; in the present work 
the generalized free volume, not its logarithm, is taken to be a power 
series in w. Further this quantity is introduced by means of a simple 
modification of the expression for the free energy used in the cell theory 
of Lennard-Jones and Devonshire (1937). 

The theory is incomplete in that it does not yet provide an explicit 
form for the above-mentioned function of density, and in that it lacks a 
completely rigorous deduction from the fundamental equations of 
statistical mechanics. The principal purpose in publishing this work is to 
describe the properties of the above parameter in the hope that more 
interest will be stimulated in hole theories of fluids, and in particular, 
that a rigorous development of the behaviour of this parameter for real 
fluids will result. 


§ 2. Toe THEORY 


In studying the radial distribution functions for argon obtained from 
x-ray diffraction experiments (Hisenstein and Gingrich 1942), one is 
impressed by two facts—that the position of the first coordination peak is 
practically independent of the density, and that the number of nearest 
neighbours (proportional to the area under the first coordination peak) 
varies almost as does the density. Referring these results to a lattice 
theory where the molecules are to be assumed to lie on or near lattice 
sites, such a theory should produce a lattice dimension which varies very 
little with density, the change in the density to be accompanied by a 
corresponding change in the number of vacant sites. The early theory of 
Lennard-Jones and Devonshire was unrealistic in these respects. 

In taking account of vacant sites, the hole theories mentioned in § 1 
obtained a free-volume per cell whose logarithm was linearly dependent 
onw. In some cases this was a result of assuming a spherically symmetric 
potential in each cell, as if the nearest neighbours were smeared out 
uniformly over a spherical surface ; and in others this dependence was 
assumed. It is probable that such an assumption is related to a 
spherically symmetric potential. Such a smearing out of nearest neigh- 
bours is possible related to a physical picture of a fluid in which the 
molecular motions are sufficiently large and random so that a central 
molecule is always prevented from approaching its cell wall. This is 
very well illustrated by the cellular potential of Lennard-Jones and 
Devonshire which becomes infinite on an arbitrarily defined spherical 
surface which plays the part of a cell wall, and also by any theoretical 
modification which simply reduced this potential everywhere by the 
fraction of occupied sites. An equally plausible physical picture, justified 
by the high degree of short-range order at high densities, as exhibited 
by x-ray diffraction experiments, is one in which the emphasis is placed 
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on the lattice structure, whereby this structure is preserved, at least 
locally, even when vacant sites are produced ; that is, when a nearest 
neighbour site is vacant, the other nearest neighbours are to a certain 
extent excluded from this vacant site on account of the tight packing. 
On the same reasoning the central molecule will also be excluded, but 
nevertheless its free volume in the cell will have been increased in the 
neighbourhood of this vacant site. Such a picture represents the extent 
to which this present theory departs from the concept of spherical 
symmetry in the cell. The above idea thus leads naturally to a free 
volume which is linearly dependent on w, and it was assumed at first that 
the free volume of the cell was the cell volume multiplied by w. When 
all sites are occupied the free volume could be taken as that of Lennard- 
Jones and Devonshire, but this is such a relatively small quantity that 
it was neglected, though in principle there would be no difficulty in 
including it. The potential at the centre of the cell is taken to be 
proportional to 1—w, the fraction of occupied sites. 
The Helmholtz free energy per molecule, f, was taken to be 
f=(1—@)b(x) — kT In (w/yx), oe ae CL) 
here the terms in f depending on temperature only are not included as 
they play no part in the minimizing process and in the calculation of the 
pressure ; w is the fraction of unoccupied sites ; «=1/a?, where ais the 
lattice constant ; ya is the reciprocal of the cell volume (for argon the 
lattice structure is taken as face-centred cubic for which y= 1/2); 
(1—w)(x) is the potential at the cell-centre assumed to be proportional 
to the fraction, 1 — w, of occupied sites. The number density, 7, is clearly 
given by 
n,=(l—w)y2, Pe ree ns ee ee) 
enabling (1) to be written 
f= (y/yx)p(a) — kT In {(ye—m)/[(y2)?j- (8) 
This expression enables x to be determined by minimizing f with respect 
tox. For argon, (x) was taken to be a Lennard-Jones type of potential, 


corrected for the attraction of more distant neighbours than the nearest 
(Lennard-Jones and Devonshire 1937), i.e. 


h(x) = 6e{(r9/a)12 — 2-4 (r/a)%}. Sr astis SCL) 


The results were promising in that a was always close to 79, and the 
number of nearest neighbours, 12(1—w), agreed reasonably well with 
those obtained experimentally by Eisenstein and Gingrich (1942). An 
equation of state was evaluated and the correct behaviour was observed. 
However the liquid densities at saturation vapour pressure were too near 
the crystal density for good experimental agreement. It was then that 
an expression for the free volume which is quadratic in w was introduced, 
purely as an assumption. The following expression is now taken as a 
starting point 

fH(1—e)f(@)—kT in (ow + w2)/yx}, + + (8) 
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the coefficient of w? being taken as unity, while for convenience the 
arbitrarily introduced parameter a was taken as the coefficient of w. 
It was hoped that « might be a function of temperature only, and that 
therefore any necessary constant of proportionality in the expression 
for the free volume would also be temperature dependent only, and 
therefore in the minimizing process and in the evaluation of the pressure 
it would be superfluous. This proved to be wrong, but it will be shown 


Fig. 1 
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© Argon, liquid or vapour at saturation vapour pressure ; © argon, 7')= 1-083, 
rE eS i rn a 
liquid ; x argon, 7)=1-25, liquid; @ argon, 7')>=2-276, gas ; + argon, 
Ty=4:167, gas. C is the critical poimt for argon. (krypton, 
T)=1-597 and 2-767, gas. 


later how the procedure may be made consistent. Thus, assuming a to 
be dependent on temperature alone, f may be minimized with respect 
to a, using (2), and an expression for the pressure obtained, resulting in 
the following two equations finally, 


fr -1=(2)(ovonery . .  . .  ) 


Clomptiteres() 


This pair of equations enables w to be evaluated by eliminating x and 
using experimental values of pressure, density and temperature. For 
argon (Michels e¢ al. 1949, Din 1956) the values of In (1 —«) thus obtained 
are shown in fig. 1, plotted against V,/V, where V,=7,3/\/2, and V is 
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the reciprocal of the number density ; ¢« was taken to be 120k, and 
ry 3834, Lennard-Jones (1924). Also shown in this figure are a few 
points for krypton, derived from the experimental data of Beattie ef al. 


(1952). The table gives some numerical values. The reduced temperature 
Ty) =kT'/e is used. 


State (19/v) ie In (1 — «) 
Fas 0-051 2-276 0-65 
" 0-051 3-526 0-65 
Saturated vapour 0-103 1-159 0-60 
Gas 0-103 2-276 0-60 
“ 0-103 3-526 0-60 
5 0-112 4-167 0-60 
x 0-340 2-276 0-41 
> 0-340 3-526 0-41 
Liquid 0-468 1-233 0-32 
Gas 0-473 2-276 0-31 
» 0-608 2-276 0-22 
s 0-608 3-526 0-23 
Liquid 0-649 1-083 0-191 
- 0-647 1-25 0-195 
Gas 0-677 2-276 0-182 
. 0-677 3-526 0-191 


Although the original assumption is not borne out, fig. 1 shows that « 
definitely appears to be a function of the density only, and to be 
independent of the state. It should be emphasized that this behaviour 
of ~ is an experimental fact, however eqns. (6) and (7) may have been 
derived. A way of removing the inconsistency that in deriving (6) and (7), 
a was assumed dependent on 7’ only, and proves to be a function of 
density only, will be indicated later. 

If the functional dependence of « on the density were known, and 
assumed to hold everywhere, an immediate application would be to the 
derivation of an equation of state valid for both liquid and gas. A first 
attempt in this direction has been made by determining empirically from 
fig. 1 the expression for In (1—«). An obvious method would be to use a 
‘least squares’, method, but it was decided not to do this as the 
experimental errors appear to be greatest in the high density region. The 
expression = In(1 — a) = In 2f(1— w/v) +O(1—vp/v)R/(1 +6) . . . (8) 
was used, as it is clear from eqns. (6) and (7) that a>+—1 as 1,>0; 
also from fig. 1 it appears highly probable that «+0 as vv). 6 was chosen 
equal to 0-4 so that the curve of (8) would pass through the value of 
In(1—a) at the critical point. The equation of state was then derived 
by assuming (8) solving (7) for a at a given density and temperature, 
and finally by substituting this value of x in (1). Figure 2 shows several 
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isotherms thus obtained. The critical isotherm, 7')=1-3 (1-256 experi- 
mental for argon) is compared with experimental points. The agreement 
is very good for such a simple expression for In(1—a). The isotherm 
T,=1, yields v)/v=0-670 (experimental 0-708) for the liquid and 0:04 
(experimental 0-04) for the vapour at a saturation vapour pressure of 
14 atmospheres (experimental l2atm) ; the compressibility of the 


Os 
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a, b, c and d are isotherms for 7')=1, 1-258, 1-3, and 2-276 respectively. Fora 
the horizontal line giving the saturation vapour pressure, according to 
Maxwell’s law of equal areas, is shown. x denotes experimental 
points for the critical isotherm of argon at 151-9°K ; © denotes experi- 
mental points for argon at 0°o (1) =2:276) ey 
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liquid at s.v.p. is calculated to be 0-0009 per atmosphere (experimental 
0-0006). An isotherm at 7')=2-276 (for argon, 0°c) is shown to be in 
excellent agreement with experiment. 

It should also be mentioned that the number of nearest neighbours, 
calculated from (2), taking the value of w obtained when solving (7) using 
(8), is in good agreement with the values obtained by Eisenstein and 
Gingrich (1942), from x-ray diffraction patterns. The corresponding 
value of a, the lattice constant, always differs from 7) by less than about 
8%, a result which is consistent with experiment. 


§ 3. IMPROVEMENT OF THE THEORY 


The difficulty mentioned, that « is a function of the density, raises the 
question as to the consistency of the theory. To avoid this it must be 
assumed that « is a parameter which, analogously to w, must be evaluated 
by minimizing the free energy. To make this possible the free volume 
should have another factor, such that its complete expression is 


m(a,T’)(aw + w?)/(yx), No eg ae te I), 


the nature of the function m not being known at present. If it is assumed 
that « must be determined by minimizing the free energy with respect to 
a, then as Of/dx is zero, « contributes nothing to the pressure. 

The physical significance of « is not known, but it may be that it is a 
sort of structure factor like w, such that it determines the arrangement 
of pairs of vacancies. It is possible, for instance, that in an arrangement 
where two contiguous nearest neighbour sites are vacant the local 
structure tends to break down, leading to a smearing out of nearest 
neighbours and making a free volume, similar to that of Lennard-Jones 
and Devonshire, the appropriate one. Work is proceeding on this aspect. 
It is worth mentioning that the free energy of the fluid could be calculated 
by integrating the equation of state from zero density, where the free 
energy is know, to any desired point on an isotherm, and thus the value 
of m could be obtained. 


§ 4. CONCLUSION 


A relatively crude theory has been described, whose main virtues are 
that it indicates a new approach to a hole theory, and that it has been 
the avenue for the discovery of an important property of real fluids, 
leading to a useful equation of state depending on one constant (b of (8)) 
over and above those of the Lennard-Jones potential. Further empirical 
work is indicated in endeavouring to find central cell potentials, (x), 
for more complicated fluids such as carbon dioxide, nitrogen etc., and 
in determining a more satisfactory empirical form for In(1—«), which 
would probably improve the agreement on the compressibility of liquid 
argon. On the theoretical side, much remains to be done in elucidating 
the physical significance of the hypothecated expression for the free 
volume stated in (9). 
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ABSTRACT 


A method is proposed whereby the light emitted by liquid scintillators may 
be increased by application of a suitable a.c. field. It is calculated that by 
this means the light output may be increased amply to allow minimum 
ionizing tracks to be recorded photographically without the possibility of 
instability arising. 


§ 1. INTRODUCTION 


Various drawbacks in the bubble chamber method of obtaining photo- 
graphs of nuclear particle tracks make a consideration of other methods 
desirable. In this investigation two main points are considered : firstly, 
techniques to increase the ionization along the track of a nuclear particle ; 
secondly, methods of observing ionization. 


§ 2. INCREASE OF IONIZATION ALONG THE TRACK OF A NUCLEAR PARTICLE 


Two distinct thresholds occur in a dielectric as the a.c. field potential 
is increased : the lower threshold is at about 42°% of the breakdown voltage 
when ionization commences in the dielectric; the higher threshold is when 
breakdown of the dielectric occurs (Bute et al. 1954). For high frequencies 
this ionization does not occur randomly in the dielectric but occurs along 
the tracks of ionizing particles and indeed breakdown streamers develop 
from these (Prowse and Lane 1956). 

Although for gases, consideration of negative ion mobilities indicates 
that fields of microwave frequency are required to confine the increase in 
ionization to a 10-3 em diameter track ; for liquefied gases, fields of fre- 
quencies of a few Mc/s should be satisfactory. 


§ 3. PHOTOGRAPHY OF IONIZATION 


For the observation of the ionization direct photography of the ‘ scintilla- 
tion’ light, emitted during decay of excitation and recombination of 
ionization, is possible. The liquefied noble gases are known to be good 
scintillators having scintillation efficiencies comparable with Nal-Tl 
(Northrop and Nobles 1956). For liquid argon it is estimated that the 
ionization along the track must be increased by about 10° in order to 
obtain photographs of minimum ionization. Using fields greater than 
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the breakdown threshold, inception of breakdown occurs when the number 
of ions have increased by about 108 (Gould and Roberts 1956). Hence 
the above required magnification can be obtained without breakdown of 
the dielectric and will occur in about 9 cycles of the high frequency electric 
field. Alternatively, using field strengths less than the breakdown 
strength, but in excess of the ionization threshold, a slower, but completely 
stable method of magnifying the light output may be obtained. 


§ 4. PRacTicaAL ASPECTS 


Ionizations occur in liquid argon and liquid hydrocarbons in fields of 
100 kv/em and 400kv/em respectively (Williams and Stacey 1957, 
Bute et al. 1954). Large chambers will be possible since Tesla coils and 
resonance transformers have been developed which are capable of producing 
radio frequency fields at several million volts (Craggs and Meek 1954). 

Heating of the dielectric volume is expected due to the presence of the 
electric field. However, comparison with dynamic losses of liquid hydrogen 
bubble chambers indicated that the loss expected for the electrical method 
is 10-* that obtained with bubble chambers. 

The velocity of the particles can be determined from the degree of 
ionization. The initial ionization produced by the nuclear particle will be 
magnified by a factor 2” where 7 is the number of ionizing collisions, related 
to the number of cycles of the field. Ifthe field is uniform, then throughout 
the chamber n will be the same, and the intensity of light emitted along 
a track will be nearly proportional to the initial number of electrons. 


§ 5. ADVANTAGES OF THE METHOD AS COMPARED TO BUBBLE CHAMBERS 


The method proposed has certain definite advantages over the use of 
bubble chambers : 


(a) Counter control is possible. 

(b) Cycling rate can be very high. 

(c) No illumination or mirror system is required. 

(d) Experimentally the method is adaptable for the use of various 
liquids in one chamber. 

(¢) Comparatively low dynamic loss is a definite advantage when large 
chambers are required. 


(f) Thermal gradients of the liquid only affect stopping power through 
density, unlike in bubble chambers where they also affect bubble nucleation. 
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ABSTRACT 


Experimental evidence is presented which lends direct support to the 
Stroh mechanism for crack formation. After a magnesium oxide crystal 
was compressed about 3%, it was observed that minute cracks had been 
formed in the region of a kink band generated by non-homogeneous plastic 
strain. These cracks lay on conjugate {110} slip planes and not on the usual 
{100} cleavage planes. They extended in the form of [001] direction slits 
running right through the crystal. 


§ 1. INTRODUCTION 


Srrow (1954, 1955, 1957) has analysed the nucleation of a crack in terms 
of the stresses associated with edge dislocations piled up against a barrier. 
On the basis of this analysis it is to be expected that there will be three 
factors which govern the formation of cracks in a solid. First, the density 
and distribution of the obstacles responsible for pile-up which are already 
present in the material; second, the ease with which new barriers can be 
formed during plastic deformation ; and finally, the ability of the material 
to relieve the stress field due to a pile-up by further local plastic flow. 

That grain boundaries in polycrystalline materials can be effective 
barriers has been demonstrated by Gilman (1957) using zine bicrystals. 
As far as single crystals are concerned, however, little is understood of the 
nature of the barriers. Sessile dislocations formed as a result of dislocation 
interaction are not sufficiently stable to support a large pile-up (Stroh 1957). 
Also, Stokes (1957) has shown that a fixed array of small angle boundaries 
(formed by annealing a dimple in a zine single crystal (Washburn 1956) 
will not support a sufficient pile-up of dislocations to produce a crack 
during simple shear fatigue. 

On the basis of the Stroh theory, it is the shear stress forcing dislocations 
into the head of the piled-up group which leads to crack nucleation. So 
far as crack growth in a single crystal is concerned, there are then two 
possibilities depending upon the nature of the external stress field. Ina 
tensile stress field, providing the Griffith energy and the crack velocity 
criteria (Mott 1948, Gilman 1956) are satisfied, microcracks will extend 
and cleavage failure will occur catastrophically. In a compressive stress 


y+ Communicated by Dr. A. H. Cottrell, F.R.S. 
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field the microcrack will remain stable. For this reason it seems advisable 
to study the problem of crack formation in single crystals during compres- 
sion. 

In the work reported here, cracks have been found to originate in kink 
bands, produced as a result of non-homogeneous compression, in single 
crystals of magnesium oxide. Crystals of this material have the advantage 
of being transparent and exhibit stress birefringence under polarized 
illumination, thereby enabling the internal stress pattern to be examined 
following plastic strain. 


§ 2. EXPERIMENTAL TECHNIQUES 


The specimens used in this investigation were cleaved from commercial 
purity (96-99% MgO) magnorite crystals purchased from the Norton 
Company. The bulk crystals were slightly amber in colour, and contained 
a number of small boundaries, otherwise there were no microscopic flaws. 

Specimens, having approximate dimensions 10 mmx 2mm x 2mm, 
were deformed in a conventional hard testing machine arranged for 
compression, after which they were examined on the metallograph. 

A successful etching solution was developed for revealing individual 
dislocations in this material. The solution consisted of: 5 parts of 
saturated ammonium chloride, 1 part of concentrated sulphuric acid, and 
1 part of distilled water. When the crystals were immersed in this solution 
for }-hour at room temperature, conical etch pits were produced. Etch 
pits on two crystal surfaces exposed by cleavage matched one another 
exactly as shown in fig. 1, Pl. 40, thus confirming the selectivity of the 
re-agent. 


§ 3. EXPERIMENTAL OBSERVATIONS 
3.1. Plastic Flow in Magnesium Oxide at Room Temperature 


By combining metallographic and etch pit techniques, the slip system 
in magnesium oxide at room temperature has been confirmed to be over 
{110} slip planes in (110) slip directions, as previously quoted by Schmid 
and Boas (1935). This is a typical slip mode for crystals having the rock 
salt structure. 

Residual stresses within these transparent cubic crystals could be 
detected with polarized illumination, so that the slip surfaces were 
delineated when viewed edge on (Nye 1949). The etch pit distribution 
pattern confirmed that the stress birefringence resulted from concentrations 
of edge dislocations on the active slip planes. 

Deformation at room temperature led to the appearance of many slip 
bands of widely varying band width. These slip bands were characterized 
by being exceptionally straight, similar to those produced on sodium 
chloride for example. Only at very high magnifications (x 1000) could 
individual slip lines be resolved within them. At this magnification there 
was a complete absence of cross slip. 
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The inherent plasticity of magnesium oxide was readily demonstrated 
by the following experiment. Particles of 46 mesh silicon carbide were 
sprinkled from a height of only three inches. Rosettes, of the kind studied 
by Gilman and Johnston (1957) on lithium fluoride crystals, appeared on 
etching in the ammonium chloride solution. These rosettes are due to slip 
dislocations generated under the impact of the tiny carborundum particles 
(of approximate weight 3-6 x 10-4 g). Loading a crystal which had been 
previously sprinkled in this way resulted in the expansion of the dislocation 
half loops, and as fig. 2, Pl. 40 shows, each slip line passes through a 
rosette. A more detailed study of the behaviour of individual dislocations 
in magnesium oxide will be published at a later date. 

Recent observations indicate that the mechanical properties of mag- 
nesium oxide crystals are sensitive to their surface condition (Lad et al. 
1958, Gorum et al. 1957). In particular, workers at the University of 
California and elsewhere have shown that freshly cleaved crystals are more 
ductile than those which have been standing in air for a length of time. 
In a study of crack formation under a three-point bending load carried out 
in our laboratories, it was found that specimens always plastically deform 
to a certain extent prior to fracture, irrespective of the time for which the 
surfaces had been exposed to the atmosphere. While the plastic strain 
cannot necessarily be detected macroscopically, it is readily visible under 
the microscope. 

Although it is recognized that the state of the surface is an important 
variable in the fracture behaviour of magnesium oxide, in the work reported 
below the surface condition was essentially constant. Specimens were 
allowed to stand in air for six days following cleavage from the bulk crystal 
before being deformed. 


3.2. Crack Formation 


Figures 3 and 4, Pls. 41 and 42, are photomicrographs taken with polarized 
illumination on a specimen after compression in the testing machinejy. 
It will be noticed that according to the stress birefringence striations in 
fig. 3 the crystal lattice was kinked in this region about an [001] axis of 
rotation, with the plane of the kink parallel to the (110) plane. 

Careful examination of the external surface in this kinked region enabled 
two kinds of cracks to be distinguished. One type was characterized by 
its short length and the fact that it was situated at the boundary of the 
kink band. It lay parallel to the [110] direction and thus to the trace 
of the kink band, as indicated by the arrows in fig. 4 (e). This type of 
crack, for reasons to be given in the discussion (§ 4), will be referred to as a 
“Stroh-crack’ to distinguish it from the other type, which will be called a 


7 This particular specimen had buckled during compression so that the 
strain was distinctly non-homogeneous. These photomicrographs were taken 
in the most severely compressed region of the specimen. The crosshead move- 
ment was equivalent to an overall strain of 3°%. 
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‘secondary crack’. Secondary cracks lay for the most part parallel to the 
[110] direction and were situated within the kink band, as shown in fig. 3. 
Most of the secondary cracks also had a short segment parallel to and in 
line with the Stroh-cracks, all delineating the boundary (A—B) of the kink 
band. It is significant that neither the Stroh-cracks nor the secondary 
cracks lay along the normal (100) cleavage planes, instead they were formed 
on the conjugate (110) and (110) slip planes. One of the secondary eracks 
in fig. 3 (top, right) however, had extended a short distance along the (010) 
plane at one tip. 

A most interesting observation was made when the metallograph was 
focused beneath the surface of the crystal. It was found that the cracks 
did not exist merely on the surface but continued right through the speci- 
men. They were, in fact, tiny slits extending along the [001] direction 
normal to the crystal surface, as illustrated in fig. 5. The series of photo- 
micrographs in fig. 4 were taken at different levels beneath the surface 
using polarized illumination. These show that the two Stroh-cracks 
extended all the way from one surface to the other. On the other hand, 
the secondary crack did not, it was tapered and stopped within 20 , of the 
far surface. 


Fig. 5 


i SLIP PLANES 


Crystallographic configuration of the slits in fig. 4. 


The formation of these slits in compression offered a rare opportunity 
to observe this phenomenon, for once having formed in a region where the 
local stresses were tensile there was no tendency for the crack length to 
increase beyond this region, since the macroscopic stress was compressive. 
However, if Stroh-cracks nucleate under a macroscopic tensile stress in this 
material brittle failure may be catastrophic and such tiny slits will be 
difficult to observe. 

A survey of the whole crystal surface indicated that not all of the cracks 
extended beneath the surface. The one shown in fig. 6 (@), Pie43 tor 
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example, was only a surface crack. This photomicrograph was taken on a 
region quite separate from that of figs. 3 and 4. Figure 6 (6), Pl. 48, was 
taken on the same area as fig. 6 (a) after an attempt had been made to 
reveal the dislocation distribution near the kink band using the etching 
solution referred to in § 2. In the case of the severely strained crystal 
used here, however, the dislocation density was so high that individual 
etch pits could not be distinguished, instead a uniform stain resulted. 
Nevertheless, in regions of extremely high residual stress a localized attack 
did occur giving rise to the ‘cracks’ shown in fig. 6 (b). The interesting 
fact was that each ‘crack’ coincided exactly with the region of stress 
concentration revealed using polarized light in fig. 6 (a). Presumably 
they were the result of a form of stress-corrosion and indicated the position 
and orientation of Stroh-cracks which would have formed next had loading 
been continued. Furthermore, they all lay on a straight line as did the 
Stroh-cracks of figs. 3 and 4, and this line, which presumably delineated 
the length of the associated kink band, was perpendicular to the ‘slip’ lines 
running into the ‘cracks’. 


§ 4. Discussion 


There can be little doubt that the cracks discussed in connection with 
figs. 3 and 4 all form as a result of the piling up of dislocations associated 
with one slip system only. An important question then arises, which of 
the two possible slip systems is actually responsible, the one with its slip 
planes parallel to the cracks in fig. 3 (having a slip direction [110]) or 
perpendicular to them (having a slip direction [110])? In other words, are 
the cracks the result of slip on a number of planes meeting a common 
crystallographic barrier (the kink band) or slip on a single plane meeting a 
number of well-spaced barriers? Although we have no direct evidence 
on this point, the most self-consistent picture evolves when the first of these 
two slip systems is considered responsible. Indeed it is difficult to under- 
stand how slip confined to a single plane could develop a number of inde- 
pendent cracks, since as soon as one crack forms, further shear across the 
plane would preferably enlarge the existing crack rather than generate new 
ones. Assuming the former case to be correct the configuration of the 
cracks and the directions of the slits with respect to the planes of crystalline 
slip represents a direct confirmation of the Stroh analysis. When edge 
dislocations are forced together at the head of a pile-up, the incipient 
crack so formed les perpendicular to the slip plane concerned, this is 
illustrated diagrammatically in fig. 7. The Stroh-cracks of fig. 4 have 
an orientation in exact agreement with this mechanism. 

It is interesting to consider the nature of the barrier and the eventual 
crack development. In our opinion, dislocation pile-ups were caused 
by the kink bands. These kink bands were geometrically necessary for 
the specimen to conform to the grip constraints and were a direct result 
of the non-homogeneous strain. Mott (1951) has considered the develop- 
ment of kink bands within a crystal. Initially dislocations of opposite 
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sign pile against a common barrier. This barrier may be the effect of their 
own stress field interaction or due to imperfections created previous yasa 
result of dislocation interaction on intersecting slip planes. Once a 
localized lattice curvature has started, its influence spreads to neighbouring 
planes causing further dislocation pile-up. The end result is a kinked 
region, shown in fig. 8 (a), across which dislocation pile-ups of opposite 
sign face one another. Further deformation forces more dislocations into 
this region which becomes more severely kinked and thereby more effective 
as a barrier. 


Fig. 7 


Formation of a Stroh-crack by the coalescence of edge dislocations (after 
Stroh 1957). 


Referring to figs. 3 and 4, we have seen that the Stroh-cracks and the 
tips of the secondary cracks are all situated accurately on a straight line. 
This line marks one boundary of the kink band (see fig. 8) which was the 
crystallographic barrier responsible for dislocation pile-up. 

Stroh (1955) has already considered the case of crack nucleation due to 
piled-up groups on several slip planes. He shows that the equation for 
nucleation of a stable microcrack is: 


nboy = 12y 


where o, is the applied shear stress over the slip plane; b is the Burgers 
vector; y is the surface energy; and n is the total number of edge disloca- 
tions involved. It must be emphasized that it is not necessary for the n 
dislocations to be concentrated on a single slip plane; they can be distri- 
buted over a number of different slip planes. In his calculation the shear 
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stress field around each piled-up group acts on the adjacent ones to assist 
the applied shear stress in forcing the leading dislocations against the 
barrier. 

By asimilar analysis it can be shown that pile-ups of opposite signs distri - 
buted as in fig. 8 (a) around a kink band will also assist the applied shear 
stress in nucleating cracks. Indeed the analysis applies for all pile-up 
distributions formed as the result of an applied shear stress of a given sign. 
It is for this reason that the kink band is a preferred site, since fewer dis- 
locations per slip band are required for crack nucleation. 


Fig. 8 
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(a) The distribution of dislocations around a kink band. (4) The position 
and orientation of the Stroh and secondary cracks around the kink band. 


Once the Stroh-cracks have formed perpendicular to their associated 
slip plane, they then run into the compression stress field of the dislocation 
pile-up immediately below (see fig. 8 (b)) and can extend no further in this 
direction. This accounts for the comparatively short length (~5 y) of 
the stable Stroh-crack component. They can, however, spread rapidly 
along the length of the edge dislocation lines, thereby forming a [001] 
direction slit. A slit is most likely to start from the surface where less 
surface energy is required for nucleation. Such was the case for the crack 
illustrated in fig. 6 (a). 

The Stroh-cracks also increase their length by switching through 90° 
and advancing into the region within the kink band, as was demonstrated 
in fig. 3. This can be explained in terms of the normal stresses around the 
pile-ups. In between two pile-ups of opposite sign tensile stresses may 
assist one another and cause the Stroh-crack to run into the kink band to 
form a secondary crack as in fig. 8 (6). There is, however, the additional 
possibility that the cracks may open up due to the residual stresses in the 
erystal when the compressive load is released. 

In terms of these observations, if a single crystal could be deformed 
homogeneously, deformation bands would not form and fracture would be 
delayed. Preliminary observations indicate that this may be the case. 
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Specimens, where particular care was taken in the alignment, were uniax- 
ially compressed by 5%. These specimens were examined carefully for 
surface cracks but none were detected. It is possible, of course, that a 
much more severe uniaxial compression would eventually lead to kinking 
and thence to cracking. Following this line of approach in the future it is 
hoped to learn more of the mechanism whereby kink bands are developed 
and thence cracks are formed in these materials. 


ACKNOWLEDGMENTS 


The authors wish to acknowledge the assistance of Mr. K. H. Olsen 
throughout this investigation. They are also grateful to Dr. J. J. Harwood 
of the O.N.R. and to Dr. F. J. Larsen, Director of Research, for continued 
interest and permission to publish these results. The authors wish to 
thank the Office of Naval Research for sponsoring the work on which this 
paper is based. 

REFERENCES 


Gruman, J. J., 1956, J. appl. Phys., 27, 1262; 1957, Fall Meeting A.I.M.E., 
Chicago. 

GILMAN, J. J., and Jonnston, W. G., 1957, Dislocations and Mechanical 
Properties of Crystals (Lake Placid Conference Report), p. 116. 

Gorum, A. E., Parker, EK. R., and Pasx, J. A., 1957, Paper presented at 
Fall Meeting, American Ceramic Society (Basic Science Div.). 

Lap, R. A., Stearns, C. A., and Det Duca, M. G., 1958, Acta Metall. (to be 
published). 

Mort, N. F., 1948, Engineering, 165, 16 ; 1951, Proc. phys. Soc. Lond. B, 64, 729. 

Nye, J. F., 1949, Proc. roy. Soc. A, 198, 190. 

Scumip, E., and Boas, W., 1935, Kristallplastizitat (Berlin : Springer). 

Sroxgs, R. J., 1957, O.N.R. Report, University of California (to be published). 

Srrou, A. N., 1954, Proc. roy. Soc. A, 223, 404 ; 1955, Ibid., 232, 548 ; 1957, 
Advance Phys., 6, 418. 

WASHBURN, J., 1956, J. of Metals, 8, 189. 


(a7 26am 


Polarization Measurements of Gamma Radiation from 
Sulphur-34; 


By G. G. Sxutst and P. S. Fisner§g 
Department of Physics, University of Birmingham 


[Received April 21, 1958] 


ABSTRACT 


Polarization measurements have been made on the gamma cascade in *4S. 
The experimental results are consistent with a spin sequence 2*, 2+, 0* with 
a mixture of MI and 1-5% of E2 in the first transition. The theoretical 
analysis of polarization—direction correlations with mixing included is 
discussed in some detail. Shell model assignments for the gamma cascade 
are given with possible explanations of the mixing ratio value obtained. 


§ 1. IyTRODUCTION 


THE angular correlation of the gamma cascade in *48 has been studied by 
Handler and Richardson (1956) and by Fisher and Shutej| (1957). Both 
these measurements assigned spins and parities 2*, 2+, 0* to the levels 
shown in fig. 1, the first transition being an M1, E2 mixture. The 
theoretical and experimental values of the A, coefficient (see I) were, 
however, only just within acceptable statistical limits and the polarization 
measurements reported in this paper were made to check the spin sequence 
and mixing previously found. The general techniques used in polarization 
correlation experiments are well known but special problems arise in the 
study of *4S because of the short half-life of the parent °4Cl (323 min) 
and because of the presence of abundant annihilation radiation. 


§ 2. POLARIZATION—DIRECTION CORRELATION MEASUREMENTS 


The polarization correlation for a gamma-ray cascade is given by 
Biedenharn and Rose (1953) in the form 


Wa, d)=, W,(8, b) + 72.W,(8, p). 


Here @ is the angle between the direction of emission of the cascade gamma- 
rays and ¢ is the angle between the electric vector of the gamma-ray 
entering the polarimeter and the normal to the plane contained by the 


+ Communicated by the Authors. 
{1956 Travelling scholar from the University of Melbourne. 


§ Now at Palmer Physical Laboratory, Princeton University, New J ersey. 
|| Referred to as I in the following paper. 
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gamma-rays of the cascade. _», is the efficiency for detecting one gamma- 
ray in the polarimeter and measuring the direction of the other. W,(8, $) 
is the corresponding correlation function. 1, and W,(6, 4) are the corres- 
ponding quantities for the reverse measurement. 


Fig. 1 
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The decay scheme of 34Cl to 348. The shell model assignments to the second 
and first excited states of 34S are included. References to this decay 
scheme are given in I. 

As pointed out by Bishop and Yorba (1955) all the necessary informa- 
tion about the polarization of a selected gamma-ray may be expressed 
by stating the intensities W(47, $7) and W(4z, 0) in the linear polariza- 
tions parallel and perpendicular respectively, to the plane containing the 
cascade gamma-rays and we call p= W/(47, $7)/W(g7, 0) the maximum 
degree of plane polarization. 

To analyse this polarization, the polarimeter relies on Compton scattering. 
Gamma radiation incident on the polarimeter tends to be scattered in the 
plane perpendicular to its electric vector. This can be seen from the 
differential cross section given by the Klein—Nishina formula for scattering 
of an incident quantum of energy ky through an angle 6 with an angle % 
between the direction of polarization of the incident gamma ray and the 
plane of scattering 

do Ck | aK : 
0 - Sole +p 2 sin? 5 cos? |. 

Here & is the energy of the scattered quantum 

ke 


ie, 1+ (ky/mC2)(1 —cos 8) 
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The number of gamma-rays scattered into the horizontal plane, defined 
by the two gamma-rays of the cascade, will be given by 
N,dQ=H dQ= W(}z, 4m) x do,_. + W(}7, 0) x do, _,, 
and similarly the number in the vertical plane will be 
N, dQ=V dQ= W(hn, 4m) x do,_,,7 + W($7, 0) x do, _o. 
The ratio R=do,_,,/do,_, is a measure of the analysing efficiency and 
it follows that H/V=(p+R)/(1+pR). R depends on both the geometry 
of the polarimeter and the energy of the gamma-ray (Metzger and Deutsch 
1950), while » depends only on the theoretical correlation. For quanta 
of energy of 1-2 Mev, R has its maximum value for a scattering angle 
O= 77°. 
§ 3. APPARATUS 
The apparatus for the polarization—direction measurements was based 
upon the angular correlation apparatus described in I. For the polariza- 
tion studies one of the simple sodium iodide detectors was replaced by an 
analyser for plane polarized radiation based on the polarimeters used by 
Metzger and Deutsch and by Bishop and Yorba. 


Vertical Tungsten 
detector alloy shield 
(Nal) 
Plastic 
scatterer 
(NE101) = 


Horizontal 
detector 
(NaI) 


Schematic view of polarimeter showing the arrangement of the detectors and 
the shield. 


The polarimeter is shown in fig. 2. Radiation scattered from the 
plastic scintillator may enter either of the two side detectors (one horizontal 
and the other vertical) which are well shielded by 3 in. of tungsten alloy 
(density 18g cm~*) from the direct radiation. The plastic scintillator 
(NE101), 13 in. in diameter and 2 in. long was attached to a 9514 photo- 
multiplier. The side detectors were standard 14 in. diameter and 1 in. 
long Nal(T1) crystals mounted on 6260 photomultipliers and positioned 
parallel to the scatterer as shown in the figure, 
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To form the horizontal triple coincidence output, H, one gamma ray 
of the cascade was detected by the directional counter and the other was 
scattered in the plastic scintillator. The outputs of these two detectors 
formed the fast-slow system (of resolving time 10-8 sec) discussed in the 
preceding paper (I). The output from this combination was then placed 
in slow coincidence with the single counts from the horizontal detector. 


Block diagram of the electronics. The meaning of the symbols used is as 
follows: A=linear amplifier; D=delay ; D.C.=directional counter ; 
F.C. =fast coincidence ; F.S.=fast-slow coincidence ; H.C.=horizontal 
counter; H.T.=horizontal triples; H.V.=high voltage supply ; 
K=kicksorter; M=ammeter; P=potentiometer ; P.S. = plastic 
scintillator; R=ratemeter; S=scaler; 8.C.=slow coincidence 3 
V=voltmeter ; V.C.=vertical counter ; V.T.=vertical triples. 


P.M. 3 D 
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Similarly the vertical triple coincidences, V, were formed by having the 
fast-slow output in coincidence with the vertical detector independently 
of the horizontal one. The resolving time of all slow units was 10~ sec. 


Fig. 4 
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Spectra of *48 in the plastic scintillator. The top curve shows the output. 
gated by that from a side counter. The bottom curve is ungated. 
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Provision for making random runs was made by switching a delay into 
one of the fast arms. A block diagram of these pulse height selection 
and coincidence circuits is shown in fig. 3. 

Energy selection of the gamma ray detected in the polarimeter was 
achieved by the method of Bishop and Yorba. The spectrum in the 
plastic scintillator was gated by the scattered gamma ray pulses from the 
side counter being considered. As the scattering angle is relatively well 
defined, the scattered quantum has a defined energy and the gated spectrum 
in the plastic resembles that in a Compton spectrometer. Figure 4 
shows both the gated and ungated spectrum of #48 in the plastic. Despite 
the large amount of annihilation radiation present it is possible to select 
the 1-16 and 2-10 Mev radiation from *48 quite well by this method. 


§ 4. OBSERVATIONS 


As in the angular correlation experiment, I, 94"Cl} was formed according 
to the reaction *°S + *He > *4#"Cl+p+n by bombarding a natural sulphur 
target with an external beam of 40 Mev He?* ions from the Nuffield 
cyclotron. The target which was about 2 mm thick, adhered to a copper 
backing and was covered during bombardment by a thin tantalum foil 
(0-0003 in.) cooled by an air blast. Bombardments of 0-5 wa for 20 min. 
provided sources of convenient strength. No additional activities 
comparable in intensity with that of the Cl activity were observed after 
the lapse of a few minutes from the end of bombardment and therefore 
it was possible to use the bombarded sulphur directly as the source without 
chemical extraction of the Cl activity. The sulphur was scraped from the 
target and placed in a thin walled (0-002 in.) brass cyclinder 1/16 in. in 
diameter and 2 in. long, extending from a 2 B.A. threaded brass rod. 

As the half life of °4"Cl is 324 min a new source was prepared every 
half hour and the initial source strength was adjusted to 0-2 mc before 
sealing and placing in the target holder centred between the directional 
counter and the polarimeter. This initial activity was measured by a 
ratemeter attached to one of the side counters. Such a source gave a 
real to random ratio of 6 to 1. 

The pulse height analysers of both the directional counter and the 
polarimeter were set to cover the energy range 0-7 to 2-5 Mev. These 
settings were decided by preliminary runs which indicated that as the 
detection process involved triple coincidences, the overall efficiency was 
so low that both gamma rays had to be included in each detector. This 
does not destroy the advantage of energy selection in the polarimeter 
because the 2 to 1 ratio of the energy of the gamma rays in the cascade 
(see fig. 1) causes the efficiency, ,, for detecting the low energy gamma ray 
in the polarimeter and high energy gamma ray in the directional coun- 
ter to be much greater than that for the reverse measurement, 7p. 
ee eee ey ee ee 

+ The ‘m’ means we are referring to the isomeric level. 
3D2 
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These efficiencies were determined by using a planimeter to measure the 
relative areas of the gamma ray profiles on the spectra obtained in the 
polarimeter and in the directional counter. It was found that , = 979. 

Observations were made in the positions 0=}7 and 6=7. In the first 
p has its maximum value. As the gamma rays do not define any plane 
when @=7, H should be equal to V and so the actual counting rates obtained 
can be used to normalize the @=47 runs. As the mean energy of the 
gamma ray entering the polarimeter is 1/10 x 2-10 + 9/10 x 1-16 or 1-26 Mev 
which is also the mean energy of the gamma rays in the Ni cascade, the 
latter was used for the normalization runs. 

In the 0=4z position 10 000 triple coincidences were recorded at an 
average rate of 2 counts a minute. In the 6=7 position, higher counting 
rates were easily attained with the ®°Ni and 40 000 triple coincidences 
were recorded at an average rate of 15 counts a minute. No changes 
were required to the energy settings for the ®°Niruns. The only changes 
were to the source and the angle setting of the directional counter. 

The result obtained was 

A/V =(H/V),/./(4/V),, = 0-842 + 0-023 
where the quoted error is the standard deviation obtained only from the 
actual counts used. The advantage obtained in using simultaneous 
detection of H and V is that with the easy normalization using ®°Ni, the 
number of cyclotron bombardments needed to obtain this accuracy was 
halved. During the course of the experiment about 100 targets were 
bombarded. 

The polarimeter was calibrated by using the gamma ray cascades in 
SONi and 1%Pd with the following results 


R=1-88+0-12 at 1-25 Mev 
and 
R=2-40+ 0-09 at 0-57 Mev. 


The predictions of the Klein—Nishina formula for point geometry are 
2-40 and 4-55 respectively. A preliminary run on the gamma ray cascade 
in *8Sr gave satisfactory agreement with theory. 


§ 5. THEORETICAL ANALYSIS 


The theoretical analysis of the polarization correlation with the first 
gamma ray mixed has been given by Biedenharn and Rose as 


W (8, é) = W498, ?) as nW (8, ) 
where now W,(6, ¢) is the polarization correlation when the polarization of 
the mixed radiation is measured and W,(6, 6) when the polarization of the 
pure radiation is measured. Now W,(6, ¢) can be reduced to the form 


Wy(8, 6) =e + 8a + Ba 
where 6? is the mixing ratio and w, and w,, are the correlation functions 
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for the pure cascades and w,,, is the interference term. If invariance 
under time reversal is valid, 5 is real. The functions may be written 


y= 1/(20, +1)(— 1yot4 2497, (L,I, L,d,; Tov) 
XZ, (LaF Lod 13 J pv)P, (LL, ; Of) 
yy = 1/(2F, +1)(— 1 pots 2-47 (L' LJ 1; Jv) 


XZ (Lod Lady; J qv)P,(L' L's; Of) 
and 
W441 = 1/(2J,+1)(- Lets SZ (Ly LJ, ; Jor) 


XZ (LJ LJ; Sv)P, (LL ; 94). 
The Z, coefficients have been tabulated in a convenient form by Sharp 
et al. (1954). The function #, is the equivalent for the polarization 
correlation of the Legendre polynomial in standard angular correlation 
theory and is given by Biedenharn and Rose as 
P (LL; 0¢)=P,(cos 8) +(—1)™ 

(v—2)! 2n(v+1)L(£4+1) 

(v+ 2)! v(v+1)—2L(L +1) 
where o(L) is the intrinsic parity of the multipole (0 for electric and 1 for 
magnetic). The function #, for the interference term is given in this 
formulation by 


cos 2¢P,?(cos @) 


2 ne te win (Y—2)! 
PLL’; 6¢)=P,(cos 6)+(-—1) heey 


x (L’—L)(L’ + L +1) cos 2¢P,,?(cos 6). 
For W,(6, 6) where the polarization of the pure radiation is measured, 
the only change to the above is that all the 7, are replaced by 7, (LL, ; 04). 
Thus from these formulae the value of p= W(4z, 47)/W(4z, 0) can be 
calculated for a given spin sequence and mixing parameter (cf. § 2). 


§ 6. COMPARISON WITH EXPERIMENT 


Before the correlation above can be compared with the experiment, 
the corrections for solid angle have to be applied to the theoretical value 
of p. The procedure given by Rose (1953) was considered sufficiently 
accurate for this experiment. The attenuation factors as given by Rose 
were evaluated as Q,=1, Q,=0-89 and Q,=0-70}. To apply these 
factors to the polarization correlation, the associated Legendre poly- 
nomials in this correlation were broken up into a series of Legendre 
polynomials. 

P,?(cos 0) = 2 — 2P,(cos @) 
P,2(cos 0) = 2+ 10P,(cos @) — 12P,(cos 4). 
+ If an angular correlation function is A)+A,P,(cos 4)+A,4P,(cos #) and 


this is attenuated to A,+A’,P,(cos 0)+A’,P,(cos 8) then Q=1, Q,.=A's/A2 
and Or AIA: 
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With the small but significant efficiency of each detector for the. 
abundant annihilation radiation resulting from the annihilation of the 
positrons there is obviously a small isotropic contribution to the correlation. 
This was determined experimentally as in I. The efficiency of the 
polarimeter and directional counter for this radiation was measured by 
setting one on the annihilation peak and leaving the other unaltered. 
The rise in counting rate at 9=7 then gives the efficiency of the unaltered 
detector. The isotropic component was estimated as 16+ 1% in W(8, 4). 


Table 1 
J, a 6 5’ (H/V)' a) H/V | Mixture type 
1+ 1+ Ot —2:094 | 1:37.| —0-477. | 1:37 Ml... E2 
I es 0° —2-094 | 1:42 | —0-477 | 1-44 Ml E2 
Va I= Os —2:094 | 0-73 | —0-477 | 0-73 El M2 
Le as Or —2-094 | 0-70 | —0-477 | 0-69 El M2 
a+ 1+ 0+ 7:73 0-82 0-451 | 0-97 Mi £E2 
2 1 0 7:73 0-87 0-451 | 1-03 Mi £E2 
2+ i= OF 7:73 }-22 0-451 | 1-038 El M2 
DiS 1+ Or Tis 1-15 0-451 | 0-97 El M2 
3+ leg 0+ — 13-12 0-83 0-710 | 1-03 M2 ~~ E38 
37 1+ Ot — 13-12 0-78 0-710 | 0-97 M2. E3 
3t 1% Ot —13-12 1-21 0-710 | 0-97 12. M3 
3 i= O+ | —13-12 1-28 0:710 | 1-03 E2>~* M3 
a+ at 0+ 0-126 | 0-82 — — Mi —E2 
2a a7 0+ 0-126 | 0-87 —_ — Mi-_—-— «E2 
2+ a= Ot 0-126 | 1-22 — — El M2 
27 2+ Or 0-126 |) 1-15 — — El M2 
3+ 2+ 0+ —0-767 | 1-19 i a Ml. E2 
oe 2 0+ —0-767 | 1:24 — _ Ml E2 
3+ 2- Ot —0-767 | 0-84 — — E1 M2 
ab a+ o+ —0-767 | 0-81 — — E1 M2 
3 3 Wy) —0-055 | 0-72 = — MI -. K2 
3+ 3t O+ —0-055 | 0-75 — — Ml &E2 
37 3t 0+ —0-055 | 0-31 os a= El M2 
3t 3- O+ —0-055 | 0-36 — — El M2 


| 


This table lists all the levels schemes consistent with the angular correlation 


experiment (I). The value of H/V listed is to be compared to the experimental 
value of 0:842 + 0-023. 


In I the spins of the states of #48 allowed by the directional correlation 
experiment are tabulated, together with suggested values of the mixing 
parameter, 6 for the first transition. The corresponding values expected 
for the ratio H/V (cf. §2) are found from the theoretical prediction of v) 
corrected as above and from the experimentally determined R (cf. § 4). 
These values are given for all possible parity combinations with the 
selected spin sequences in table 1. The cascades which seem consistent 
with the observed H/V of 0-842 + 0-023, i.e. which lie between the values 
0-89 and 0-80 are in heavy type in table 1. To obtain the values for the 
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cascades 3, 2, 0 which gave no real solutions for 5 when the value of 
A,=0-332 was used, the value A,= 0-300 which is one standard deviation 
away, was used. The mixing ratio used for 2, 2, 0 is a mean of that 
obtained by Handler and Richardson (1-:7°%%) and in I (1:4%). The 
mixture considered in the table is always that between the lowest order 
multipoles allowed by the spin and parity values. 


§ 7. Discussion 


In the arguments presented in I for assigning the spin sequence 
2*, 2*, O* to the levels of *48S, it was assumed that the beta decay of 
3amCl(3+) to the first and second excited state of 348 was allowed. The 
ji-value for the decay to the first excited state is however just consistent 
with first forbidden decay and the schemes allowed by the polarization 
experiment have been examined with this possibility in view as well. 


Table 2 

hee i gid, ely Beta decay to J, Beta decay to J, ‘ 
ee Lee Or: Allowed 2nd forbidden 
ea tesa (3: Ist forbidden Ist forbidden 

2 (team wigs Allowed Allowed 

Fe as 0 Ist forbidden Ist forbidden 
Bee UI es Allowed lst forbidden 
Sue ea Oe Allowed. Ist forbidden 
See ee Oy Ist forbidden Allowed 


In this table the beta decay from the 3+ level in 34”Cl to each level has been 
classified for each parity and spin assignment allowed by this experiment. 
Experimentally, it is known that the decay to J, is allowed and that to J, is 
either allowed or Ist forbidden. 


The evidence is summarized in table 2. It is known experimentally 
(see fig. 1) that the ratio of the intensity of the crossover to the first gamma 
ray is 7/15. In table 3 are presented the theoretical values (Moszkowsk1) 
for the mixing ratio and the above intensity ratio as well as the 
experimental values. 

It seems reasonable to assume on the evidence presented in tables 
2 and 3 that the only acceptable cascade is 2+, 2*, 0*. This sequence 
requires only the assumption of slight mixing (1}+3)% to bring 
experimental and theoretical correlations into agreement. Both the 
mixing and relative intensity of the crossover seem favourable. Beta 
decay to the 2* levels from the 3+ level in *4”"Cl would be an allowed 
Gamow-—Teller transition agreeing with the ft values. Both the 2* assign- 
ments are in reasonable agreement with even—-even nuclei systematics. 
In view of this evidence we must conclude that the slight difference 
between the observed and calculated values of A, (see 1) is not significant. 
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Morinaga and Bleuler (1956) have assigned the jj-coupling shell model 
schemes to the levels in 348, as shown on fig. 1, from the beta decay rates 
to these levels from #4P and %4Cl. The scheme given applies to the 2 
protons and 4 neutrons outside the djj. shell which closes at *8Si, and 
these particles will then be in the 8,2 or dy. state. The second excited 
state, and the 3+ isomeric level in #4Cl are assumed to have either four 
S12 and two dg particles or six dy, particles, whereas the first excited 
state and the 34P ground state, should consist of one to three s,,. particles 
and five to three dy. particles. If the M1 transition between the second 
and first excited state were due solely to the transition of a single particle, 
then as it would be dsj>81)2 OF 8yj2>dg/. it would be * 1-forbidden’ (Mayer 
and Jensen 1955, or Moszkowski 1955). 


Table 3 


Mixture | Theoretical | Cross-| Theoretical 
Ji te A 6° type mixing over ratio 


2 eke Or 6-0 x 101 Ml E2 76x 10-4 K2 Lb6éx10-* 
mie mee oe G0 oe LO* Ml K2 Cie an {age M2 EXO 
2 Ziel) 6x 10-2 | "Mie 2 TOO * E2 L-6x 207! 
nis APR A ir 6x 10>") MT 2 76x 10-4 M2 iio 
Soe be, BO 7x 10-4) |e M2 ES Z1oclO4 M3 3:4 x 10-2 
Re eee Ong >OxX10=) EP M2 1-0 10-% M3 ao 107 
en Cy it whley 5:9x10-* | EP «M2 Oc Les K3 19> 1057 


In this table the experimental mixing ratio 8? (the ratio of the intensity of 
the higher order multipole to the lower) is given in column 2. This is to be 
compared to the theoretical single particle estimate in column 4. Column 6 
gives the theoretical value for the ratio of the crossover intensity to that of the 
lower order multipole. The experimental value for this ratio is 7/15 (see fig. 1). 


As the mixing ratio is 20 times larger than the single particle estimate 
and the crossover intensity is 3 times larger than expected on this model 
for 2+ 2+ 0* (see table 3), we find that we can explain these values as shell 
model phenomena by introducing the 1-forbiddeness as discussed above. 
The alternative description furnished by the Bohr and Mottelson model 
is that there will be collective enhancement of the E2 radiation from the 
second excited state. The values quoted above seem reasonable either 
as an estimate of the reduction in the dipole matrix element or as an 
estimate of the collective quadrupole matrix element. On the whole, as 
the evidence for the shell modei classification of the states seems quite 
sound it would seem as though the mixing experienced in #48 is a shell 
model phenomenon rather than a collective one. 
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ABSTRACT 


Experimental results obtained by studying the low temperature creep 
behaviour before and after a small stress change are interpreted in terms 
of barriers due to intersecting dislocations. Using the data available on 
the temperature dependence of the flow stress in conjunction with these 
experimental results, estimates of jog energies are derived and compared 
with values obtained from flow stress data alone. Relative dislocation 
widths for various metals are also deduced. From the theoretical analysis 
relating dislocation width to stacking fault energy, tentative values are 
suggested for the stacking fault energies of Au, Ag, Al, Ni, Pb, Zn, Cd, Co and 
Pt, relative to that of Cu. These values are in reasonable agreement with 
previous estimates where available, thus confirming the hypothesis that 
low temperature creep is controlled by the proposed mechanism. 

The significance of these results is discussed in terms of the exhaustion 
and work hardening theories of low temperature creep. It is shown that 
the experimental observation that the apparent activation energy is 
insensitive to stress and pre-strain can be explained quite naturally on 
theories based on stress-aided activation and that it does not invalidate 
such theories. 


§ 1. INTRODUCTION 


THE temperature dependence of the flow stress and transient creep in 
f.c.c. and h.p.c. metal crystals have been interpreted in terms of a process 
of dislocations cutting through a ‘forest’ of intersecting dislocations 
(Cottrell 1953, Seeger 1954, Friedel 1956). Seeger has developed an 
expression which can be applied to both the flow stress and the creep data, 
and which involves the jog energy and the width of a dislocation as para- 
meters. The present paper describes experiments designed to determine 
the jog energies and relative dislocation widths in a number of metals. 
Values of stacking fault energies can be deduced from these results, and the 
general agreement obtained between these results and previous estimates 
confirms the correctness of the proposed mechanism. 


§ 2. THEORY 
The application of rate process theory to a dislocation model for creep 
gives the creep rate ¢€ as: 


€= NAbyyexp (— U/kT) er tA se 0 


* Communicated by the Authors. 
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where NV is the number of dislocation elements per unit volume contri- 
buting to creep, A the area swept out by a dislocation element. after 
thermal activation, b the Burgers vector, vy an atomic frequency and U 
an apparent activation energy. Experimental studies of the variation of 
creep rate with time indicate that € may be expressed as: 


oe CE Cee aE yas et (2) 


with the parameters C and n depending on the test temperature, the 
applied stress, the defects present in the metal and their distribution. It 
has been shown experimentally that for many polycrystalline metals the 
form of (2) applicable to low temperatures is: 


Ce ee Se > Re CO) 


(Smith 1948, Wyatt 1953, Olds 1954). In the region where (3) is obeyed, 
i.e. in logarithmic creep, it follows from (1) and (3) that: 


Ges cP loo ee), oy one 4) 
where v= NAb», /C. 


In recent years U has been measured directly by comparing creep rates 
before and after a small abrupt change in temperature (Sherby e¢ al. 1957, 
Landon et al. 1957, Lytton et al. 1957, + Haeszner and Krause 1957). U 
is found to be insensitive to stress and pre-strain at a particular temper- 
ature. The significance of these results will be discussed in §5, but it 
should be stated in passing that the present authors believe that no 
experimental evidence exists which conflicts with theories of low temper- 
ature creep based upon stress-aided activation. Before U can be expressed 
in terms of the applied stress a realistic model for the barriers impeding 
flow in low temperature creep must be chosen. 

Experimental studies of the flow stress as a function of temperature 
have been interpreted, in the low temperature region, in terms of barriers 
due to intersecting dislocations (Cottrell 1953, Seeger 1954, Friedel 1956). 
The barriers operative in low temperature creep are also thought to be 
these ‘forest’ dislocations. In f.c.c. and h.c.p. metals the dislocations 
lying in close packed planes are dissociated into partial dislocations with 
a stacking fault ribbon between them. (Cottrell 1953, Read 1953). 
Before such dislocations can cross each other the stacking fault ribbon 
must become constricted at the point of intersection. Part of the energy 
U, required for the process of intersection is due to the energy of the con- 
strictions and part is due to the energy of jogs or kinks produced as a 
result of the intersection. In the presence of an opposing internal stress 
and an applied stress the activation energy, U, will be equal to Uy minus 


+ These papers by Dorn and his coworkers will be referred to in future as 
Dorn et al. (1957 a, b and c). 
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the work done by the difference of the applied and internal stresses in 
moving the dislocations across each other, so that 


U=U,-1,.6.D.(r—7,). eter ee ad (irae 25). 


where 1, is the average distance between ‘forest’ dislocations, thy the 
activation distance, i.e. the distance travelled by the moving dislocation 
in surmounting the barrier, 7 the average resolved shear stress in the slip 
direction, 7, the average opposing internal shear stress resolved in the slip 
direction and v (=1,bD) is called the activation volume. ; 
In the absence of a more complete analysis of the crossing mechanism 
the process of intersection is thought to be as shown in fig. 1; fig. 1 (a) 


Fig. 1 


(a) 


Intersection of two dislocations showing the increasing degree of constriction 
at the intersection as two extended dislocations are forced together by 
an applied stress. 


shows the dislocation ribbons when they are just touching, figs. 1(b) and 
(c) show the configurations obtained by increasing the applied stress. 
Consideration of such diagrams indicates that 


D~d,+ds, 


where d, and d, are the effective widths of the two dislocations (see 
Appendix (A)). 
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From (4) and (5) we obtain: 
kT logvi=Uy—-v(r—t,). « . ww te (6) 


Suppose that the applied stress 7 is changed by a small amount Ar at a 
constant temperature; then 


kT A log vt = —vAr. ae ne eee ee ees) 


provided the stress change is small, so that v and +, may be regarded as 
constant. Multiplying (7) by 7 gives 

eerie tNUOS IAT. “dos oe (8a) 
If the applied tensile stress is o, the fractional charge in the resolved 
shear stress (Ar/7) may be equated to the fractional change in the applied 
tensile stress (Ag/c). Equation (8a) becomes 

ne RTA (log pt)o) Aon) eo oe es 8B) 
Hence by changing the applied stress slightly during a creep test, the 
quantity vr can be studied. Experiments designed to give relative 
values of D for different metals and to test the basic mechanism of logarith- 
mic creep are described in the following section. Such experiments were 
first performed by Wyatt (1953) on polycrystalline copper and discussed 
by Cottrell (1953) in terms of barriers due to intersecting dislocations. 


§ 3. INCREMENT TESTS ON POLYCRYSTALLINE PURE METALS 
3.1. Experimental Method 

Bright-drawn wires of pure metals (Johnson, Matthey) were heat- 
treated to obtain recrystallization with a small grain size. Electropolished 
creep specimens of diameter approximately 0-5 mm and gauge length 5cm 
were prepared and loaded into self-aligning friction grips on a constant 
stress creep machine. The stress constancy was obtained by use of a 
linear spring. Apart from the controlled change in stress, the specimens 
were crept under conditions of constant stress. The strain sensitivity 
was 2x 10-5. Temperature control, in these short time tests, was obtained 
by use of constant temperature baths; liquid nitrogen, solid CO, in ace- 
tone and ice/water mixture. A stirred water bath in a large Dewar vessel 
was satisfactory at room temperature. The temperature remained 
constant to within +0-25°c for several hours. Figure 2 indicates the 
experimental procedure. 

Each specimen was crept at a low stress o for five to ten minutes and 
the creep curve recorded photographically. A small change Ac was made 
and the creep behaviour recorded for a further ten minutes. This pro- 
cedure was repeated at a series of increasing stresses throughout the range 
of logarithmic creep. Each specimen gave a value of v7 at each stress 
level used. Successive specimens were crept at the same stresses but with 
increments of different magnitude. Alogvt was obtained from a plot of 
strain against log time on which the curve subsequent to the increment is 
moved along the axis of logt until its slope is equal to that of the initial 
curve (see fig. 3). 
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3.2. Hauperimental Results 


Figure 4 shows the experimental results for Au, while table 1 sum- 
marises the results obtained to date. 

In each case v7 was insensitive to stress, pre-strain and, with one 
exception, independent of the size of increment or decrement used. This 
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Illustration of the creep tests performed on one specimen. 
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Method of obtaining J log t. 
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exception was Ni in which the value of v7 increased significantly with 
increment at the lower stresses used. The reason for this difference is 
not understood. In this case the curve of v7 against Ac/o was extrapolated 
to Ac/o=0 and the value of v7 so obtained quoted in table 1. In the case 
of Cu and Au where measurements were carried out at several temperatures 
a dependence of v7 on temperature was found. 


Table 1. Values of v7 from Creep Increment Tests on 
Polycrystalline Pure Metals 


meta | TEP | angen) | Steers | specimen 
Cu 300 5-0+0-5 6 to 18 10 
7 Cu 80 4-2+40-4 8 to 20 4 
Au 300 Le 4-4+0-3 5 to 12 10 
Au 80 . 2:°7+0°3 5 to 12 2 
Pt 300 9-0+1-0 8 to 15 2 
Ni 300 6-7+0-6 14 to 28 9 
Co 300 5640-5 35 to 70 9 
Ag eC 300 4-5+0-5 6 to 13 10 
Zn 80 1:0+0-2 3 to 6 e 
Al 80 0:9 + 0-2 2 to 5 2 
Pb 80 1:340-5 1 to 3 I 
Cdt iS 80 1-4+40:3 single cryst. 1 


7=0-3 


7 Taken from Glen (1956) 


§ 4. Discussion or RESULTS 
4.1. Dislocation Distributions 
Equation (6) may be written 
Uy+vry=07r+ U sur + kT log vt. ed ahd We) 
The results of Dorn and his coworkers show that U is insensitive to stress 
and pre-strain (see $5 for full discussion) and the present results indicate 


that vz is also independent of stress and pre-strain. It follows from (9) 
that v7, is likewise independent of strain, since U, is a constant. 
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We may write 
T= aGb/l, Se tn ee (LO) 


where « is a constant which may depend on the geometry of the dislocation 
arrangement, /, is a dislocation length or separation corresponding to the 
dislocations responsible for 7,, and @ is a shear modulus. Hence 


vty =bDlyr,=0Gb2DI|l,. . . . . . . (Al) 


On current theories of work-hardening of f.c.c. metals (Seeger et al. 1957), 
J, is considered to be the average distance between dislocations or groups 
of dislocations on the primary slip system. (These dislocations will be 
referred to as the ‘primary ° dislocations.) This distance J, is therefore 
not necessarily the same as /;, the distance between ‘forest’ dislocations. 
It appears from the strain independence of v7, however that al,/l, is a 
constant. A similar result has also been deduced by previous workers 
from flow stress measurements (see below). The interpretation of this 
result is not clear. If the dislocation arrangement remains similar during 
the deformation, but changes only in scale, « will be constant, and on this 
interpretation the constancy of al,/l, implies that the density of ‘forest’ 
dislocations increases proportionately to that of the ‘ primary ’ dislocations 
with increasing deformation. (This will be referred to as the ‘ Cottrell— 
Stokes’ law.) It is possible however that 7, represents the long range 
stress due to the same ‘ forest’ dislocations as are responsible for 1, in which 
case 1, may be equal tol,. In the absence of a proper analysis of the inter- 
action between intersecting dislocations it is not possible to decide between 
these and other possibilities. However, whatever the correct explanation 
may be, since al,//, is constant, we may write 


eee ey ont ee eoeemiae. Bes 91(E) 


where a, is a constant (for further discussion see § 6). 


4.2. Relationship of the Present Results to Measurements of the 
Temperature Dependence of the Flow Stress 


The flow stress of a metal is considered to be the sum of two contri- 
butions: a temperature independent part due to the internal stress 7, and 
a temperature dependent part due to ‘forest’ dislocations. Using 
Seeger’s theory of the temperature dependence of the flow stress (Seeger 
1954) and neglecting Friedel’s modification (Friedel 1956), the flow stress 
is given by 


t/7g=1+ U,/vrg+(kT[ur,)In(ée/NAby) . . . (18a) 
for T'<T,, where T= U,/k log (é/NAbv)) and 
7=7, for ite = Se (13 b) 


This is indicated in fig. 5. 


P.M. 35 
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By measuring the relative flow stress at a series of temperatures Adams 
and Cottrell (1955) and Cottrell and Stokes (1955) showed that in single 
crystals of both copper and aluminium 7/7, is independent of the strain. 
It follows from (13a) that v7, is constant, as pointed out above. Cottrell 
and Stokes (1955) also found the same value of 7/7, for polycrystals as for 
single crystals. In the following sections it will be assumed that this 


result is true for all other metals. 


4.3. The Energy for Jog Formation 


From the temperature dependence of the flow stress it is possible to. 
obtain values of U,/v7,; putting 7’=0°K in (13a) leads to the expression : 


To[Tg=1+ Uy/e7,. ree ee 
By rearranging (6) we obtain 
vr (U forg+l)=ortkTlogvt. . . . . ~ (15) 


The frequency term v can be determined for Al, Cu and Ni from the results 
of Dorn ef al. (see §5). The value obtained (10%sect where w= 10-5 + 1-5) 
is roughly the same for all three metals. This value is assumed for the 


Fig. 5 


TEMPERATURE 


Temperature dependence of the flow stress shown schematically. 


other metals studied ; it is unlikely to be greatly in error and since kT log vt 
is small compared to v7 the uncertainty in v is not important here. Using 
(14) and (15) values of both Uy and vr, can be derived. Unfortunately 
flow stress measurements are not yet available for all the metals studied. 
The available values are listed in table 2. Values derived from sparse 
data are indicated by an asterisk. The values of U,/vr, for the low 
stacking fault energy metals (Cu, Au and probably Ni) ane somewhat 
uncertain for two reasons: (1) 7, is taken to be the value of 7 in the 
temperature range 400 to 600°K, where apparently 7 has little, if any, depen- 
dence on temperature ; however, the true value of 7, may be a little lower 
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(Seeger 1955a). (2) The curve of 7 vs 7 is more complex for Cu than is 
shown in fig. 5 (Adams and Cottrell 1955); an average extrapolation of 
this curve to 7’'=0°K gives the value quoted in table 2. The possible 
errors in the values of both vr and of U)/v7, lead to uncertainties in the 
values of U, and vr,, indicated in tables 3 and 4. In table 3 values of U ‘ 


Table 2. Values of (1+ U /v7,) from Flow Stress Measurements 


Metal 14+ U,/e7, Source 
Cu 1-2+0-08 Adams and Cottrell (1955) 
he 1:35+0-1 Seeger (1954), data of 


Andrade and Henderson 


; end Seeger (1954), data of 
Ni Eee U lS Andrade and Henderson 


Ni 1-20+ 0-05 Haasen (1958)+ 


Seeger (1954), data of 


: 2 
au seat? Boas and Schmid 
Lander and Howardt 
Dp 2. : 
Pb 2:-44+0-4 (1936) 
Seeger (1954), data of 
2. ae fo) 2 
Bs poate Schmid and Fahrenhorst 
1-40 Cottrell and Stokes (1955) | 
= 2.5 Seeger (1954), data oft 
Rosi and Mathewson 
Co (1-2) Assumed 
Pt (1-35) Assumed 
Ag (1-2) Assumed 


+ Data kindly supplied by Dr. Haasen prior to publication. 
t Limited data. 


are compared and found to be in reasonable agreement with approximate 
values given by Friedel (1956), whose values were derived from (13a) by 
making plausible assumptions for Naby, and inserting the experimentally 
determined value of 7. A further comparison may be made with a 
theoretical estimate as follows. A study of the crossing process shows that 
in addition to two constrictions, two lengths ~b, and b, of dislocation line 
must be formed (b, and b,) are the Burgers vectors of the intersecting 


352 
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dislocations). These will be kinks or jogs depending on the orientation and 
Burgers vectors of the dislocations considered (Cottrell 1953). The total 
increase in dislocation length is, to a first approximation, the same for all 
possible orientations and Burgers vectors. The total energy will therefore 
vary with orientation only through the energy of formation of the 
constrictions, which may depend somewhat strongly on orientation. 
Thus the values of U, derived above probably represent the energy 
required to form constrictions in two dislocations of average orientation, 
plus the energy to form, for example, two jogs, i.e. Uy is nearly twice 
the energy to form a constricted jog in a dislocation of average 
orientation. Thus the experimental values quoted in table 3 can be 
compared directly with the theoretical values obtained by Friedel 
(Friedel 1956). 


Table 3. Estimates of U, obtained from: (i) Combination of Creep and 
Flow Stress Data, (ii) Flow Stress Data alone, (ii) Friedel’s Calculation 


Us (ev) Uy (ev) T (ey 
Metal ee : From creep and | From flow ae fe 1 
: flow stress data | stress data apg Ye 
300 1:0+0:5 0-7 2-6 
Cu aise Si femtomash nba > : 
80 0-7 + 0-35 0-7 1-9 
300 1:3+0-5 0-8 2-1 
Au — 
80 0-8 + 0:35 0:8 1-1 
Ni 300 1-5+0-6 0-8 2:9 
Cd 80 1-140:3 0-8 0-3 
Zn 80 0-7+ 0-25 0-6 0-3 
Al 80 0-5 + 0:35 0-3 0-2 
Pb SO ; 0-85 + 0:50 0-3 0-25 


If the values of the dislocation widths given in § 4.4 of this paper are used 
in Friedel’s equation, the estimates listed in the final column of table 3 
are obtained. In general these calculated values are of the same order of 
magnitude as those determined experimentally. It should be emphasized 
that these experimental values of the jog energy 4U, are of the order of 
+ to lev. ' 

4.4, Interpretation of v7, 


Using the values of tables 1 and 2 estimates of v7, are obtained and are 
shown in table 4. For Co, Pt and Ag values of 1+ U, /v7, are not available. 
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Comparison of the plastic behaviour of these metals with Cu, Ni and Au 
suggests the values which are given in brackets in table 2 and which will 
be assumed in the subsequent discussion. In § 4.1, it was shown that 
vr,=a,Gb°D=a,Gb(d,+d.) . . . . . . (16) 
where d, and d, are the effective widths of the two dislocations. In 
Appendix (A) it is shown that the volume average (D) of the activation 


Table 4. Values of (a) vr, obtained from Creep and Flow Stress Data 
(6) Relative Dislocation Widths (1) and Estimates of Stacking-Fault 


Energy (y) 


Most probable values 


ee a ad ene 
(ergs/cm?) 

Cu 300 4-95+ 0-8 4-5 8-57 407 

Cu 80 | 3540-55 aes 6-1+ ‘Bey 45 <57 < 80 

Au ~ 300 3°85 + 0-6 2-7 7942-4 - 24<33<47 

Au SO | 2-154+0-3 2-8 3-941-2 50 < 65 < 90 
Ni 2 300 | 5:6 +0-7 7:8 6342-0 64<95< 140 

Cd 8O | 0-4140-11 2:0 11 0-5 > 130 

Pb 80 | 0-6 £03 0-75 2541-3 25 <50< 140 

Zn sero ® 0-48 + 0-14 3:9 1-0+0-45 > 250 : 

Al oa: 0-61 + 0-25 2 i 0-74.04 > 200 
prose 300 b3 + 1-0 5:20 18-5 + 6-9 22 < 26 < 36 

Pt * 300 | 7-2541-4 6-2 16 t2e7, 56 <69<110 
ae 300 3°8. + 0-8 2-9 72+ 2:6 rj 26<35<58 on 


+ Assumed. The uncertainties quoted represent the probable experimental 
error ; no estimate of the uncertainty due to the approximation used in the 
analysis relating dislocation width and stacking fault energy has been included 
in the limits shown for y. 


distance can be expressed in terms of the width (d) of a dislocation line 
which makes an angle of 45° with its Burgers vector. Assuming that in 
f.c.c. metals both dislocations lie in close-packed planes and are therefore 


extended, “¢ 
Dee aod 2d. 
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In h.c.p. metals the forest dislocations are assumed to lie in pyramidal and 
prismatic planes and are assumed to be unextended, whence 


Diop, =O. 
Defining N =d/b we obtain from eqn. (16) 
v7, {Gb'=a,2N (f.c.c.metals) . . . . . (17 a) 
while v7 (Gb3=a,(N+1) (h.c.p. metals) . . . (17 b) 


In $4.1 it was shown experimentally that «, is a constant independent 
of pre-strain. The constancy of vr, with respect to pre-strain indicated 
that «, is independent at least of the scale of the dislocation distribution. 
If the assumption is now made that «, is independent of the geometry of 
the dislocation distribution and of the temperature, relative values of V 
for different metals can be calculated from the given experimental data. 
Such values of NV are shown in table 4, where a value of NV =8-5 has been 
assumed for Cu at room temperaturet. This assumption is based on 
calculations due to Seeger and Schoeck (1953) who showed that the widths 
may be considerably different for dislocations in the screw and edge 
orientations. Assuming a value of 40ergs/em? for the stacking fault 
energy, these authors obtain for Cu Neage=12 and Negcrey=5. If these 
values are accepted, a reasonable estimate for a dislocation of average 
orientation in Cu is N=8-5; while no great accuracy can be claimed for 
dislocation widths derived in this manner, these values are generally in 
agreement with estimates derived from theoretical considerations, work- 
hardening characteristics and twin boundary energies. This agreement 
lends support to the view that the assumed mechanism for logarithmic 
creep is correct. 

It is of interest to calculate the magnitude of «, from eqn. (17a). In- 
serting the values applicable to Cu at room temperature, «, is found to be 
approximately 0-1. If 7, is due to the internal stress of the ‘primary’ 
dislocations, it would seem reasonable to expect « (in eqn. (10)) to be about 
1/27; hence it follows from eqn. (12) that l,~/,. Thus it appears on this 
interpretation that in the regions where the activation is occurring the 
densities of dislocations on the primary and secondary system are approxi- 
mately equal. If the dislocations on the primary system occur in piled- 
up groups, as expected for example on Seeger’s theory of work-hardening 
(Seeger et al. 1957), lL, ~nl;, where l,, is the spacing between groups and n 
is the number of dislocations per pile-up. In this case the density of dis- 
locations of the secondary system is ~n times that of the primary system, 
a rather unlikely result. It seems that there is at present no satisfactory 
explanation of this experimental result (see however § 6). 


4.5 Stacking Fault Energies 


Estimates of the stacking fault energies (y) can be obtained from the 
observed values of the dislocation widths. 


+ In these calculations we have used average shear moduli for f.c.c. metals 
and 1/s,, for hexagonal metals. 
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Seeger and Schoeck (1953) and Seeger (1955 a) have carried out detailed 
calculations of dislocation widths as a function of the stacking fault energy 
and elastic constants. The calculations are based on the Peierls model of 
a dislocation and take the crystal anisotropy fully into account. Full 
numerical calculations have been completed only for Cu, Co, Ni and Al; 
however estimates of the dislocation widths as a function of stacking 
fault energy in the other metals studied can be obtained from the analysis 
of Seeger and Schoeck by an approximate method. The approximation 
involved concerns the parameter y (see Seeger and Schoeck) for which the 
isotropic value is taken, instead of the value obtained by full consideration 
of the anisotropy. An estimate of the error introduced by this approxi- 
mation was made by comparing the values obtained by the approximate 
method for Co, Cu, Niand Al with those obtained by Seeger and Schoeck 
using the full analysis. The maximum difference was found} to be 20%. 
In table 4 we have shown for each metal the value of the stacking fault 
energy which has to be assumed to obtain agreement between experimental 
and calculated estimates of the dislocation widths. Taking into account 
the errors in the experiments (which are reflected in the limits of stacking 
fault energy given in table 4) and the uncertainties in the calculations, 
these values can be considered to be in satisfactory agreement with previous 
estimates where available. Thus Cu, Ag and Au all appear to have 
similar low stacking fault energies, as expected from theoretical con- 
siderations (Seeger 1955a, b). 

The value of 95 ergs/cm? for Ni is in good agreement with the value of 
90 ergs/em? estimated by Haasen (1958) from an interpretation of the 
temperature dependence of the onset of stage IIT of the work-hardening 
curve for single crystals. The comparatively small (but rather inaccurately 
determined) value for Pb was at first sight surprising, but is confirmed by 
the fact that observations in this laboratory indicate that annealing twins 
occur rather more frequently in Pb than in Al, which is a typical metal of 
high stacking fault energy. For Al, Zn and Cd only lower limits could be 
determined, but it is clear that these metals have high stacking fault 
energies. The value for Co appears to be rather higher than expected 
from data on the martensitic transformation (Seeger 1955a), but it is 
quite possible that the constant «,, assumed to be the same for all metals, 
may in fact have a different’ value for the hexagonal metals. 

A further consequence of assuming «, to be independent of temperature 
is the apparent dependence of dislocation width on temperature in Au and 
Cu (see table 4). It is possible that this dependence is due to a significant 
entropy contribution to the free energy of an extended dislocation. This 
contribution has hitherto been neglected. Both the entropy due to the 
stacking fault and that due to the partial dislocations should be taken into 


+ In these calculations values of the shear moduli called Gy, by Seeger and 
Schoeck have been used, except for Pt for which the value for the polycrystal- 
line metal had to be assumed. These values are not the same as those shown 
in table 4, which are used to derive N from v7,. 
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account. Assuming isotropic elasticity, the approximate width of an 
extended dislocation is given by 


go," 4 8a 

d= Lr {mo-(Sr+ at) pete Mean 
where K’ is a function of Poisson’s ratio, 6, is the Burgers vector of the 
partial dislocation, y, and S, are the energy and entropy of a stacking 
fault, and @S,,/dd represents the rate at which the entropy of the twe 
partial dislocations (S,) varies with their separation (d). It is not clear at 
present which of the two entropy contributions would be the more im- 
portant. However, it is expected that for large d, S, is important, while 
for very small d, 0S,/dd may make the larger contribution. Since the 
Debye temperatures for both Au and Cu lie between the temperatures at 
which creep measurements were made, (S,+0S,/éd) is probably tem- 
perature dependent. However, the average value of (S,+0S,/dd) 
between 80 and 300°K is found to be 0-114ergs/em?/°K for Au and 
0:0544 ergs/em?/°K for Cu. These values are somewhat larger than the 
value of S, for Co of 0-0228 ergs/cm?/°kK estimated from the martensitic 
transformation. It is hoped that experiments at present being carried 
out on the temperature dependence of vz, for Co will lead to a clarification 
of this point. 


§ 5. THE RELATIONSHIP BETWEEN THE PROPOSED INTERPRETATION 
AND ForRMAL THEORIES 


The discussion in the preceding sections deals with the interpretation of 
the activation energy for logarithmic creep in terms of the ‘forest’ har- 
dening. This interpretation taken by itself does not lead to a prediction 
of the dependence of creep rate on time. Two formal theories have been 
put forward in explanation of logarithmic creep. (a) The exhaustion 
theory (Mott and Nabarro 1948, Smith 1948, Cottrell 1953, Mott 1953); 
(6) The ‘ work-hardening’ theory, (Orowan 1947, Mott 1953, Seeger 1954). 

On the former theory, the metal is thought to contain dislocation barriers 
of different heights and during the course of creep, the dislocation elements 
hindered by low barriers move first. If no replenishment occurs they 
become exhausted with a resultant increase in activation energy with time. 
The theory leads to logarithmic creep with an activation energy given by 


U=KkT log vot 


where vy is the atomic frequency. 


On the work hardening theory, the activation energy is taken as pro- 
portional to the creep strain, whence 


U=kT log vt 
with 
v= vo( Nb? BO)/kT a, 
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where NV is the number of elements per unit volume, / the average distance 
between forest dislocations, B the jog energy, @ the strain hardening co- 
efficient and o, the applied stress. 

The two theories are similar in the sense that for both there is an acti- 
vation energy characteristic of any given time during the creep curve, and 
that this activation energy increases gradually with time, in one case be- 
cause the elements of lower activation energy become exhausted, in the 
other because the internal stress increases gradually and therefore the 
activation energy, which is assumed stress dependent, also increases with 
time. 

On the exhaustion theory, the activation energy is independent of stress 
and pre-strain and since vy~ 10sec"! (as the activation is expected to. 
take place in a small volume around the constriction, see (§ 2), the variation 
with time will be small, possibly too small to detect even in a long time 
creep test (see below). On the work-hardening theory, it is not possible 
to determine whether v should be strain dependent or not. Inserting 
reasonable values for the parameters, v is again found to be large, and U 
is therefore insensitive to time. 

Dorn and his collaborators (1957 a, b, c) have found the apparent acti- 
vation energy U by direct measurement. These workers claim that U is 
independent of pre-strain, stress, time and creep strain. With regard to. 
the dependence of U on time and creep strain, it should be emphasized that 
such a dependence is very difficult to detect experimentally. For consider 
the case in which the maximum creep rate measurable is €,, while the 
minimum observable creep-rate is ¢€,. Let the times at which these rates 
occur be ¢, and f,, while the corresponding activation energies are U and 
U+AU. From eqn. (1) we obtain 


log (€,/€3) =AU/kT, 


where AU is the change in activation energy necessary to account for the 
change in creep rate observed. A study of the results of Dorn et al. (1957 a) 
shows that these workers used €,/é,~ 10. Thus the expected change AU 
is approximately 5k7', which, judging from the experimental error, would 
be very difficult to detect. A detectable change in U, say 12k7’, involves 
measurement of €,/é, up to 10°. Since ¢,/t, would be 10°, tests over periods 
much longer than hitherto are required. Hence the claim of Dorn e¢ al. 
that the apparent absence of a dependence of U on time or creep strain 
invalidates the present theories of logarithmic creep cannot be accepted. 
In fact it is possible that the results of Dorn et al. (1957 a) shown in table 5 
do indeed indicate a dependence of U on time, although it is not clear 
whether the trend is significant. Values of v obtained by using eqn. (4) 
are also shown; the average value of vy~10105+1° is in good agreement: 
with the predictions of the creep theories above. 

With regard to the dependence of U on applied stress or pre-strain, it Is 
clear that similar arguments apply. Whatever the dislocation distri- 
bution, the activation energy for those dislocations contributing to creep 
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at a given temperature is essentially constant over the range of time or 
creep rates so far considered. In terms of eqn. (5) for U, this means that 
the quantity »(7—7,) must be constant, independent of pre-strain or 
stress. This quantity has a very simple physical interpretation. Consider 
the dislocation loop held back by two intersecting ‘forest’ dislocations as 
shown in fig. 6. The equilibrium form of the dislocation loop is given by 


7-7, = Gb/2R = (Gb/1,) sin «9 ~ Gbao/L,. ere EY 


if x» is small, where R is the radius of curvature of the dislocations and a 
the angle indicated. Hence, for f.c.c. metals, 


u(t — Ty) = GbF2N a. eee re ey 


It follows from the above arguments that irrespective of the absolute 
values of the applied stress, internal stress, and density of ‘forest’ dis- 
locations, only those dislocations for which all these factors combine to 


Table 5. Dependence of Activation Energy on time for Polycrystalline 
Al at 78-8°K (Dorn 1957 a) 


Time to Total , 0 one 

Temp. change (sec.) | strain | (Acal/mole) S10 v(sec~+) 

600 0-18 4-56 10-2 

600 0-20 4-40 9-7 
1 500 0-21 5-02 10-9 10105=10 

3 600 0:22 4-88 10-0 

4 500 0-23 5-27 11:3 

10 800 0-24 4-72 9-5 


give a value of ag, say a,, appropriate for creep at the particular tem- 
perature, will contribute to creep. Dislocations with angles a, greater than 
x, Will effectively have taken part in the instantaneous extension, those 
with o<«, have no chance of moving in a measurable time. Inserting 
experimental values in eqn. (5) «, is found to be ~0-5 degrees for Cu at 
80°K and 6-0 degrees for Al at the same temperature. (It should be noted 
here that w, is the additional angle corresponding to the excess of the 
applied stress 7 over the internal stress z,; if the major part of T, is essen- 
tially due to local elastic interaction between the mutually intersecting 
dislocations, the ‘primary’ dislocation may be bowed out considerably to 
overcome this local elastic interaction; the additional curvature corre- 
sponding to «, will then be relatively small.) 


Or 
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It should perhaps be emphasized here that the experimental results 
discussed in §4 show that v7 and vr, separately remain constant during 
creep and are independent of stress eS pre-strain. On the formal creep 
theories only the difference of these quantities is expected to be constant ; 
the constancy of vr, is therefore an additional result which is not under- 
stood but which presumably has its origin in the geometry of the defor- 
mation (see however § 6). 

Dorn and his coworkers have also determined U as a function of tem- 
perature for Cu, Ni and Al. Their results on polycrystals are shown in 
fig. 7 together with the results obtained by Haeszner and Krause on work- 
hardened Al. In the low temperature region (below 47',,, where 7’,, is the 
melting temperature) a rough proportionality to k7’ is observed ‘at the 


Fig. 6 
(1-75) 


. (eS 


Equilibrium form of a dislocation loop under an effective stress (7—7,). 


lowest temperatures as expected from the formal theories of logarithmic 
creep with deviations from linearity near $7',,. These departures are 
possibly due to impurity, particularly in the case of Ni, as a similar but 
more pronounced effect has been observed in an aluminium alloy containing 
3-2atm% of Mg in a solid solution (Dorn ef al. 1957¢). Also, if a stress 
aided activation is considered, the activation energy will depend on the 
elastic moduli. In a change of temperature test the difference in creep 
rate observed is caused by the difference in probability of a thermal 
fluctuation and by a change in activation energy due to the temperature 
dependence of the elastic constants. Therefore the value of U derived 


from 


U € ] 9 
— =—log— /A(— Che ann eee Ca) 
z= los 4 (7) 
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should be corrected to take into account this second contribution. Basinski 
(1957) has indicated the necessary correction for Al. It is small at 80°K 
but at room temperature it substantially reduces the value of U obtained 
from eqn. (20). It is considered that the low temperature values (circa 
100°K) are little in error and from these values we have determined a 
value of v~1011+1'5 sec! (see fig. 7). 

It appears then that, contrary to the view expressed by Dorn et al. the 
two formal theories mentioned above both lead to an interpretation of all 
the existing data on low temperature creep of polycrystalline metals. It 
is not possible at present to distinguish experimentally between the two; 
the experimentally determined value of v(~ 10" sec) can be accounted 
for on either of the two theories. 


§ 6. CONCLUSIONS 


Experimental results on low temperature creep have been interpreted in 
terms of barriers due to intersecting ‘forest’ dislocations. Such an inter- 
pretation leads to reasonable dislocation widths, jog energies and stacking 
fault energies, and is consistent with the current interpretation of the 
temperature dependence of the flow stress. The hypothesis that logarith- 
mic creep is controlled by the proposed mechanism appears to be correct. 
It is also shown that the current theories of logarithmic creep account for 
all the experimental facts and that the experiments of Dorn and his co- 
workers fit in well with these theories and the proposed stress-aided 
activation mechanism. 

The interpretation used in this paper shows how dislocation widths, jog 
energies and stacking fault energies can be derived by combining measure- 
ments of logarithmic creep and flow stress. At present flow stress data are 
not sufficiently extensive to enable accurate values to be deduced except 
in a few cases. More detailed flow stress and increment tests are required, 
and it is likely that the numerical values presented in this paper may 
have to be revised. Hence the numerical estimates should be considered 
as tentative. 

From a theoretical point of view the process of intersection of dis- 
locations also requires further detailed attention. It is usually assumed 
that the internal stress 7, is due to the dislocations on the ‘primary’ 
system (see Seeger 1954). The experiments indicate that on this inter- 
pretation the density of dislocations on the ‘forest’ systems is proportional 
and at least equal if not greater than the density of dislocations on the 
‘primary’ system. It is therefore essential to consider the elastic inter- 
action between the moving dislocations and the ‘forest’ dislocations. 
This elastic interaction can be very large, for example in cases where 
Cottrell-Lomer locks are formed, and in fact larger than for dislocations 
on parallel slip planes, as intersecting dislocations can actually touch. 
The present authors would like to suggest, that contrary to the views held 
by Seeger et al. (1957), the ‘internal stress’ 7, is due mainly to the elastic 
interaction between the ‘primary’ and ‘forest’ dislocations. On this 
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interpretation the activation volume v and the ‘internal stress’ 7, depend 
both on /,, the distance between ‘forest’ dislocations, and the ‘ Cottrell— 
Stokes law’ follows. The constant , will then depend simply on the 
distribution of the ‘ forest’ dislocations. It is not unreasonable to suppose 
that the distribution of ‘forest’ dislocations is more or less random and 
similar for all metals, at least for those with the same crystal structure. 
This is the basic idea which we believe justifies the assumption made in 
this paper that «, has the same value for all metals. It should be empha- 
sized that these ideas are at present tentative; they are presented here in 
order to show how this assumption of constancy of «, might be supported. 
It is hoped to discuss this model for the flow stress of a metal in another 
paper. 

The experiments also suggest that the stacking fault energy may vary 
considerably with temperature. If true, this must be taken into account 
if the variations of mechanical properties with temperature are to be 
explained satisfactorily. It is hoped that further increment tests which 
are to be carried out in this laboratory on single and polycrystals, and 
which are to extend over wider ranges of stress, time and temperature 
than hitherto explored, will help to answer some of these and other 
questions. 
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AP PeHONGTy Ts Xa 


Introduction 


For a given pair of dislocations, D is equal to the sum of effective widths 
d, and d, of the gliding and ‘forest’ dislocations, i.e. the sum of the widths 
of each dislocation resolved in a direction perpendicular to the moving 
dislocation line. Now d,=d,°, and d,=d,°x some angular function, 
where d,° and d,° are the widths of the two dislocation ribbons. Then 
D=d,+d,=d,°+d,°x angular function, The quantities d,° and d,° 
depend on the relative orientation of the dislocations to their Burgers 
vectors, while the angular function depends on the relative orientation of 
the two dislocations. The average value (D) of D in a crystal of normal 
size will be the average value of D for all possible orientations and possible 
Burgers vector pairs. 
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Value of D for f.c.c. Metals 


If the isotropic approximation is valid the variation of dislocation width 
with orientation for dislocations lying in close packed planes is given by : 


d(ys) =d(1— K cos 2%) 


where K =2y/(2—v), % is the angle hetween the dislocation line and its 
Burgers vector, d is the width of a dislocation line inclined at 47 to its. 
Burgers vector and »v is Poisson’s ratio (Read 1953). In terms of the 
quantities defined in fig. (8) 


d,°=d(@)=d{1—K cos 2(8—,)}, 
d,°=d(¢) = dil — K cos 2(¢— po}; 


and 


d, =d(6), dy=d(d) sin 6/sin 4. 


Fig. 8 


We ™ Ald) 
; / Lok ‘aang dislocation 


/ —— ‘forest dislocation. 


Intersection of two extended dislocations. Illustration of the angles used to 
define the orientation of a ‘ forest’ dislocation relative to that of the 
‘primary ’ dislocation. 


If the ‘forest’ dislocation lies in a close packed plane 
then 


p= - x is | _ PGB) vl8, $) dp dd / ib (0, +) db db 


where the summation is over all possible Burgers vector pairs, ” 1s the 
total number of such pairs and p(6, ¢) d¢ dé is the relative probability of 
an interaction between dislocations of this relative orientation occurring. 
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Thus p(0, 4) is the ‘relative cross section’ of a ‘forest’ dislocation whose 
orientation relative to the gliding dislocation is defined by 6 and ¢. The 
‘forest’ dislocations have an average length J;, and over a volume large 
compared to 1,3 all possible orientations relative to the gliding dislocation 
may be assumed. The length of a ‘forest’ dislocation of orientation 
defined by 6 and ¢ resolved parallel to the normal of the primary slip 
plane is 
1,sin B=1,sin ¢ sin «. 

The ‘relative cross section’ of such a dislocation is/,sin ¢ sin «/l,;sin a=sin¢. 
Therefore 


ue [d{1 —K cos 2(0—6,)} sin 6 + d{1 — K cos 2(¢ — po) sin 6} ] db dé 
D= > 0/0 
Nin el eee” 
| | sin 6 dé dé 
oJ 0 


Sy od. =2d 


Nin 


Hence D,,, =2d= 2Nb. 
Value of D for h.c.p. Metals 


In metals with a hexagonal close packed structure it is assumed that the 
‘forest’ dislocations of importance in this context are those in pyramidal 
or prismatic planes with Burgers vectors for example corresponding to 
close packed directions in the basal planes of the type [1210]. Such 
‘forest’ dislocations are thought to be unextended. Therefore, 


d, =d(8) while d,=b sin 6/sin 4. 
By repeating the calculation of the preceding paragraph with 
D(8, 6) =d(8) +b sin 6/sin } 


the expression D,, ,. =d+b=(N + 1)b is obtained. 


h.e¢.p. 
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ABSTRACT 


Thermal and optical effective masses are defined and expressions for them 
obtained for metals with arbitrary Fermi surfaces. The ratio of these two 
effective masses is shown to depend on the shape of the Fermi surface. It 
follows from experimental values of the optical constants and the electronic 
heat capacity that the Fermi surface is not spherical in Cu, Ag and Au and 
that it contacts the surface of the first Brillouin zone to an appreciable degree 
in Cu and Ag. Similar evidence indicates a non-spherical Fermi surface in 
Li, Na and K. A simplification of the theory of the optical constants of 
metals is given which fits the experimental data for Na, K, Rb and Cs satis- 
factorily. Values of the plasma frequencies and widths of the plasma 
resonances of the alkali metals are calculated. 


§ 1. InTRODUCTION 


THE electronic heat capacity of a metal is related to its band structure 
and Fermi surface in a comparatively simple manner. The same is 
true for the optical constants when relaxation effects can be ignored. 
Accurate measurements of the heat capacity (Corak et al. 1955, Rayne 
1956 a, b) and of the optical constants (Schulz 1954, 1957) have recently 
been carried out for Cu, Ag and Au. Data are also available for the 
Alkali metals (Parkinson and Quarrington 1955, Roberts 1957, Wood 
and Lukens 1938, Ives and Briggs 1936, 1937a,b). In this paper, we 
attempt to ascertain just what information about the band structure or 
Fermi surface can be obtained from comparison of the results of the 
optical and thermal experiments. This aim is readily achieved for the 
noble metals. For the alkali metals, however, it is first necessary to make 
a detailed analysis of the experimental values of the optical constants 
(Ives and Briggs 1936, 1937a, b), and this, in turn, requires a suitable 
simplification of the general theory of the optical constants of metals. 
The analysis of the optical constants also leads to values of the plasma 
frequency and the width of the plasma resonance for the alkali metals. 


§ 2. THERMAL EFFECTIVE Mass 
The electronic heat capacity of a metal is y7’, where (Seitz 1940, p. 151) 
y=4n ne.) 1 1) 


{+ Communicated by the Author. Guggenheim Fellow, 1957-58, 
t Permanent Address, 
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here 
n(€y) = (2/(2n)*) [A S,|¢/AR|* eae ane (2) 


is the density of states at the Fermi energy, e,. By noting that 
v=h'de/dk is the velocity of an electron, n(e,) may be written as 

n(ep) = (2/(2r)*h)S,(l/v)yp, - . « «+ « « (8) 
where S,, is the area of the Fermi surface and 


(1/0)p= Sp |dSy(1/2). Berd ole gyn 2) 


It is customary and convenient to present experimental values of y in 
terms of an effective mass. We define as the thermal effective mass m, of 
the conduction electrons 

nim = yiy(free), 8. for ae ee fd) 
where y(free) is the value of y for free electrons of the same density as those 
in the metal. Available experimental values of the thermal masses of 
the monovalent metals are listed in table 1. The expression for m, which 
follows from eqn. (1), (3), (4) and (5) is 
NG RE SH Oa )K LIU 2 ce ee te (6) 
where S,,° is the area the Fermi surface would have if it were spherical. 
The radius of this hypothetical Fermi sphere is k,, S,,°=47k,”, and the 
volume such that it holds the correct number of electrons per atom. 


§ 3. OpricaL ErrectiveE Mass 


The dielectric constant of a cubic metal with an infinite relaxation time 
is (Seitz 1940, pp. 649-651) 
D=)P' 4D OF tert he eet a) 


w,\2 2 ai = 1 I 
(Fie fetes 2 es ee AWE 8 
s " =) i i 2 ap) ja s x =) Oo 


262 ee af il 
ee SS 3 Lerch, 
eae mo) 


§ 


where 


and 


d(w—w,,) a ee) 


for light of angular frequency w. The second term in (8) is the contri- 
bution to D’ of the acceleration of the conduction electrons by the electric 
field. The usual expression for w,” is 

w,2=4rN ,e2/m,, pO reeks (10) 
where JN, is the number of conduction electrons per unit volume and m,~" 
is the average inverse effective mass of the conduction electrons, 


Me (42 Nf?) |PREV Ze moe LT) 


The integration in (11) is over the occupied region of the conduction band. 
An expression for m, more useful for our present purposes follows from 


3F2 


764 M. H. Cohen on Optical Constants, 


(11) by partial integration, 

M, STK S| Sale > ey reer ev 
where », is the average velocity at the Fermi surface. The third term in 
(8) is the contribution to D’ of the atomic polarizability as modified in the 
metal. In both this term and in D", the sum over r includes all states of 
wave number k which are occupied, the k-integration is over the occupied 
region of the sth band, and the sum over s includes all partially or fully 
occupied bands. That is, we have «,> ey and e,<€, in both (8) and (9). 
The quantity (1/m*),., is the contribution of the rth band to the effective 
mass of an electron in the sth band, 


1 = 1 oy it 
(3), =h 2V 26,= m + > (=) . . . (13) 
where the summation is over all bands except the sth. The explicit 
expression for (1/m*),,, is 
1 Si2h ee 
(as) cauialae se ee eae 


where 
Poe Alt, Ve lene Avs Pee aoa ee 


ns” 


In eqn. (15), u,, and u,, are the periodic parts of the Bloch functions 
p,, and #,,. They are normalized to unity over a unit cell, the domain 
of integration for P,,. Finally, we note that D” arises entirely from the 
absorption of light via interband transitions which occur when w equals 
any w,,, hw,,=«,—€, Local field corrections are not included in (8) and 
(9). 

At frequencies sufficiently less than the interband absorption edge, w? 

much less than any w,,.”, eqn. (8) simplifies to 

D'=1+47Nap—(w,/w)?, . . . . . - (16) 
where 

=e 2 4 1 1 

“=F Gap 2 (as), an3 
Here N~t is the volume of the unit cell and «, is the total atomic polariz- 
ability. Local field corrections are’ not important if «, is sufficiently 
small. We see that in this limit, D’ has the simple form A—B/w? as a 
function of w. 

It is possible that the dependence of D’ on w= A— B/w? even when w 
exceeds the minimum value of w,,.. Let us consider the case where inter- 
band excitation of conduction electrons takes place and not that of core 
electrons. Define w*(k) as the minimum optical excitation energy 
for a conduction electron of wave number & and band index c. Suppose 
now, that the spread of values of [w*(k)}* is very large, its comeanrnviie 
much greater than w? and its minimum much less for the pertinent range 
of w. Let us replace the energy denominator (w?—w,,2) by w? for all 
W,~<w? and by —w,,? for all w,,>w?. This approximation is good if 
those k for which w*~q give a relatively small contribution to D’, as 


(17) 
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occurs when w* has a wide oe of values. We obtain in this way 


dea 
pia m* or Dee 


w,” wee ae Y | 
—-—< + —-—, | r(— rae 
w? w? ae k 2dr (=). j he) 


where a is the contribution of the core bands to a,. When the spread 
in values of w* is sufficiently large, the third term in eqn. (18) depends 
only slightly on w, and the fifth term depends on w nearly as w-2. The 
real part of the dielectric constant is again of the form A — B/w?, 


D’=1+4rNa, — (w/w)? . Pe ee ee LO) 
The definitions of x, and w, become apparent upon comparison of (18) and 
(19). We note further that a, is greater than the core polarizability a . 


tel Tia Ea 


Table 1. Experimental Values of Thermal and Optical Effective Masses 


Metal m,/m Mo/m 
Cu py 38a- 0-01 1-45 + 0-06 = 
Ag 0-96 + 0-01 0-97 + 0-04 = 
Au 1-16+0-02 0-98 + 0-04 = 
Li 2°32 0-98 & 
Na 1-6 (a), 1-22 (0) 1-01 + 0-02 < 
K 1-1 1-08 + 0-02 <= 
Rb _ 1-08 + 0-08 < 
Cs —_ 1-02 + 0-02 < 


(a) Parkinson and Quarrington (1955). 
(6) Roberts (1957). 

We have just shown that D’ is of the form A — B/w? if (1) w is sufficiently 
smaller than the minimum frequency at which interband absorption 
occurs or if (2) interband absorption occurs, but the absorption frequencies 
are spread over a sufficiently wide range. In these two cases, it is possible 
to define an optical effective mass my 

WN seeO) > mena es 8 teen 20) 
where Met ihOo) = 47tN 67) isg 2 a eye AL) 
is the plasma frequency of free electron at the density of the conduction 
electrons. In case (1), m, is identical to m,. In case (2), however, it is 
only possible to say that m, is less than m,. The fifth term in eqn. (18) 
has the effect of removing the contributions of those states for which 

we<@ from m, as it appears in w,? (eqn. (10)). Since these levels lie 
above the conduction band, their contribution must increase the effective 
mass, or the removal of their contribution decrease it, according to the 
sum rule (eqns. (13) and (14)). We conclude, therefore, that 


jy Silty 2 8 > (a SESS Oc (22) 


is true if the conditions for cases (1) or (2) are met. 
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Experimental values of the optical masses of the monovalent metals 
are given in table 1. The values for Cu, Ag, and Au are taken from the 
infra-red measurements of Schulz (1954, 1957) at wavelengths well beyond 
the absorption edge and hence relate to case (1), m=m,. On the other 
hand, the values for the alkali metals relate to case (2), my<m,. The 
value for Li follows from the ultra-violet-transparency results of Wood 
and Lukens (1938) after a correction for the core polarizability is made. 
The values for the remaining alkali metals are derived from a least squares 
fit of the data of Ives and Briggs (1936, 1937a, b) to the form A— b/w? 
for D’, as discussed in § 5. 


§ 4. SHAPE OF THE FERMI SURFACE 
We now show that certain gross features of the Fermi surface can be 
derived from measurements of the thermal and optical effective masses. 
We obtain 
(ny tis (Sahat) Ds LTO pao es ae eee eee 
for the ratio of thermal to average effective mass by combining eqns. (6) 
and (12). The quantities v, and (1/v),, are both proportional to integrals 
of functions of wave number over the Fermi surface. The Schwarz 
inequality applies to the integral of the product of two such functions, 
a and 6b; that is, 


| {as,ab |*< [asyar[agy, . . . . . (2) 


holds for any two functions defined on the Fermi surface. The equality 
pertains only to the case a proportional to 6. If we substitute v1? for 
a and vl? for 6 in eqn. (24), we obtain 


9 Cl lp) ct 1.) Bi a anes, eae ee 
which is merely the relation between the harmonic and arithmetic 
mean. By recalling that m, > mz, (eqn. (22)) and using eqn. (25), we see that 

M4| 1g 2 Me = (S| See) a8 ee ee) 
The seond equality in (26) holds only if v is constant. Then the Fermi 
surface is spherical, and m,=m,. 

Let us now restrict ourselves to the monovalent metals, in which the 
Fermi surface is most probably singly connected. If the Fermi surface 
does not intersect, or contact, the surface of the first Brillouin zone (BZ), 
it follows that S,,>8,,° and 

m/mM,2m/m,>1 (ocontact), . . .. . (27) 
The second equality in (27) holds only if the Fermi surface is spherical. 
Thus, the thermal mass exceeds the optical mass for a non-spherical 
Fermi surface which does not contact the faces of the BZ. 

If the Fermi surface does intersect the faces of the BZ, the intersected 
area is not a part of the Fermi surface because the states there have 
energies lower than the Fermi energy. If the area of contact is large 
enough to compensate for the departures from a sphere elsewhere in the 
zone and for the excess of v,,¢1/v),, over unity, then m,/m_, can be less than 
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unity. Certainly an experimental value of m,/m,, significantly less than 
unity must indicate contact. Our conclusions for the various cases of 
contact and no contact are summarized in table 2. 


Table 2. Interpretation of Values of m,/m, 


Experimental : 
ea Interpr etation 
m/Mq>1 Anisotropic Fermi surface and little if any 
contact. 
m,/m_21 Nearly spherical Fermi surface or anisotropy 
with appreciable contact. 
m,/m, <1 Considerable area of contact. 


To clarify the preceding statements about the effect of contact consider 
a Fermi sphere of radius R separated by h from the surface of the BZ. 
Suppose that contact with the zone boundary is provided by a cylindrical 
collar of radius 7, as in fig. 1. Ifr and A are much smaller than R, then 


A simple example of contact of the Fermi surface with a zone boundary. 


the change in R is required to keep the volume constant may be neglected 
in calculating the area. Thus S, is given by 

S,=S,) + 2arh —ar*. MS oo (28) 
We see that S,, increases until r=h and then decreases as 7 increases. 
When + equals 2h, S,, equals S,°. We designate this situation as the case 
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of critical contact. The area of contact must in general exceed the 
critical area before m, becomes less than m, in order to overcome the effects 
of departures from the sphere elsewhere and the excess of v,(1/v), over 
unity. On the other hand, eqn. (28) probably leads to an over-estimate of 
the critical radius. 

The Fermi spheres of the monovalent metals do approach close enough 
for our simplified description of the contact process to be applicable 
(h/R=0-11 for face-centred cubic and 0-13 for body-centred cubic). If we 
compare the experimental values of m, and m, for the noble metals (table 1) 
in the light of the results given in table 2, we arrive at the following con- 
clusions. The critical area of contact is appreciably exceeded in Cu. 
Hither the contact is probably about critical in Ag, or the Fermi surface is 
nearly spherical. Finally, there is little if any contact in Au but 
decided anisotropy. 

Pippard (1957) has obtained a value for the radius of contact from the 
interpretation of his anomalous skin effect measurements in the following 
way. He first found a closed Fermi surface which gave the correct 
variation of surface resistance with orientation. When the volume 
enclosed was adjusted to hold one electron per atom, the surface was found 
to project beyond the zone face. The Fermi surface was then slightly 
deformed so that it fitted into the zone. The deformation argued for by 
Pippard leads to a radius of contact 0-7 of our critical radius, which is in 
good agreement with our rough estimate. Previous arguments for con- 
tact, or near contact, in Cu were based on the ordinary transport properties 
(Klemens 1954, Jones 1955, Chambers 1956) for which the theory is in a 
rather less satisfactory state than for the optical, thermal, and anomalous 
skin effects and in which the anisotropy of the Fermi surface is beclouded 
by anisotropy of the scattering processes. 

There are two difficulties which prevent comparably detailed state- 
ments on the shapes of the Fermi surfaces of the alkali metals. First, the 
optical constants have been measured only in a wavelength range where 
interband absorption occurs so that my<m,. It may happen that m,/m, 
exceeds unity while m,/m, is less than unity, therefore, a value of m,/m, 
significantly greater than unity does not preclude either contact with the 
zone faces or a nearly spherical Fermi surface. Secondly, electronic specific 
heats have been measured only for Li, Na and K. There is, however, 
a low temperature phase transformation in Li and Na of martensitic type 
(Barrett 1955, 1956, Rayne 1954). It is conceivable that gradual trans- 
formation could lead to an effective heat capacity of the form y7'+ AT? 
and that measurements taken on heating would yield a spuriously large 
value of y. It would be necessary to obtain information on the degree of 
transformation of the specimen before the specific heat data could be 
used safely. On the other hand K, Rb and Cs do not show this trans- 
formation from body-centred cubic to close-packed structure. It would 
be of great interest, therefore, to have accurate values of y and of the 
optical constants in the infra-red for these three metals. Nevertheless, 
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if one uses the existing specific heat data unquestioningly, the indications 
are that the Fermi surfaces of Li, Na and K are anisotropic, particularly 
that of Li. 
_ It is not entirely legitimate to compare the values of m, with the values 
of my in table 1 in that the optical data were taken at room temperatures 
primarily and the thermal data all at helium temperatures. In the 
absence of suitable low temperature optical measurements, we must rely 
on band structure calculations of the dependence of effective mass on 
lattice constant for an estimate of the errors involved. The calculations 
of Brooks (1953) indicate a negligible change in the effective masses of 
the alkali metals between room and helium temperatures except possibly 
for Li. Corresponding calculations for the noble metals are unavailable. 
Our considerations have thus far been confined to the one-electron 
approximation. However, it is now known that electron correlation 
destroys the equality of m, and mg, for the electron gas (Pines 1954, 1955, 
Gell-Mann 1957, Sawada et al. 1957, Brout 1957), which certainly has a 
spherical Fermi surface. Our previous conclusion that a spherical Fermi 
surface implies equality of m, and m, is incorrect. On the other hand, 
our most important conclusion, namely that m, exceeds m, only when the 
Fermi surface contacts the zone boundary, would be strengthened if the 
effect of electron correlations were to cause m, to exceed mp, for a spherical 
Fermi surface. It is well established that the long-wavelength plasma 
frequency is given by eqn. (21) (Bohm and Pines 1953, Pines 1955, Sawada 
et al. 1957, Brout 1957) so that m, is the free electron mass m. Gell-Mann 
(1957) obtains m/m,=|1+0-083r,(—Inr,—0-203)+...] for the leading 
terms in an exact expansion of m/m, in r,, the radius (in units of the Bohr 
radius) of a sphere holding one electron. Relevant values of 7, may not 
necessarily lie within the as yet unknown radius of convergence of this 
series. Nevertheless, we can estimate the effect of the next term on the 
range of r, in which m,/m exceeds unity very roughly as follows. We 
choose f in Pine’s result for m/m, (Pine’s 1955, eqn. (8°3)) so that it agrees 
with Gell-Mann’s exact result at high densities. This requires 8 = 0-816r,1? 
and yields a range of 1-2 Sr, < 6-0 in which m,/m>1. Normal conduction 
electron densities are such that 1:8<7,<5-6. Thus, in the absence of 
more precise results, it seems reasonable that electron correlation tends 
to cause m, to exceed m,. We maintain, therefore, that a value of m, less 
than that of m, indicates contact. This view is supported by our results 
for copper. 


§ 5. OpricaAL CONSTANTS AND PLASMA FREQUENCIES 
OF THE ALKALI METALS 
5.1. Real Part of the Dielectric Constant 
Ives and Briggs have measured the complex refractive index, n,="—‘k, 
of sodium (1936), of potassium (1937a) and of rubidium and caesium (1937 b) 
asa function of wavelength A in the range 2536 A to 5780 4. They attempted 
to fit their data by neglecting the atomic polarizability and interband 
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absorption and found good agreement for Na, rough agreement for K, and 
very poor agreement for Rb and Cs. We have found that the values of 


D! =n? 12 


derived from the measurements of Ives and Briggs can be fitted quite 
accurately to the form A—B/w*. Figure 2 displays the experimental 
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Real part of the dielectric constant versus wavelength squared for Na, K 
Rb and Cs. Values of «?—n* calculated from the experimental data 
of Briggs and Ives (1936, 1937 a, b) are indicated by circles or triangles. 
The hordes lines represent the least squares fit of these values to the 
equation «?—n?=—A+ Br. The scatter appears rand cce 
possibly for Rb. : pra ae A 
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values of n?—«x? vs ? together with straight lines representing the least 
squares fit to the form — A + [B/(27c)?],® obtained by assuming no error 
in A. The agreement is excellent except perhaps for Rb, where there are 
small systematic deviations from the straight line instead of scatter about 
it. In all cases, the agreement is considerably better for the form 
A — B/w? for n?—«? than for the form 1— B/w? for —x? used by Ives and 
Briggs so that the deviation of A from unity must be regarded as 
significant. 

We interpret this agreement as indicating that eqn. (19) is valid and that 
the conduction electron excitation energies are spread over a wide range 
which includes the working range of Ives and Briggs. The values of 
47Na, obtained on this interpretation from the least-squares values of 
A-—1 are given in table 3. Values for B have already been presented in 
table 1 in the form of values for the optical mass. For both quantities the 
probable errors indicated in the tables are statistical and not experi- 
mental. 


Table 3. Values of 47Na 


Metal 47 N x(a) 47Na, (0) 
Li 0-02 — 
Na 0-13 0-10 + 0-08 
K 0-23 0-28 + 0-05 
Rb 0-28 0-50 + 0-05 
Cs 0-38 0-37 + 0-03 


(a) TKS values of « (Tessman et al. 1953). 
(6) Present analysis of the data of Ives and Briggs (1936, 1937 a, b). 


Values of 47Na, are also included in table 3, for comparison. These 
were calculated by using the Tessman, Kahn and Shockley (1953), or 
TKS, values of the polarizabilities of alkali ions in alkali halides for the 
core polarizability «». In principle, a, exceeds «) by the value of the 
conduction electron polarizability at these frequencies, which should be 
relatively small. We should expect, therefore, that the two sets of 
entries in table 3 agree quite closely except in so far as the core polariz- 
ability differs in the metal and in an ionic crystal (changed energy denomi- 
nators and matrix elements) and in so far as the local field changes the 
value of A. The two sets of entries do agree to within the statistical error 
for Na, K, and Cs, a strong confirmation of the present interpretation of 
the optical constants. In the case of Rb, on the other hand, 47No, is 
considerably larger than 47Nay. This and the poorer fit for Rb suggest 
that the conduction-electron excitation energies may not be so widely 
spread out in Rb as in the other alkali metals. 
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The value of the optical mass of Li given in table 1 is derived from the 
wavelength at which the metal becomes transparent (Wood and Lukens 
1938) that is, at which D’=0. This occurs when 


w=w2=w,2/(1+407No,), » - - + + + (29) 
according to eqn. (19), so that 
2 
_ wp?(free) 30 
ii= ETOH ibe bo ; (30) 


If we use the TKS value of a,» for o (see table 3) and the result of Wood 
and Lukens that w,=w,(free), we obtain 0-98 for the value of m/m for Li. 

Existing calculations of the excitation energies support the contention 
that the lowest excitation energy varies over a wide range in Na, K and 
Cs. One can estimate very roughly that the range of A* [=w*/c] is 780A 
to 88504 in Na (Von der Lage and Bethe 1947, Howarth and Jones 1952, 
Ham 1955) 10404 to 270004 in K (Callaway 1956) and 13004 to 16000 4 
in Cs (Callaway and Haase 1957) from energy band calculations. These 
are to be compared with the experimental range of 2536 A to 5780 4 (Ives 
and Briggs 1936, 1937a, b). Similarly one can estimate that the range of 
A* is 7004 to 41004 for Li (Glasser and Callaway 1958) which is to be 
compared with the wavelength of transparency of 1550 4 (Wood and Lukens 
1938). 


5.2. Imaginary Part of the Dielectric Constant 


The optical absorption of Na, K, Rb and Cs is shown in fig. 3. in the 
form of D" =2nx vs wavelength squared. There are two origins of non- 
zero values of 2nx, relaxation effects and interband absorption, and a 
complete theory must include the anomalous skin effect, as well. Ives 
and Briggs have already noted that relaxation effects would be far too 
small to account for the observations, even for Na. The anomalous skin 
effect is also too small. We are left with interband transitions as the 
mechanism responsible for the optical absorption of the alkali metals 
(Butcher 1951). 

The optical absorption of these metals increases systematically with 
atomic number (fig. 3). We can shed light on the origin of this by putting 
eqn. (9) for D” into the form 


477202 2 . 
Dpi= ek Sere ee Y Fe 1 
InPitae > Bay | SpePor* Pro(Yre) + - + + (31) 


In eqn. (31), we have assumed that only conduction electrons undergo 
interband transitions, dS,, is an element of the surface in k-space on 
which 

¢,(R) — €(k) =he, 
and v,, is defined as 


Vy =h|d€,/0k — 0€,/OR|. 


The area in k-space on which energy can be conserved, S,,, has an almost 
geometric significance and should not vary strongly from one metal to the 
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next. The similarity of the optical masses of the four metals suggests that 
%¢ Should also not differ widely. Therefore, the variation of P,,. with 
atomic number should be principally responsible for the differences in 
the optical absorption. 

Fig. 3 


Qnk 


O 10 20 3© 


2 9. 
A X10'in cm? 


Imaginary part of the dielectric constant versus wavelength squared for 
Na, K, Rb and Cs. Values of 2n« calculated from the experimental 
data of Briggs and Ives (1936, 1937 a, b) are indicated by circles, squares, 
or triangles. The lines are of no theoretical significance. Note the 
extreme scatter for K. The absorption is presumably associated with 
interband transitions (Butcher 1951). 


This indicated increase of P,,, with atomic number can be understood by 
representing the Bloch functions as linear combinations of ortho- 
gonalized plane waves (OPW). Each Bloch function is then a linear 
combination of plane waves and of the atomic orbitals to which the plane 
waves have been orthogonalized. One expects the matrix elements 
between the plane-wave parts to be roughly independent of wave number, 
zero, in fact, if a single OPW is an adequate representation of each Bloch 
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function. Therefore, the matrix element increases with atomic number 
primarily because of the increased number of atomic orbitals present in 
the wave function, which couple both with the plane-wave parts and 
among themselves. 

One would expect that P,, is a decreasing function of w. Indeed, P,, 
is inversely proportional to w in the nearly-free-electron approximation. 
The dependence of v,, on frequency may well be slight ; it does not depend 
on w in the nearly free electron approximation. The area S,, increases 
from zero at the absorption edge to a maximum value and then decreases 
to zero when w exceeds the maximum value of «,,.. One expects, there- 
fore that D” should increase linearly just beyond the absorption edge, 
reach a maximum value before S,,, reaches its maximum, and finally 
decrease very approximately as w~?. The data of Ives and Briggs show 
the behaviour expected on the high frequency side of the maximum. 
The rough estimates of the absorption edge given earlier are 41004 for 
Li, 8850 4 for Na, 270004 for K, and 16000 4 for Cs and are derived from 
various band-structure calculations. It would be of great interest, both 
intrinsically and in relation to the band structure calculations, to obtain 
experimental values for these cut-off wavelengths by extending the 
observations into the infra-red. 


5.3. Plasma Frequency and Line Width 


Mott (1954) has pointed out that the plasma frequency is the frequency 
of a spontaneous polarization wave in the metal (Pines 1955). A spon- 
tanous oscillation is one for which Maxwell’s displacement vector and 
therefore the dielectric constant vanish (Fréhlich and Pelzer 1955), 
For the alkali metals D is complex, and so must the plasma frequency 
be. That is, a plasma oscillation has a finite lifetime through energy loss 
via interband excitation of the conduction electrons. 

If we are to find the complex frequency 


W=WptiT OL tis “nena Mee on ee Oe) 


which is the root of 
Die Ooo seg ee eS) 


we must know the behaviour of D in the upper half of the complex frequency 
plane. We do know this for D’= A — B/w? but not for D’. Values of the 
latter could in principle be obtained by analytic continuation of D” into 
the upper half plane. However, it has proved sufficiently accurate to 
substitute for D” in eqn. (33) its value at the frequency at which D’ vanishes, 
w,A?. The resulting values of @, and IT are listed in table 4 together 
with the w (free) for comparison. Table 4 shows that w,(free) exceeds 
Wp, a consequence of the departure of my and 1+47Noa, from unity, as 
well as of the finite width, P. These results for hw, are to be compared 
with the earlier ones of Ferrell (1956) and of Pines (1956) based on a less 
complete analysis, 
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Table 4. Energies and Widths of the Plasma Resonances 


Metal hwp(free) hwp AY 
(ev) (ev) (ev) 

Li 8-0] 8-01 = 
Na 5:98 5-62 0-01 
K 4.32 3-67 0-10 
Rb 3-85 3-01 0-13 
Cs 3°42 2-91 0:30 
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A Method of Detecting an Observer Bias 


By R. M. TENNENT 
Physics Department, University of Leeds 


[Received March 23, 1958] 


In a recent experiment, electrical pulses derived from a cosmic-ray 
detector were displayed on a cathode-ray oscilloscope and photographed. 
The pulse height spectrum was found by projecting the developed images 
on to squared paper and estimating their heights by eye. In a particular 
case the pulse magnitudes were estimated to the nearest } volt and, as 
the estimations were difficult to make and fatiguing, it was considered 
necessary to analyse the measurements in order to determine whether, 
when making a particular reading, the observer was influenced by the 
previous reading. The analysis was performed by finding the distri- 
bution of runs of a particular reading, that is, the numbers of times the 
particular reading occurred once only, twice consecutively, three times 
consecutively etc., as the pulse heights were read from the film. 

Consider a set of M readings, of which N are of a particular type x (say). 
If the probability of obtaining a reading of type wz is f, the fraction of such 
readings present in the total (i.e. V/J/) is an estimate of f and approaches 
f when J is large. In this case, the expected distribution of runs (that is 
consecutive readings) of x, containing one, two, three etc. x's is given by 
the terms of the series : 


Mf(i—f)1+f+fe+fet+...). eee eee 


As an example of a distribution of this sort, the runs of even numbers 
occurring in the first 7500 digits in a table of random numbers (Fisher and 
Yates 1938) were analysed. Here N/J/=0-495 was taken as the best 
estimate of f and the observed and expected distributions are shown 
plotted in fig. 1. 

However, when the readings taken from the film referred to above were 
analysed in the same way it was found that the agreement with theory 
was poor. For example, the distribution of runs of 2 volt pulses is shown 
in fig. 2 with the theoretically expected distribution, calculated from eqn. 
(1), drawn as a full line. Clearly the agreement is poor, (a y2 test gives 
the probability of such a poor fit as less than 0-03). 

The observed distribution of 2 volt pulses can be understood if we assume 
that the probability of the observer’s calling a particular pulse a 2 volt 


= 
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pulse depends on the value of the previous pulse. In particular, with the 
distribution shown, the observer is biased towards recording a 2 volt 
pulse when the previous pulse recorded was also a 2 volt pulse. As a 
check, the same readings were analysed again for runs of 2 volt pulses, but 
taking the readings in an order different from the order of film reading. 
A distribution was obtained which agreed with that predicted by eqn. (1). 
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Number in Sequence 


Distribution of sequences of random numbers which are even. 


Suppose we now define two probabilities; f, is the probability of 
recording a reading ‘x’ after a dissimilar (‘not x’) reading and Ie is the 
probability of recording a reading ‘x’ after a similar reading. Neither of 
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these probabilities is equal to the fraction of events which is finally 
observed, and which we now call f,(=N/M). It ean be shown that f,, i 
and f, are related by the equation 


fo=aAlith—f) 


Fig. 2 


Number of sequences 


1 2 3 4 5) 6 Y 


Number in sequence 
Distribution of sequences of 2 volt pulses compared with theoretical distribution 


for an unbiased observer (full line) and a biased observer (dotted line). 


and the expected distribution of runs of x containing one, two, three etc. 
x’s is given by the terms of the series 


Mi (La fa(L fo fo°+fe° + ebe...-.), 8 ee mem 2) 
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when J is a large number. An estimate of /, may be made from the 
distribution of sequences by measuring the slope of the logarithmic plot, 
or by using the relation: 


NY Leet) aM cs 3 oh aca. 5. (8) 


where A is the total number of runs. In a similar way, an estimate of the 
probability j, may be made remembering that the expected distribution 
of runs ‘not ‘x’ events is similar to eqn. (2) above, with (1 —/,) substituted 
for fg; and that the total number of ‘a’ sequences is equal to the total 
number of ‘ not x* sequences, whence 


(een eed te et ae net ong (4) 


Using this analysis on the results displayed in fig. 2 we have M=4611; 
N=1907; fp=0-414; A=1070. We find, using relations (3) and (4), 
jf, =9-396; f,=0-439; and the expected distribution is shown as a dotted 
line on fig. 2. The fit is good, a y* test gives the probability of a worse fit 
to be 0-7. We conclude then that the distribution obtained can be 
explained by an observer bias such that the probability of recording a 
2 volt pulse is increased when the previously recorded pulse was also a 
2 volt pulse. 

Having got thus far, the only conclusion possible is that the film should 
be read again, more slowly, or with some better measuring technique. 
The method of analysing the results described has revealed an observer 
bias. However, in some experiments, results of this kind might be useful 
even though biased since the method of analysis does yield estimates of the 
two probabilities f, and f,. For instance, neon tube counters (Conversi 
and Gozzini 1955, Gardener et al. 1957) are known to be left in a super- 
sensitive state for some seconds after a discharge has occurred. If during 
an experiment such a counter were pulsed frequently, the counter would 
have two probabilities of firing, depending on whether it was in a super- 
sensitive state (following a previous discharge) or in its normal state. 
Thus an analysis of runs of discharges would show whether the effect was 
serious in the experiment and, if it was, would enable an estimation of the 
true probability of discharge to be made. 
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The Specific Heats of Gadolinium and Terbium 
between 0-2°k and 6°K 


By N. Kurti and R. 8. SarRatat 
Clarendon Laboratory, Oxford 


[Received June 4, 1958] 


THE specific heats of most rare earth metals have so far been determined 
down to 15°K (Griffel et al. 1954, 1956, Skochdopole et al. 1955, Gerstein 
et al. 1957, Jennings et al. 1957) and, in some cases, down to 1-5°K (Parkinson 
et al. 1951). Measurements at lower temperatures seemed to us desirable 
for two reasons : 

(a) They might permit a reliable separation of the electronic and spin- 
wave specific heats, 

(6) The determination of the hyperfine coupling anomalies might give 
information about the electronic states in these metals (cf. Heer and 
Erickson 1957, Arp et al. 1957) and about their use for the ferromagnetic 
nuclear orientation method (Grace et al. 1955, 1957, Khutsishvili 1955). 

Ingots of gadolinium and of terbium made from Specpure materials and 
weighing about 20 g each were used in these experiments; we are very 
grateful to Mr. A. R. Powell, F.R.S., of Johnson, Matthey and Co. for the 
loan of these samples. The metal ingot, with its heater and carbon resis- 
tance thermometer (470 ohm, 4 watt ‘Old’ Speer resistor, cf. Nicol and 
Soller 1957) was glued by means of GE 7031 varnish to a crystal of cerium 
magnesium nitrate weighing 10g. This assembly was mounted on a thin 
walled German Silver tube in a damagnetization cryostat, with manganese 
ammonium sulphate pills serving as thermal shields and getters. Cerium 
magnesium nitrate is eminently suitable for calorimetry below 1°K since 
its entropy of magnetization and cooling power is of the same order as that 
of the more usual cooling agents while its specific heat in zero field is very 
small (Daniels and Robinson 1953). 

Some difficulty was caused by the strong hysteresis heating in the metal 
when its magnetization dropped appreciably below the saturation value ; 
as a result the specimen actually warmed up as the field was reduced below 
a few kilogauss during demagnetization. In the case of Gd this difficulty 
was overcome by carrying out some measurements in a residual external 
field of 650G. For Tb, where this heating effect was marked even in 
high fields, the metal ingot had to be separated geometrically from the 
cooling salt while remaining in thermal contact with it. The calibration 
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of the thermometer below 1°, carried out in a separate experiment, was 
based on the constancy of H/T’ down to H=2kG during the isentropic 
demagnetization of cerium magnesium nitrate, 

Figure | shows the results for Gd. The curve is calculated from the 
data of Hofmann e¢ al. (1956) which, in turn are based on the analysis of the 
results of Griffel ef al. (1954). Our experimental results, which are in 
tolerable agreement with the curve above 4°x, deviate markedly at lower 
temperatures. The scatter at the higher temperatures cannot be explained 
by random errors; more experiments are needed to establish whether, as 
in other rare earth metals, thermal hysteresis is the cause. It is not likely 
that the excess specific heat is due to oxide or nitride impurities because 
they would have to be present in a concentration of about 2°%, more than 
ten times the amount one would expect. Since the entropy content of the 


Fig. 1 
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anomaly, S/R= 0-036, is less than 2% of the total spin entropy S/R=In 8, 
and is within the error given by Hofmann et al. (1956) for the magnetic 
entropy computed from experiment, the anomaly could be of magnetic 
origin. But the possibility of an austenitic transition is not excluded and 
more accurate measurements of the specific heat and of the temperature 
dependence of the saturation magnetization would be needed to decide 
this question. 

There is no anomaly due to hyperfine coupling down to 0-2°x. None 
would be expected, since the effective field of the 4f electrons at the nucleus 
should be about the same in the metal as in the salts of Gd?* and in the 
latter (Low 1956) the nuclear specific heat at 0-2°K would be C/R~10 
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only. An extension of the specific heat measurements to lower tempera- 
tures might give information about possible s-electron contribution. 

The results for Tb, given in fig 2, show below 1°K a rapid rise which is 
likely to be due to hyperfine coupling. The curve represents the sum 
of the lattice, electronic, spin-wave and nuclear contributions. For the 
lattice and electronic terms we used the values given by Jennings et al. (1957). 
The spin-wave term was determined from our results between 4°k and 7°K, | 
and the nuclear term from our results below 1°x. The curve fits the experi- 
mental data reasonably well except between 1°K-4°K where we find an 
anomaly of the same order and, possibly also of the same origin as in Gd. 
The nuclear term O7?/R = (24-8 + 1-2) x 10-% is slightly smaller than the 
value (28:2 + 0-6) x 10-3 which paramagnetic resonance measurements by 
Baker and Bleaney (1955) predict for Tb?+ in the ethyl sulphate. Since 
the electronic moments are fully aligned both in the metal and in the salt, 
one would expect the same value; the small difference could possibly be 
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The specific heat of Tb. 
————  C/Rx10?=25/T?+1-27 +0-05973 + 2-5 732, 


accounted for by a slight reduction of the mean value (r-) for the 4f 
orbits in the metal as compared with the salt. Any s-electron contri- 
bution, e.g. by conduction electrons would tend to increase the nuclear 
term. 
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REVIEWS OF BOOKS 


Proceedings of the Rehovoth Conference on Nuclear Structure 1957. Edited by 
H. J. Liexry. (Amsterdam: North-Holland Publishing Company.) 
[Pp. xvi+ 614.] 45 guilders (90s.). 

Tue value of the formal proceedings of a conference is largely determined by 

the quality of the review papers and, in this respect, the Rehovoth Conference 

was well served by a series of contributions which should answer some of the 
prayers of the non-specialist and clarify the ideas of those more familiar with 
the subjects under discussion. 

The proceedings were devoted equally to nuclear models, electromagnetic 
effects, parity violation and nuclear spectroscopy experiments and, in the 
main, they are intelligible and happily free from unnecessary detail. The 
section on parity violation, which should attract the greatest interest, contains 
discussions of the theoretical implications and summaries of the confirmatory 
experimental evidence. 

The work is enlivened by the occasional, brief, outstanding description of a 
special topic—Racah’s account of the seniority quantum number springs to 
mind—and also by the penetrating wit of the humorous Daily Bulletins of the 
Conference. 

In spite of the exceptional speed of publication, the quality of production 
maintains the high standards of the North-Holland Company. V.J.E. 


Surface Chemistry. Theoryand Applications. By J.J. BrkeRMAN. (Academic 
Press Inc.) [Pp.x+501.] $15.00 


THIS is a second edition of the author’s book “ Surface Chemistry for Industrial 
Research ” published ten years ago. The author expresses the hope that it 
will be used for the teaching of surface chemistry in universities as he feels 
that a unified account of this wide subject is long overdue. The reviewer 
has found in the book clear and simple accounts of branches of this subject 
which are perhaps little known to most physicists, such as the principles of 
foam formation and electro-capillarity. As a book of reference and an intro- 
duction to what is known in this wide field of physical chemistry, the book will 
be very useful indeed. On some subjects, however, it does not go very far 
or give up to date references. The discussion of surface phenomena and 
crystal growth, for example, makes no mention of such concepts as dislocations, 
kink sites and so on, and the section on friction does not mention some of the 
most recent work. None the less it contains a clearly set out body of informa- 
tion which is not easily available anywhere else. ; SNM: 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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Internal Slip zones, low mag. 11 x 100. 
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Notch effect, 20x 1000. 
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Linking of fissures at fracture stage, tapering section, c. 40 x 500. 
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Peak effect,” 20 x 500. 
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Combined notch and peak effects, 20 x 1000. 
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Slip and fissures parallel to surface, 20 x 1000. 
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Crossed slip and fissures, 20 x 1000. 
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Block movement beginning, 20 x 500. 
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Block movement enhanced, 20 x 500. 
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Block movement with disintegration, 20x 500. 


Fig. 13 


Ww. A. WOOD Phil; Mag. Ser. 8, Vol. 37 Pla37e 


Fig. 14 


Fig. 15 


Crazy slip, 25 x 500. 
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Irregular cracking at large amplitude, 25 x 500. 
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Ship-less cracking at large amplitude, 25 x 500. 


Fig. 18 
ee re 
Sa ‘ ~ 
Sd 
ae 
f 4 
ie OE, 
ee 
nee 
Bee pk a eS 
eg 
hi GOO we 
ad e 
all ’ bes el 
ft mage lta ial Ye * F 
Bee 
° 
ie + b- =. ie ania ape 
* * 
: 7 Mijn? ws pine, wat! Mp 5° # i een hos 
rv - im a tae ee alll KS 
ie * 
* x 
* 
i, OP 
€ a ied % 
ego ’ = x 
Er espa OM, « * Si ing al 


‘ Balling-up’ of slip zones after heating, x 1000. 


Fig. 19 


Cavity formation by heating, 20 x 1000. 
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Fig. 20 


Before and after bending of section, 15500. Slip zones with notch at head 
open into cracks (which extend into silver plating above). 


Fig. 21 


Before and after bending of section, 15 x 500. Slip zones with notch at head 
open into cracks (which extend into silver plating above). 
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Fig. 1 


Slip lines originating from ‘rosettes’ in a lightly strained magnesium oxide 
crystal (x 50). 
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Fig. 3 
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Three secondary cracks lying within a kink band. One tip of each crack is 


situated along the kink band boundary A-B (polarized illumination 
x 300). 
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Fig. 4 


(d) 
75 


hotomicrographs taken at different levels: beneath the crystal 
te The ee crack in (g) is included in fig. 3 (top, right). 
Line A-B represents the boundary of the kink band. (e) and (f) were 
photographed on internal surfaces cleaved after compression.  (¢) and 

(g) were printed on reversal of the film (polarized illumination x 400). 
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Fig. 6 


Effect of etching on a kinked region. (a) before etching (polarized illumina- 
tion) : (b) after etching (normal illumination x 300). 


